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Abstract 

, Fffects of rotation on onset of convection during plane-front directional solidification of Ph <?n 
and the pseudobinary system mercury cadmium telluride fH^ PH tm 

1 Pb p S " h r ^"stuCie/by means' of “ °" de " dmiC 

, corporattng Conolis and centrifugal accelerations into the momentum eauation of Cnriell *» 
of we find that under realistic processing conditions, a large degree of X5irario?£n £ 
ac teve using modest rotation rates for both Pb-Sn and mercury cadmium telluride (He Cd Tet 

phsn, 

51 “ 

F ? r H 8i-xCd x Te, we have also studied the nonrotating case. The key differences are due to 
nraai^ 

r? T^u Se Where den . s,,y de P ends monotonically on temperature and cot^osirion for Ha 

rotation irfMbhluhema^con'vrctionkjnTsingl^^ binarv 

tastabilUy 0 ^ ^“o^cy-driven convection, and has 

At large growth velocities, the plane-front interface between liquid and solid becomes 
Tone F f S H CC d t0 3 instability and solidification occurs dendritically with a mushv 

rotot- f dei l dntes a n d interdendritlc fluid separating the solid from the melt. For the’ Pb-Sn svsten? 
rotation substantially suppresses the onset of convection in the mushy zone and in the overlvina 

defects.^ d ' ng ° Pen thC Pr ° miSe th3t r ° t3tion Can su PP ress heckling and other macrosegregation 
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1. Introduction 

SESS^iS^^SSSaS 

steady and one-dimensional plane-from case inStablhties can cause depanures from the nominally 

- mXuS?» ^ Ch 

aa- 1 — “■ stfa^^.’Bssftasss 

composition. When an alloy ^soHdified* hv* coolhin "f d . e P ends .on both temperature and 

gradient overcomes the stabihrinTteLeratu IfmadiVn, ?L!? e , COndlll ° n ^ thls adverse 

aut^re^we^thLm^onmt^w^ri/H^^^tOr ^rpholo^ r ^l^^ S ^^'^”^^ eC ^^ 

that the buoyancy force does not senShlv aiSr^kl^ ™^ h i g cal or convective instabilities, and 
which occurs at higher wavenumbers than does rh^h 116 ™ 0 " ^ ? nset . of morphological instability, 
was reviewed by Glicksman et S ( QjSf IS? &?T y ' d 2^ Subsequent work 

ssrvssss 

(1981) studied the efferrc off im//fi suggesled rotation might reduce freckling. Ridder eial 
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Sned *** °" ? acrose 8 re Sation in directly 

solidification in a low-gravity environment^ ^1™^ ^f^nQQOU 1 nSdUCed ^ 

"o^ a ^vi; nderC ° nditi0nS WhCTe Ihe ™ 8 niLe d „S y 1^XrS 

oppo^d^lo^n'tacfforces suc^s'p^ssijre^vfscours^re^^and^ 6 ^^' 0 ^ 2C ‘ °" the ««““• As 

surfaces of a fluid element body force? £t n„ ?>! 1 ’ and surface tension, which act on the 

strengths are proportional to the local fluid density Tn a . flui< ^ e . le ment. Hence, their local 
forces, there are other "pseudo-body forces" which manifest < themf r f Vltatl0 r a ^- a - nd magnetic body 
(centripetal and Coriolis) when the reference frame to whfoh rh^f^v! VCS * S fi ? dtl ° us accelerations 
steady rotation relative to an inertial Sme ( a™hh f ^ md motlon 15 referred undergoes 

will not be concerned, manifests itself if the rotation? 2, ^ fic . ut mus ^deration, with which we 
frame is unsteady.) AZughS fi^me relative to the inertial 

case of the gravitational acceleration^ thev h™t 1 ° l corTes P° nd dir ectly to forces (as in the 
associated with body forces and can have^vnamnV?! SamC mathematlcaI form as accelerations 
gravitational and magnetic forces discussed abov^Gree^fa? ioS S ? Q r al ! y a l P rofound as the 
that modification of the graSonal field SS T 1968) * In , hght of this ’ ™ d the fact 
always be possible, the poSh^ ™ gnetic ™ d may not 

undergoing directional solidification is of interest PP h of convectl °n in a liquid 

The e^mf.^ d «“«*' by Schulz-DuBois (1972). 

and Weber et al. (1990) have shnJn Lt rll; S p and H ^ llaweI1 0984), Muller (1990), 

macrosegregation in binary aNoys d^aion^ly^^difi^MnHer^i 1 ^ 108 / 11 ^ red ““ ."!« degree of 
Kou and Kou et al. studied the effect of srpaHv^mtaH CC !j? nder Plane-front or dendntic conditions. 

Sn-Pb alloys. Sampt and Hellawd Sdere? SS * 7T* °" freckle Nation in 

studies S 'show* 1 | a h^t ^he^egree^ macro^gi^ation^fr^kT 

reduced by rotarion. macrosegregation (freckles and srnations) can be significantly 

fluidMnTems^a^eference^rame^ladnTwbh 0 ab ° U ' 3 flXed , axis can influen “ motion of a 

measured from the axis of rotation u ic thp lrv^i fi.„H Ration Here, r is the position vector 
frame, and p L is the local X“of the HquM d Vd ° C " ) ' rc,a " Ve '° Ihe noninertial «&■*«<* 

convwt^n^av^been^estn^cled^ to'ca^s 6 ^^^^^'.!’? °" ‘ hC ° nSe * of buoyancy-driven 
unbounded layer of a single-comDonent finiH rh^n s °bd | fi c ation occurs For a horizontally 
(1955), and Niiler ^and^ Bifshonn K h . ' Chandrasekhar (1961), Chandrasekhar and Elben 
axis 4 sig„i fl clm,; d i„ B h‘S h ,°^ 

buoyanVdnVen^Salem^ ^.rjos^S^^f ^ r i' a '' ( ' 9M > h ave compS 

maimained between teeXf t^JJtaW b ^' a ‘ gradient 







ssg^giggg?^ 

g^s^a^SSSSSS 

sssifeslpsps 

buoyancy-driven convection in a roStingfluid 6 P ° tentiaI ° f ^ C ° ri0liS acce,eration ^ suppress 
convective and nrorpholo^^^ the C ° riolis acceleration on 

s£8S 

shown in the 1992 University of Illinois Ph D Di w™* “ft 8 , 18 * a " d Pl ots showing results are 

T>a^Sm2mn 80 iC |992") gini d^" g> ’ ‘ n 3 pape . r co-author^khDr. Oaeki n^mUu^lcS 

2. Effect of Rotation on Plane-front Solidification of Pb-Sn 

of Coriell e a ml'nTuoaTcumfwere mneraT^ h" d c T pared our results 10 those 

mmmmMm 

individually computes .he ^pSa^ « which 

pres^^ou^prinf^afres^l^fn terns oTstal^^ybou^dariMT^he^V^'^pi^ie^^ ^' ^v^'* ^ 

growth rate and cjs the bulk concentration of tin. ^e critical value of r“ fc W J?7 u l ’? the 

»»«^-js5S£rSSS^S 





c£ff?f Cted il y rotation >. whereas as Q 0 increases, the value of V, at which the onset of instabilitv 

max, mum re afv= stabilization by rotation occurs near the locafminimum!“nd that for (l - 

Coriefut at SiSh “"“" lr “ 1 ° n of Sn is increased more than a hundredfold. Although 
conen et at. (1980) noted a local minimum in the stability boundary near V = 1 ,, c-r-i in ,iL 

nonrotating case for the largest gravitational acceleration considered, they offered no explanation 

Ld tam Gr -"of: |? C “ v/lTm 1 ™ c ° nsec i“ en « of the fact that as v(-> 0, the concentration 
u W~V/(k PL0 D) vanishes. As the temperature gradient is independent of r anH 
is stabilizing, the critical value of Coo must ultimately increase as V, 4 0 (Of Se the weS 

^“1° S ,°'T C - cann0 ‘ cxceed 100 >- As cl must initiallyVecmaS wh^LUasina V 

there must be a local minimum on the convective branch before c* can increase to ir»in nn .u 1 ’ 
morphological branch of the stability boundary. This nonmonot^i c Sa tf "• on V can 
be interpreted in terms of the existence of fouf critical values of V, fOT?ertairvalues"of c '™" 
1990) hat f ° r suffklently lar S e V 1 . c - ultimately increases on the morphological branch'fDavrs 

in the^ridcTwa^ue^fc^afv* 0 "! 0613 ' 113 !^ 6 a, J ow ? ow,h ra,es < more tha " “ hundredfold increase 
in me cnucai value ot at V, = 1 4 seer 1 can be achieved by rotating the laver at 500 mm far th * 

Tay lor-l^udrm^ 

parallel to the axis of rotation in a uniformly rotating inviscid fluid is prohibited at anv nonzero 
ofShidnn r C ° US n .'f one sees . an inhibition of the onset of steady convection, with chedecree 

stead eC fjo n i^^^y^^^^^diKTl^h^facfthar^^aylt^ProuSni^tl^rem^i^Iie^o^rto 

The ^^ologic^'^stabiH^^TOcure^t'^ve^^t^wavelen^dis^so 5 die^i^tion^s^aTnmst 
perpendicular to the solidification front (i.e., aligned with the axis of rotation 1 Henee th*» r 1 ■ r 
accejmtion does not sensibly affect the morphological instability. ‘ ’ th ConoIls 

the Ph ^ re cfc. Vely m ? de l t rotatlon rates required to significantly inhibit the onset of convection in 
t e Pb-Sn system make the proposed method an interesting candidate for a program of laboratory 
experiments. Experiments might be conducted using ! completely fSU^^n^caTSSd 

tnnitahl/ S**? °",/ he v k ertlcal and mounted axisymmetricaHy on a routing horizontal 

turntable. If the liquid at the top of the rotating mold is in contact with a gas or vacuum it will 
have a nearly_parabo oidal free surface, on which the elevation above ^ thfpdm o^th? axis of 
? tatl i° n i S AZ ~ Although this configuration is consistent with rigid-bodv rotation it 

also leads to a nearly paraboloidal solid-melt interface. For high lotttiol? rateso? Ce mold 

tih?rnH!?f S h W1 i n !f d t0 si s nificant radi al variations in the solidified alloy. Also cooling at 

r^ntlnftlaiwf 0 ^ 17 i eadS t0 3 fadial tem P erature gradient, which in turn leads to significant 
entnfugal effects. Experiments of the type proposed above were initiated some time ago bv 

Smion PnVate COmmumcanon) for the crystallization of ammonium chloride from aqueous 
nqut SSgf SUPPreSS b - y - y -«ven conv^ticn in the meSd SUSSES 






3 ' Pla™e-Fronf^SoHd^ Pb-^n^ e< * UCetl Gravi,Jr ° n S,abili( > of 

solidification were studied for dS^S? h S,ab i li,y of Piane-front 
suppressed to a far greater degree than woSld taexu^ed on buoyan p- d ™n convection is 
suppression by the individual effects of reduced b ,,n < £ earher w . ork . in which 

(Oztekin and Pearlstein 1992) was demonstrated ^ ^ ° ne et a ' 9 ^) anc i uniform rotation 

conv^SoTta ‘.SSlKS? undeZil [h s :S‘ S a °'T"°" 7 «"** on the onset of 

effects of reduced gravity have been® considered by SieM «°^‘ "ft f ““a Th 'j"<“vidual 

combinations of the liquid-side temoerature ornHi Jn^ 11 J ' v , ^ who showed for several 
stability boundary in the V r c' plan7conS?f m^h^o." d f'T level 0 < < 1 'hat the 

the convective branch is stabilized by reducing tfi e af an! 7 , c h orlvec,,v f branches, and that 
convective branch increases as C decreased ’j t , \ V t the cntic al value of c„ on the 

gravity) with 0 L = 200 3 and fl, !' 0 100 2«> 3^ aS7rl? 92) 7^ for « = 1 <«■" 
rotation stabilizes the convective branch (i e for ariv V Ipm , (their FigUre 1} that 

branch increases with Q 0 ). For example with V yahie of c„ on the convective 

which plane-front solidification is unstable is more than ’ a critical Sn concentration above 
500 rpm than in the nonrotating case ThevTso SS ° f magnitude higher at Q 0 = 

unaffected by rotation, and that ? C WMo&al branch w°as 

shifts from the convective branch to the mnmW n • V , a u Ue °[ X 1 at whictl the onset of instability 
the nonrotating case to about 27 p sec-’ for decreases from about 40 p sec~i in 

Kh 6 .^each^conJs^o^ some^ofo^irffo^moro^ f °- r ^ = - 1 ’ 10 ~ 2 ' 1(H > and 

joined to a convective branch. The critical value on C ~ dec T eas ! n & with increasing V,) 
decreasing £, showing the inhiWto^^^S^^^^? 11 ^ branch ^creases witii 
convection. Note that for V = l n 2^-1 • red “ ced gravity on the onset of buoyancy-driven 

Sn concentration at which planeTront solidifiSion Scom l °JM 8e increase s the critical 
magnitude relative to the full gravity case We also ™te7h»fT ab l ®.f y . about four orders of 
gravity is relatively less at lower solidification n °! e u - th u t the i ?tabihzmg effect of reduced 

the morphological instability branch is not affected bv so ld,fiGation rat e. As expected, 

value of V. at which the onset 7insS wbereas as C decreases, the 

decreases irom about 40 p sec- 1 for C = 1 to about L3 f ^ °n5 6 mor P hoIogical branch 

For the same liquid-side temperature gradient fGi — 200 if rm -h r 
rotation rates (Q 0 = 5 rpm) the de?rep of ctihir ~ 20 ?. K cm *). w e find that for small 

considerably larger than th£ sum of ^individual efrecfs Xo^tionf If *5 Ca ? gravity level is 
full gravity case, the critical Sn concentration is increas^ m^re aM d re ^ U ?^ gravit >'- As in the 
higher solidification rates. For C = l(h 2 and O, = s mm ^ Wer solld tfication rates than at 

about an order of magnitude for V, - 0 5 ll sec-i C cnt |c a ^ I concen^tion is increased by 

V, =0.5 p sec-i is very smlll fo contrast to rhe’ the L nh ! bition of convection fdr 

convection by reduced y gravity is larger at hilh! tatlI ! g case » the inhibition of the onset of 
solidification rare a! wE V *„'„*£ ^'in ^ fab^i^SoT' HenCe “! e decrease »f the 

25 rpm^he erbietd by^the rotation 

™ r = e r rders of magntode 

l(H and =*o!and C / T1 - i IH ndVl 0 - n s 7 m 7 ^ V an<l = 0. t = 1 and fio = 5 ipm, 5 - 
™ effectivc,y inh?bi,s ' he 

in the critfcal ™?ue of can te acEd ^ 3 hu " dredfold increase 

^ ^ undoubted l^due 5 
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4. Stability of Plane-Front Solidification of H gl . x Cd x Te 

on temperature density does not depend monotonically 

stability of plane-front solidification when the i n ran S® of composition, we consider the 
normal case where 7 c density derend^^^^^^ !? n C ° oIed from ***>"• contrast to the 

Pb-Sn alloys considered by Coriell et al (198011 for „^P erature , ar) d composition (e.g., the 
parameters (solidification rate nominal limHH-dH? ve^? c ® rtain combinations of the operating 
concentration) there can exist a critical value of the t ^ rn P erat . ure gradient, and bulk solute 

front solidification is unstab e at SdimenSnieL , moIe fractlon (CJ Ww which plane- 
plane-front solidification at any C Is™ eJ X sSble STf^ ^ wher ?, as norm ^ case 
Moreover, when the density varies nnnl^Z • n sufflc,entl y small solidification rates, 
value of the dimensionles?solidfficatioT™?i v cf eh with temperature, there can exist a critical 
unstable for all C In Ms a£e Y ‘ ? Iane - front solidification is 

ssss 

t^Kel^ 

soIu^on^^HgdT^ cons^^bl^imponance d 

system in the fabrication of electro-optic detectors. dUC '° he Wlde use of thls pseudobinary 

liquid-side temperature gradient r^G|V/[ L." (gD^/f P tt!!? sahT^ of the dimens ionless 

SS? r f ‘'.Ww branch 

the other hand, below he aSS™ ’’"‘"'i' f °l SOrae range of a - 0n 

all wavenumbers decay an7 pTa „° u c °" vec "y' branch, disturbances of 
rate Te = 2.70 x 10-* (V> 102 u/°Z “ stable. Beyond the critical soUdification 

intersect, there is no stable range of bulk mole fraction CFmrt wu m< iiv3 0l i! 8 ' cal , branches 
in the y-C_ plane is analogous to that for r= 8 2 x itH Th. m»j, * Ihe stable region 

the right and the convective branch has shifted eliahrlv Ho.. 6 rno ri a holog , cal branch has shifted to 
higher y (y = 5 61 x 1(V4 v — ? 19 ,./ c ~ a a u ® y downward. Their intersection occurs at a 

leL ,o further enlargement of, he sabfc nt '^ S Jf, bIe regi °? is largCT - r still mom 

6.56 x. 0-3 ( g l . Va ‘ UeS ° f ^ *= fOT F “ 3 - 28 x 'O' 3 a " d 

characteristic of the'nomtafcase ir^whkth densttv varies r UrVC ^ and stability boundaries and those 
(e.g., the Pb-Sn resultsrfCorid ", , , oggSjTS "2* le ™ pera,ure and imposition 

or alters the mechanism by wWch .hfonLt i mio^^^^^ “i 110 ' 10 " of stale af fi« 

binary liquid cooled from below. ° D occurs dunn g directional solidification of a 

net density s trati fica'tio^^n * ^ he ^and^is^n r ^ or P holo § ical instability is insensitive to the 

composition outside ^relativel'y ilXfTh- 0 ' ' " m P ara, “ t re a " d 

comparison of the neutral curves and srahHitv hi„„H ° e iater * ace - This is evident from a 

(Coriell « at. 1980; OzSin S pSdn JsMh “ 0r ' he CaSe 10 ,hose for Pb 'Sn 

the da P e ^ s ,? tron fy »n the nature of 

configuration in which there exists adjacent to thl interface a subfaye?!rwWchTe posi“ve 
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^s j&aEm aBgaBS 

sp|p^slSi=i 

in thJ? C C f ST FoVH^'rvV^T 5 ha ^ e pr0f ?^ conse< i uencei for the stability boundaries 

3KSSK if ’ S °' U,e bU ' k m °‘ e fraCti ° n " SUffi ^ n,ly ' OW < C °™'' « "?980T6t!cin 

fractp &-&£. p ski?e asssssf s zzszzsg&z ° f aa: 

hnnl a ^ Ue , S Of C ” Plane ' fr0m vilification is unstable with res^c to tKe onsIi of o^staS 

™ *t ?' "•■ f ° r ^'ch plane-front solidification is stable in a region bounded above bJ the 
morphological and convective portions of the stability boundaries y the 

Qtr a t;^H Sin8 -^ ^-component fluid with an unstably stratified layer overlying or underline a stahlv 
tratified region, convective motion is typically localized in the unstably stratified lave/ hut ma v 

th^ « ed - that convectlon occurs in the unstably stratified region and penetrates into 

^*.*0985 ss5a^^sss!s, , i^ 

rhpc/ci 11 ?! 1 - 0n ^, vertlcal f em Pf rature gradient and a nonuniform vertical solute gradient Both of 

wWphrtJ C M?°r the d , 1StUrbanC t S t0 be Iocalized ab °ut the neighborhood of tlfe critical depth at 
ich the salinity gradient reaches its minimum value. In our work, when the liquid- P solid 





disturbances associate? wit^S^osdHato^Lse? temperature gradient is destabilizing, the 
confined to that thermally unstably stratified sublayer? buoyancy-driven convection are largely 

of the distuK^^ and the localization 

gradient is destabilizing. Hence the results 2 H l f Subl ? yer ln which the temperature 

v^thin^ I ayerT 0th er b * nar y or multicompmienU^itids^ni^cnfdj^^ 

dramatic qualfiativ? effea? o^^he^^^illsta Whv^n an^h!? T? temperature can have 
P^ ane “front solidification The existence nf un ^ unc * e d horizontal fluid layer 
stratification is destabilizing should also hav? imDOrtn?i?v? a sublaye J in which the thermal 
occurs in vertical Bridgman growth in amnnniLT^J r ”1 J^ quences for the convection which 
other processes used to grow Hg, Cd Te crystals fwf "lif radl } IS ’ as weI1 as in zone melting and 
have no. been observed in sim?aS?fS TodaK ' however, these eftos 

geometries, which have used equations of sEZ, convection tn Hg,.,Cd,Te for these 

ejwndence of density on temperature. The axisymmetric smndario acc ? um f or 1 the nonmonotonic 
by Ktm and Brown (1989) employ a linear equatonTf mte of »f ventcal Bridgman growth 
mole fraction in the liquid near the interface is about n ffit ;° f Il,e . form (*)• even though the CdTe 
known to depend nonmonotonically on temperantre for StS ^ ““Pnhttion and the density is 
computations of Apanovich and Ljumkis ( 1991) ^nons up to at least 0.1. The 

of 0J3. ln wh'ch density of Hg^,Cd,Te is^nde^endem ^emp^ture ^a^^mole^action 

5 ' Hgr x Cd° x f Te 0,a,i0n °" S ' abili,y ° f P,a "- Fr »n< Solidification of 

venicafS SSJSt S <*««* -*»> ™don about die 

cadmium telluride cooled from below. Since for suffidentI^m?? na ? S °i ldif ? Cation of mercur y 
does not depend monotonically on temDeratnre rh^ lcient y sn l a11 x the density of Hg, Cd Te 
mole fraction, nominal solidification rate and HnnidVvi 3 nln ® e of operating conditions (bufk C^Te 
» the liquid-solid interface them fra “bla^e Twfi'mTe™^ f?!- c,ll) . for *«■ advent 
This differs from the normal case fi e hinarv ,ii m " C r e thermal stratification is destabilizing 
such as Pb-Sn), in which th^the^ sM.fSiZTT’u? 1 !* a Mnear equation of sal; 
Hg, Cd^Te melts there can exist a critical hulk rviT* 00 i S ?^ ab ? llzin g everywhere. Thus, in 
sohditicanon is unstable at all solidification rates V bel ° W which P lan e-front 

solidification at any x„ is linearly stable for sufficiently small V S ill the non ? aI case plane-front 
can exist a critical solidification rate V. cr such that for V ?v” r Mo J eover ’ 1 for H ^i r x Cd x Te there 
for all x Our results show that modest rates of uniform ro P ane '^" on t sohdification is unstable 
onset of buoyancy-driven convection for th e f H? DiT’c" sl S ni J cant Jy suppress the 
substantially increase the critical solidification rate l^ S L en l’ and that rotati on can 

unstable for all bulk CdTe mole fractions 1 & y d Wh,Ch P ,ane - f ™t solidification is 

° S ? 1 !rS ry conv ^ ctive branc h . C For G?= To^K em'^th? J h P.[? OI ^J ol °g ic : al branch joined to an 
and 100 rpm show that above the i J stability boundaries for Q =0 25 50 

wavenumber a and plane-front sol id i f i ?aU o^' u n^mbl e Ch ^ f 31 s . t I ur f >an J Ces grow for some ’range of 
appropriate to each rotation rate shown, disturbances ero? ^ ° W the convective branch 
between the morphological branch and the convective s ? me / an | e ^ a - However, 

cl hi f °Ki X °° x ~ < x ~^> disturbances decay for all wavenumber? a*?! 3t f d W l tb eacb ro . ta don rate 
stable. Note that the morphological branch is unaffected hv mt ^ anc * plane-front solidification is 

convective branch is shifted downward The ^riricaf LTrdV^f inc, S ases - Ihe 

convective braneh is a dcGr^sing function of Q at anv cnli Aifin m ° e ^f actlon x « on the 

tnhtbttory effect of rotation on the onset of 
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there is no stable nu^ofbuik irote < ft^ri<^T feC T^S^i mQt ^ol°gH^- bran 5 :faes intersect, 
at all bulk CdTe mole fractions for sufficiently high snlidifinari lane-frontso! hdification is unstable 
rate V,“ at which the motpholoS and MnjSe taSfS n ‘“ S ' ^ Cri,ica ? solidification 
increasing rotation rate Q , The critical value V CT • S 1 i lt 1 f rsect occurs at higher x with 

n o = 100 tpm relative to Vhe nonrZl ca“ We Z^thX T* ,ha " a faC,0r of a > 

coi^^t^stobl'yat'tii'e lov^t vahre of tT ^ 

reducing tfee range of solidification rat^ 7nd’ bulk™™^? glC ?' 1S Shifted t0 the ~ Ieft - 

solidification is stable, as expected on the basis of rhf* tT°if ^?,n 10 - n ?i for whlch plane-front 
Without rotation. On the other hand ?h. Ln“,?5^h '“S* s ° f ?° rle " e,aL (1980) for Pb-Sn 

G l =50 K cm-' case, comes^nding manTncrease inZrllTlf 1 , do , wnWardS - relati ™ 10 «* 
We note that decreasing the temperamre gradient has vervKn?iff£f a S e °P erat ' n S conditions, 
a nonrotating layer, but reduces the range of stable hulk orr^ e ^ e . on .the onset of convection in 
for a layer rotating at 100 tpm The Sal vSue of v fractions by more than twofold 

G 0 = 100 tpm relative to the nonrotaS TaM almond Lf" by , mor ®,^ an . a f ?«or of five at 
factor of ten predicted for G r = 25 K cm -1 Althrmah^h^ degree of stabilization is less than the 
the critical value of V. at whfch the mSnholo] S h hc rcducuon ,n .V has very little effect on 
the nonrotating case, tor Q= 100 rpm toe critical ySmaTV**™ Stablh £ y boundaries intersect in 
the temperature gradient is reduced ^ decreases by a factor of three when 

of moiphologlctd a * cach ™"“. °f o, onset 

convection sets in via disturbances with relativelt fnnaf? ( g , Renumber) instability, while 
this point is discussed in - 1 " the foll owing section 

convection. tt1e mechan ism by which rotation inhibits the onset of 

linear equation ofstate with* rejectiono mtito ?* ^ t COoI f d [f°™ below and characterized by a 
the vertical temperature and hquid ' s ° Iid *"««*** (e.g., Pb-Sn), 
throughout the liquid layer. In the Pb-Sn svstem The nnci IZ fif a ° d des j ablll zing, respectively, 
via monotonically growing disturbance*: ac chA u r> ° • buoyancy-driven convection occurs 
pseudobinary Eg, CdT svstem .hfruxt r b * Condl “ al (1980 )- However, for the 

the liquid-soIid 8 interface there can exist a ^ublavw^r'T'il^ S ' a ? l,dng ’ Md adjacent to 
destabilizing. In this system the onset of mnv^h ayer 10 ^lch the dermal stratification is 

that the system is unstable at 'the indicated values of G^fo" afi^ wh 2 ' 00 ' disturbances - w = note 
corresponding the absence of a stab, y stratified *.%& .^cZ 

the onsefrf buoyX W «ven c2nvifon‘^c'’ura d , Pea ^ S,ein I992 > thal significant suppression of 
known Taylor-lW^^SS^ASr ,t • °", ra,eS ' and is due » > h ' well- 
rotation since it occurs at very short wavelengths' JiX£?^S* a tnstabihty is unaffected by 
solidification front (i.e„ aligned with the rotation axis). * * nearl ly Perpendicular to the 

6 ’ Pb^Sn ° f Ro,a,ion on s,abili ‘y of Dendritic Solidification of 

unde^oing^oHdlfi^atioi^by^TOlmg^om Sw^The °wo s “ lu,ion f° r Pb ‘Sn melts 

porous medium with anisotropic Dermeahilifv w / WO P base mu shy zone is modeled as a 
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to be dynamical varifbles. The ^ ^ ,iquid ’ “* taka " 
self-consistent nonlinear model of solidification ex srt’ ff?r nni d USmg a ^^odyuamically 
operating parameters (solidification rat^S^S 



solutions corresponding to' plane front'and^SeS!?' y pace l ^ c -2 >niin S 10 ‘he existence and stability of 
the Pb-Sn system ifnde r somt SSotlhe dSS.v 3 *?' *? prese "" d a " d discussed to 
solidification is linearly stable with resoect to mnmhninff ^ solution exists when plane-front 

in light of previous work on nonlinear morphological stabilftv^For a^ 5 ' ,T hi ? P°. int is discussed 
gradient, the dendritic solution is unstable belnv/a r ® g liquid-side temperature 

composition. As the imposed temped solidification rate, nearly independent of 

decreases, and the dendritic solution is stabilized decreases - the critical solidification rate 

* Coriel, , al . 

Sn in Pb used corresponds to an infinitely difute solution^ Th^™ 6 ^ P ° mt ' dlffusivit y of 
as a porous medium with anisotropic oemieabititv Thci ^ two-phase mushy zone is modeled 
the boundaries separating the SfromT S' he ,OCal P orosit y> as weI1 as *e locations of 
taken to be dynamical vanablesThe e a " d l he L mush y zone from liquid, are 

thermodynamically self-consistent nonlinear mn^ nS f° na L t S,C state, computed using a 
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temperature gradient, and bulk solute concentration} A HfSci 3tC f ? ominal liquid-side vertical 
to the existence and stability of solutions corresnonHin^ T S OI A° f ^ P^ ameter s P ace according 
presented and discussed for the Pb Sn sv^etT f nlr P lane - fen < ™ d dendritic solidification is 
exists when plane-front solidification il lS? i JJu **1°™ condltlons dendritic solution 
This point is discussed in light of previous S res P ect to morphological disturbances, 

liquid-side temperature gradient the dendritic snintin - neaf m< ^P^ogical stability. For a given 
rate, nearly IhclZnlTT^ bdow * c , ritical ^d^ation 

cridcal soHdi&adon rare decrease and ,he denc^c'Xto SS " * 
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variation of H m with y is much weaker and for r = 0 28S (C - ^ ^ C - m * 

independent of y We also note that the miichv. 7 nno v j 200 K cm -1 ) is nearly 

temperature gradient. ^ 6 ^ lc ^ ness decreases strongly with increasing 

onse^f buoyancy-dri ven^onvectirin ^we^ hav^to^fl curves jV a) correspond to steady 

inability in fhe dendritic soScS °f 

solidification rate (denoted by /) is the maximum poL on ea g ch^ SSh cul^ 
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ided wording to the stability of the dendritic solution. In region 6, the dendritic and ph£e- 



front solutions are both stable, as for T = 7 12 x 10-2 while in reoi™ n *u a j • • , . . 

wwtll 1 — I 

prcdic,ed b * 

mmm wmmm 

S«S£k«pS5S 

ouTS" gradlen ‘ ' S COnS ' S,em Wi,h the locario " and character of U,e ckStol'^o^l and a^ 

Ra e 'l T‘ ma Ifri£ 

ayleigh number R m we define for the mushy zone is proportional to AC m H m (Worster 1992), 


n sa 




"ear 0.2. Fowler (1985) defined a 

J« permeability, and sh^edTat fe^^is^ntkafvalue 1 R a ^l/?C t0n c a w 'Ti e h aSin f g D f ' ,n S don of 
dendrilic solidification is stable with respect to ™d t l: B ~ CJ J uch ‘, hat lf R < R j. then 
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of th?Kq3S velocity assume a maximum on the liquid side 
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stable, it is not surprising that we also find a one-dfmensk>^ ,S hnearl y 



References 


Alexander J.I.D. D.J. Wollkind, and R.F. Sekerka 1986 J. Crystal Growth, 79, 849-865. 

exander J.I.D. J. Ouazzani, and F. Rosenberger 1989 J. Crystal Growth, 97, 285-302. 
Antoranz, J.C. and M.G. Velarde 1978 Phys. Lett. A, 65, 377 - 379 . 

Antoranz, J.C. and M.G. Velarde 1979 Phys. Fluids, 22, 1038-1043. 

Baines, P. G. and A.E. Gill 1969 J. Fluid Mech. 37, 289-306. 

Bhattachaijee, J.K. 1988a Phys. Fluids, 31, 2456-2461. 

Bhattacharjee, J.K. 1988b Phys. Fluids, 31, 2462-2466. 

Bhattacharjee, J.K. 1988c Phys. Rev. A, 37, 1368-1370. 

Biihler, K. and H. Oertel 1982 J. Fluid Mech. 114, 261-282. 

Chandrasekhar, S. 1953 Proc. Roy. Soc. London A, 217, 306-327 
Chandrasekhar, S. and D.D. Elbert 1 955 Proc . Roy. Soc. London A, 231, 198-210. 
Chanrtasekhar, S. 1961 Hydrodynamic and Hydromagnetk Clarendon, Oxford, pp. 76- 

Conell, S.R., M.R. Cordes, W.J. Boertinger, and R.F. Sekerka 19807. Crystal Growth. 49, 13- 


Davis, S.H. 1990/. Fluid Mech. 212, 241-262. 

Fowler, A.C. 1985 IMA J. Appl. Maths. 35, 159-174. 

Fultz, D. and Y. Nakagawa 1955 Proc. Roy. Soc. London A, 231, 21 1-225. 

Giamei, A.F. and B.H. Kear 1970 Metall. Trans. 1 , 2185-2192. 

Ghcksman, M.E., S.R. Cornell, and G.B. McFadden 1986 Ann. Rev. Fluid Mech. 18, 307-335 
Goroff, I.R. 1960 Proc. Roy. Soc. London A, 254, 537-541. 

89L J ^ ’’ $ p Nandapurkar, and D.R. Poirier 1989a Metall. Trans., 20B, 883- 

Heinnch, J.C., S. Felicelli, P. Nandapurkar, and D.R. Poirier 1989b AIAA Paper 89-0626 
Homsy, G.M. and J.L. Hudson 1969 /. Fluid Mech. 35, 33-52. 

Homsy, G.M. and J.L. Hudson 1971a/. Fluid Mech. 45, 353-373. 

Homsy, G.M. and J.L. Hudson 1971b /. Fluid Mech. 48, 605-624. 

Homsy, G.M. and J.L. Hudson 1971c Int. J. Heat Mass Transfer 14, 1 149-1 159 
Homsy, G.M. and J.L. Hudson 1972 Appl. Sci. Res. 26, 53-67. 

Huppert, H.E. 1990 /. Fluid Mech. 212, 209-240. 

Kim, D. H. and R.A. Brown 1989 /. Crystal Growth, 96, 609-627 

KOU Ins2 "‘ t£K3T" “ EleC!r ° S,aS Reme " ed '"^ Ph D ' ^sis, Massachusetts 


Kou, S., D.R. Poirier, and M.C. Flemings 1978 Metall. Trans. 9B, 71 1-719. 
Maples, A.L. and D.R. Poirier 1984 Metall. Trans. 15B, 163-172. 





Masuda, A. 1978 7. Oceanogr. Soc. Japan, 34, 8-16. 

Muller, G. 1990 7. Crystal Growth, 99, 1242-1257. 

Muller, G. 1988 Crystals, 12, Crystal Growth from the Melt, H.C. Freyhardt, ed., Springer- 
Verlag, Berlin. 

Mullins, W.W. and R.F. Sekerka 1964 7. Appl. Phys., 35, 444-451. 

Nakagawa, Y. and P. Frenzen 1955 Tellus, 7, 1-21. 

Niiler, P.P. and F.E. Bisshopp 1965 7. Fluid Mech. 22, 753-761. 

Oztekin, A. & Pearlstein, A. J. 1992 Metall. Trans. B, 23, 73-80. 

Pearlstein, A.J. 1981 7. Fluid Mech. 103, 389-412. 

Poirier, D.R., M.C. Flemings, R. Mehrabian, and H.J. Klein 1981 in Advances in Metal 

Processing, J.J. Burke, R. Mehrabian, and V. Weiss, eds., Plenum, New York, pp. 277- 

Polezhaev, V.I. 1988 in Crystals, 10, Growth and Defect Structures, H.C. Freyhardt, ed., 
Springer-Verlag, Berlin, 87-141. 

Riahi, N. 1983 7. Phys. Soc. Japan, 52, 2620-2621. 

Ridder, S.D., S. Kou, and R. Mehrabian 1981 Metall. Trans. 12B, 435-447. 

Sample, A.K. and A. Hellawell 1984 Metall. Trans. 15A, 2163-2173. 

Schmitt, R.W. and R.B. Lambert 1979 7. Fluid Mech. 90, 449-463. 

Schulz-DuBois, E.O 1972 7. Crystal Growth, 12, 81-87. 

Sekerka, R.F. and S.R. Coriell 1987 in Proc. Sixth European Symp. on Material Sciences under 
Microgravity Conditions, Bordeaux, December 2-5, 1986, ESA-SP-256, 3-11. 

Sengupta, S. and A.S. Gupta 1971 Z. Angew. Math. Phys. 22, 906-914. 

Stem, M. E. 1960 Tellus, 12, 172-175. 

Trivedi, R. 1984 Metall. Trans. A, 15, 977-982. 

Ungar, L.H., M.J. Bennett, and R.A. Brown 1984 Phys. Rev. B, 31, 5923-5930. 

Ungar, L.H. and R.A. Brown 1984 Phys. Rev. B, 29, 1367-1380. 

Venugopalan, D. and J.S. Kirkaldy 1984 Acta. Metall. 32, 893-906. 

Veronis, G. 1963 Astrophys. J. 137, 641-663. 

Wollkind, D.J. and L.A. Segel 1970 Phil. Trans. Roy. Soc. London, A 268, 351-380. 

Worster, M.G. 1991 7. Fluid Mech. 224, 335-359. 

Worthem,S., E. Mollo-Christensen, and F. Ostapoff 1983 7. Fluid Mech. 133, 297-319. 
Walton, I. C. 1982 7. Fluid Mech. 125, 123-135. 

Zangrando, F. and L.A. Bertram 1985 7. Fluid Mech. 151, 55-79. 




Coriolis Effects on the Stability of Plane- 
Front Solidification of Dilute Pb-Sn Binary Alloys 

ALPARSLAN OZTEKIN and ARNE J. PEARLSTEIN 

The possibility of using steady uniform rotation about a vertical axis to suppress the onset of 
buoyancy-driven convection during solidification of a binary alloy is considered using a linear 
stability analysis. For Pb-Sn alloys, our results clearly show that the onset of convection in a 
horizontally unbounded layer can be suppressed significantly at modest rotation rates. Specif- 
ically, plane-front solidification is linearly stable at higher Sn concentrations in a rotating 
configuration than in a nonrotating one. The predicted inhibitory effects of rotation on convec- 
tion are discussed in terms of previous experimental and theoretical studies of the effect of 
rotation on the onset of buoyancy-driven convection in single-component fluids heated from 
below and in binary fluids subject to thermal and solutal stratification. 


I. INTRODUCTION 

During directional solidification of alloys, it is fre- 
quently desired to produce large single crystals with very 
low densities of macrosegregation defects and other im- 
perfections. In principle, this can sometimes be achieved 
by “plane-front” solidification, in which the melt-solid 
interface remains perfectly planar. In such a case, the 
solidification process would be steady in a reference frame 
moving with the interface, and the only spatial variation 
would be in a direction normal to the interface. How- 
ever, in real systems, several instabilities can cause de- 
partures from the nominally steady and one-dimensional 
plane-front case. 

First, the solid-liquid interface may be subject to the 
so-called morphological instability, which has been stud- 
ied extensively since the early work of Mullins and 
Sekerka. m This instability causes deformation of the 
nominally planar interface, ultimately leading to for- 
mation of a two-phase “mushy zone” of dendrites and 
interdendritic liquid. Departures from one-dimensionality 
and steadiness in the mushy zone result in nonuniform 
distribution of solute in the solidified material. 

Second, the density of a binary or multicomponent melt 
depends on both temperature and composition. When an 
alloy is solidified by cooling from below, rejection of 
solute(s) at the growing interface is potentially destabi- 
lizing if the solute-enriched liquid just above the inter- 
face is less dense than the warmer overlying bulk liquid. 
Under some conditions, this adverse solute gradient 
overcomes the stabilizing temperature gradient, leading 
to convection in the melt. This fluid motion provides 
another transport mechanism, besides molecular diffu- 
sion, for redistributing solute(s) into the bulk liquid from 
the relatively enriched region near the interface. Con- 
vection in the melt is often referred to as thermosolutal 
convection or, because the diffusivities of heat and sol- 
ute are different, as doubly diffusive convection. 
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Convective and morphological instabilities in a binary 
alloy undergoing directional solidification were First 
studied by Coriell et aL {2] using a linear stability anal- 
ysis. These authors showed that motion may occur due 
to either morphological or convective instabilities and 
that the buoyancy force does not sensibly alter the cri- 
terion for onset of morphological instability, which oc- 
curs at higher wavenumbers than does the buoyancy-driven 
instability. Subsequent work was reviewed by Glicksman 
et alP ] and Sekerka and Coriell. [4] More general discus- 
sions of the effects of convection on plane-front and den- 
dritic solidification have been given recently by Worster, (5] 
Davis, [6] Huppert, [7] Polezhaev, m and Muller . 191 

Buoyancy-driven convection in the melt has been shown 
to be the dominant factor in the formation of “freckles,” 
a macrosegregation defect deleterious to the mechanical 
properties of directionally solidified alloys. The forma- 
tion and characterization of freckles in nickel-based 
superalloys were first studied experimentally by Giamei 
and Kear . f101 Poirier et alJ ] ^ investigated macro- 
segregation in electroslag ingots, showed that convection 
in the melt results in freckling in the solidified material, 
and suggested that rotation might reduce freckling. Ridder 
et al. U2] studied the effects of fluid flow on macro- 
segregation in nominally axisymmetric ingots and showed 
that melt convection results in macrosegregation in the 
mushy zone. In a theoretical study of a binary alloy so- 
lidifying radially inward, Maples and Poirieri 131 con- 
cluded that macrosegregation results from natural 
convection in the mushy zone driven by nonuniform 
temperature and solute distributions. A recent review of 
the effects of buoyancy-driven convection on macro- 
segregation in binary and pseudobinary nonmetallic sys- 
tems have been given by Muller . [9] 

One means by which the onset of convection can be 
inhibited in a density-stratified fluid layer is to subject 
the layer to a magnetic field aligned parallel or perpen- 
dicular to the stratification . 191 For horizontal Pb-Sn layers 
solidified from below at several growth velocities, Coriell 
et aL {2] showed theoretically that the critical Sn concen- 
tration above which instability occurs can be increased 
an order of magnitude by applying a vertical magnetic 
field of the order of one tesla. This technique requires 
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that the liquid be an electrical conductor and so is ap- 
plicable to metallic alloys, semiconductors, and aqueous 
solutions. 

Macrosegregation might also be controlled by increas- 
ing or decreasing the magnitude of the gravitational ac- 
celeration or changing its direction. Alexander et a/. [14] 
and Heinrich et al . [15 ,l6] theoretically studied the effect 
of reduced gravity on macrosegregation in directionally 
solidified alloys. These authors suggested that macro- 
segregation in alloys can be reduced by solidification in 
a low-gravity environment. Muller n 71 and Weber et al. [m 
have recently discussed solidification under conditions 
where the magnitude of the body force is greater or less 
than that of normal gravity. 

Both gravity and an external magnetic field are body 
forces which act on the liquid. As opposed to contact 
forces such as pressure, viscous stress, and surface ten- 
sion, which act on the surfaces of a fluid element, body 
forces act on the mass of a fluid element. Hence, their 
local strengths are proportional to the local fluid density. 
In addition to gravitational and magnetic body forces, 
there are other “pseudo-body forces” which manifest 
themselves as fictitious accelerations (centripetal and 
Coriolis) when the reference frame to which the fluid 
motion is referred undergoes steady rotation relative to 
an inertial frame. (An additional fictitious acceleration, 
with which we will not be concerned, manifests itself if 
the rotation of the noninertial frame relative to the in- 
ertial frame is unsteady.) Although these accelerations 
do not correspond directly to forces (as in the case of 
gravitational acceleration), they have the same mathe- 
matical form as accelerations associated with body forces 
and can have dynamical consequences equally as pro- 
found as the gravitational and magnetic forces discussed 
above. 1191 In light of this, and the fact that modification 
of the gravitational field or imposition of an external 
magnetic field may not always be possible, the possi- 
bility of using rotation to suppress the onset of convec- 
tion in a liquid undergoing directional solidification is of 
interest. 

Several effects of rotation on solidification have been 
discussed by Schulz-DuBois.* 201 The experiments of 
Kou,* 211 Kou et al .,* 221 Sample and Hellawell, [231 Muller,** 71 
and Weber et a/. 1181 have shown that rotation can sig- 
nificantly reduce the degree of macrosegregation in bi- 
nary alloys directionally solidified under plane-front or 
dendritic conditions. Kou and Kou et al. studied the ef- 
fect of steady rotation about a vertical axis on freckle 
formation in Sn-Pb alloys. Sample and Hellawell con- 
sidered solidification of die transparent alloy NH 4 C1-H 2 0 
in a crucible rotating about an axis inclined between 0 
and 30 deg with respect to the vertical. Muller* 171 and 
Weber et al. have studied the effect of rotation on the 
growth of Te-doped InSb crystals and have shown that 
striations can be suppressed at sufficiently high rotation 
rates. These studies show that the degree of macro- 
segregation (freckles and striations) can be significantly 
reduced by rotation. 

There are two means by which steady rotation about 
a fixed axis can influence the motion of a fluid. In terms 
of a reference frame rotating with constant angular ve- 
locity (l about an axis, these correspond to the centri- 
petal and Coriolis accelerations, for which the terms 


p L i 1 x H x r and 2 p L il x u, respectively, are added 
to the momentum equation. Here, r is the position vector 
measured from the axis of rotation, u is the local fluid 
velocity relative to the noninertial reference frame, and 
p L is the local density of the liquid. 

Previous theoretical studies of the effects of rotation 
on the onset of buoyancy-driven convection have been 
restricted to cases where no solidification occurs. For a 
horizontally unbounded layer of a single-component fluid, 
Chandrasekhar,* 241 Chandrasekhar and Elbert,* 251 and Niiler 
and Bisshopp 1261 have shown that steady uniform rotation 
about a vertical axis can significantly inhibit the onset 
of convection, with the Coriolis-related Taylor-Proudman 
mechanism* 271 playing the dominant role. The effects of 
centripetal, Coriolis, and gravitational accelerations on 
convection in horizontally confined rotating fluids in cy- 
lindrical containers of various aspect ratios were consid- 
ered in an early series of articles by Homsy and 
Hudson.* 28-321 More recently, for a single-component fluid, 
Weber et a/.* 181 have computed buoyancy-driven flows 
equivalent to those which can be driven by a temperature 
gradient maintained between the ends of an otherwise 
insulated right circular cylinder rotating at constant an- 
gular velocity about an axis perpendicular to and inter- 
secting the cylinder axis but not passing through the 
cylinder. Three-dimensional computations, in which the 
variation of the magnitude of the centripetal acceleration 
ft x ft x r along the cylinder axis was neglected (a 
good approximation when the cylinder length is small 
compared to the shortest distance between the axis of 
rotation and the cylinder), were performed with the term 
accounting for the Coriolis acceleration 2 fl x u either 
included or omitted. Weber et al. found excellent agree- 
ment between experiment and computation when the 
Coriolis acceleration was included. 

For a binary fluid, Pearlstein* 331 has shown that the 
Coriolis acceleration can either stabilize or destabilize a 
horizontally unbounded layer, depending on the values 
of the Prandtl and Schmidt numbers, the dimensionless 
rotation rate (expressed in terms of a Taylor number), 
and the dimensionless temperature or solute gradient 
(expressed in terms of thermal and solutal Rayleigh 
numbers). Other work concerning the effect of rotation 
on doubly diffusive convection in binary fluids (with no 
phase change) has been reported by Sengupta and 
Gupta,* 341 Masuda,* 351 Antoranz and Velarde,* 36 371 Schmitt 
and Lambert,* 381 Riahi,* 391 Worthem et a/.,* 401 and 
Bhattacharjee.* 41 ’ 42 431 

In this work, we use linear stability analysis to study 
the effect of the Coriolis acceleration on convective and 
morphological instability for alloys which solidify with 
a nominally planar interface. This work serves to iden- 
tify the mechanism by which Coriolis effects affect the 
onset of convection in solidifying binary alloys. The 
analysis is illustrated by results for the Pb-Sn system. 

II. BASIC STATE AND 
LINEAR DISTURBANCE EQUATIONS 

We adopt the model of solidification used by Coriell 
et a/.,* 21 in which the Oberbeck-Boussinesq equations 
govern motion in the liquid. In a reference frame trans- 
lating with the nominally steady velocity (V,) of the 
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moving interface and rotating with constant angular ve- 
locity ft, the basic state 


0 = (0,0, -eV,) [1] 


^ ^ , ™/Tx t p L0 K L G L 

* i. — * m H : i 

k Ps^i 


1 - exp I — 


PsVj 

Plo k l 


c l = c . 


f , I " * { PsV ! 

1 + exp I 

^ V PlqDl 


m { c * k s G s 
T s — T m + — 1- 


K L 


. , V' 

1 - exp I ; 


[ 2 ] 

[3] 

[4] 


is the same used by Coriell et a/., where u, f L , c L , and 
T s are, respectively, basic state values of velocity, tem- 
perature, and solute concentration in the liquid and tem- 
perature in the solid, k l and k s are the thermal 
diffusivities of the liquid and solid, respectively, D L is 
the diffusion coefficient, p s is the solid density, p L0 
is the liquid density of pure lead at its melting point, c 0 c 
is the bulk concentration, G L is the liquid-side temper- 
ature gradient at the planar interface, V, is the nominal 
growth rate, k is the segregation coefficient, m L is the 
slope of the liquidus, T M is the melting temperature of 
Pb, e = ps/ptu - I is the fractional shrinkage, G s = 
(LV I + Gik,)/k s is the solid-side temperature gradient 
at the planar interface, k r and k s are the thermal con- 
ductivities of the liquid and solid, respectively, and L is 
the latent heat of fusion per unit volume. (All concen- 
trations used herein are weight percentages.) We note 
that the basic-state temperature distribution in the solid 
(Eq. [4]) is valid only near the interface. 

As in stability analyses of other flows subject to uni- 
form rotation. M vv - 441 the linear disturbance equations 

we use differ from those in the nonrotating case only by 
addition of the Coriolis acceleration to the disturbance 
momentum equation. (This approach neglects the con- 
sequences of density variation in the term corresponding 
to the centripetal acceleration.) The equations governing 
small disturbances in the liquid are 


V ■ u = 0 

du p s du Vp 

— - v, — + 2H + U = ot T gT L i z 

ot Pli I dZ PiQ 


[5] 


— a c8^jz + vV 2 U 


[ 6 ] 


bt l 

Ps 

<>T, ^ ( 

PsV, 

dt 

V, 

1- G l w exp 


Plo 

dz \ 

Plo k l 

9 c l 


dc ‘ ^ ( 

PsV, 

■ 

Vi 

r O c H’ exp 


dt 

Plo 

dz F V 

PlqDl 


[7] 


2 I = D,V 2 c 


L v U 

[ 8 ] 


where the dependent variables u, w, p, T L , and c L are, 
respectively, the disturbance values of the velocity, its 
z-component, pressure, temperature, and concentration 


in the liquid. Here, i 2 is the unit vector in the 2 -direction, 
g is the magnitude of the gravitational acceleration, v is 
the viscosity, a T and a c are the thermal and solutal ex- 
pansion coefficients, respectively, G c = ( k - 1 )p s c*V,/ 
(kp L0 D L ) is the concentration gradient at the planar inter- 
face, ft = ftol. is the angular velocity, and ft 0 is a con- 
stant. The disturbance energy equation in the solid is 


dT s BT S 

— ~V,— 
dt dz 


tfcV 2 Tc 


[9] 


where T s is the disturbance temperature in the solid. 

Taking the curl of Eq. [6] and the dot product of the 
result with i., we obtain an equation 


dco z p s dto. 

— - — V,— 
dt p L0 dz 


„ dw 

— 2ft 0 — = vV'ai 
dz 


[ 10 ] 


for the 2-component of the vorticity, to.. After twice tak- 
ing the curl of Eq. [6], the vertical component of the 
resulting equation is 


3 , p s d dto, 

-V 2 w V, — V 2 w + 2ft 0 — - 

dt p L0 dz dz 

= a rgS\T L + a c gV\c L + uW [11] 

The disturbance boundary conditions at the interface are 

dr)(x,y,t) 


u{x,y,0,t) = eV, 

v(x,y,0,t) = eV, 

w(x,y,0,t ) = -e 
dw(x,y,0,t) 


dx 

dy 

dypc^t) 

dt 


dz 


= -eV,V 2 ± y(x ,y,t) 


[12a] 
[12b] 
[12c] 
f 1 2d] 


r dy(x,y,t) , dT L (x,y,0,t) f dT s (x,y,0,t) 
l — ft i I- ks 


dt 


dz 


dz 


( Ps G L V,k L G s V,k s \ 

[ ]v(x,y,t) 

\Plo Kl k s j 


[12e] 


DlPl o 


V (x,y,t) + V,c L (x,y,0,t ) + 


cudyO^yj) 
k dt 


_ dc dx,y,Q,t) 

Ps(l _ k) dz 

T L {x,y, 0,t) + G L y(x,y,t) = T s (x,y, 0, t) 


fl2f] 


+ G s y(x,y,t) 


{ 1 2g] 

T L {x,y,Q, t) + G L rj(x,y,t) = m L G c y(x,y,t) 

+ m L c L (x,y,0,t) -I- T M VVl V (x,y,t) [12h] 

where rj is the interface position and is the capillary 
coefficient. 
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III. NONDIMENSIONALIZATION 
AND FOURIER DECOMPOSITION 

We scale the velocity, time, length, temperature, con- 
centration, and vorticity with k l /H, H 2 /v, H , G L H, c x , 
and k l /H 2 , respectively, where H — D L /V { is the char- 
acteristic length (associated with the destabilizing gra- 
dient of rejected solute, see Eq. [3]), and write the 
horizontal and temporal dependence of the nondimen- 
sional disturbance quantities as exp (or 4- ia x X 4- ia v Y), 
where a x and a s are the x- and y-components of the 
wavevector and a is the temporal eigenvalue. Substitu- 
tion into Eqs. [7] through f 1 1 ] yields 

<7 (D 2 - a 2 )W -—D{D 2 - ar)W = -Ra T a 2 6 L 
Sc 

- Ra s Lea 2 Ci + (D 2 - a 2 ) 2 W - Q*Dfl z [13] 

crSl z - — Dn z = n^DW + (D 2 - a 2 ) fl z [14] 
Sc 

Pr(70 t - pLeD& L = (D 2 - a 2 )® L - Wexp (-pLeZ) 

[15] 

yp 

S caC L - pDC L = (D 2 - a 2 )C L + — W exp (-pZ) 

Le 

[16] 

in the melt and 

PraQs ~ LcD&s = *(D 2 ~ a 2 )©* [17] 

in the solid. Here, Pr = v/k l is the Prandtl number, 
Sc = v/D l is the Schmidt number, Le = D l /k l is the 
Lewis number, Ra r = a T gG L H 4 /(k l v) and Ra s — 
a c gc?,H */(D l v) are the thermal and solutal Rayleigh 
numbers, respectively, f 1* = 2 Yl 0 H 2 /v is the square root 
of the Taylor number, p = Ps/Plo I s the density ratio, 
k = k s /k l is the thermal diffusivity ratio, and y = 
(1 — k)/k. The boundary conditions (Eqs. [12a] through 


[12h]) at the interface become 

W(0) = -sPr a(3 [18a] 

DW(0) = sLQa 2 p [18b] 

ftz(O) = 0 [18c] 

Scafi = -k*D& L ( 0) + k* s DQ s ( 0) 

+ Leikfp ~ k*/K - 1 /k)P [18d] 
(1 - G)P = 0 5 (O) — 0 L (O) ri8e] 

1 

(Sccr + kp)l 3 = -kC L ( 0) DC l ( 0) [18f] 

YP 

n m* 

(- <GT * + m* - l)/3 = — C L { 0) + 0,(0) [18g] 

PY 


where W , C u O z , and ® L are, respectively, the ampli- 
tudes of the disturbances to the vertical velocity, con- 
centration, vertical vorticity, and temperature in the melt, 
0 5 is the amplitude of the disturbance temperature in the 
solid, and ft (a constant) is the amplitude of the disturbed 
interface position. Here, we define dimensionless pa- 
rameters k* = G L k L H'/ (LD L ) , k* = G L k s H/(LD L ), 


m* - m L G c /G L , T% = T M V/(G L H 2 ), and G = G S /G L . 
For the far-field boundary conditions, we follow Coriell 
et al. l2] and set all disturbances to zero 

W = DW = n z = C L = ® t = 0 as Z— ► sc [19a] 
©5 = 0 asZ^-oc [19b] 

far from the interface. 

IV. NUMERICAL SOLUTION 

Our objective is to find conditions under which infin- 
itesimally small disturbances neither grow nor decay for 
a finite number of wavenumbers and decay for all other 
wavenumbers. Disturbances which neither grow nor decay 
are said to be neutral. The neutral disturbances can be 
of two types, depending on the imaginary part of cr. If 
the imaginary part of a is zero for a neutral disturbance, 
the onset of instability will be via monotonically grow- 
ing disturbances (steady onset). If the imaginary part of 
cr is not zero, the neutral disturbance will oscillate in 
time (oscillatory onset). 

In the previous section, we formulated an eigenvalue 
problem for two systems of ordinary differential equa- 
tions on two semi-infinite intervals (in the melt and in 
the solid), coupled by boundary conditions at the de- 
formable interface. For convenience, we follow Coriell 
et al. [2] and solve the problem on a finite interval [—h^h], 
where 2h is the dimensionless height of the computa- 
tional domain. With regard to the far-field boundary 
conditions at Z = we set all disturbances to zero 

W = DW = il z = C L = 0 L - 0 at Z = h [20a] 
in the melt and 

0 5 = O at Z=-h [20b] 

in the solid. In this work, we have used h = 10. (For 
steady onset, we have checked a number of our results 
using more accurate asymptotic boundary conditions ap- 
plied at Z = ±h derived following Keller's 1451 procedure 
and have found excellent agreement between the eigen- 
values computed using the two sets of boundary condi- 
tions.) Since we use Chebyshev polynomials in our 
numerical solution, we scale the liquid and solid regions 
by z 2 = (2 Z - h)/h and z t = (2Z 4- h)/h, respectively, 
so that each region lies between - 1 and + 1 . The re- 
sulting system is then solved using a spectral Galerkin 
technique developed by Zebib. [46] The problem thus is 
reduced to a matrix eigenvalue problem 

Ax + crBx = 0 [21] 

where a is the temporal eigenvalue and the elements of 
the square matrices A and B depend on a 2 = a\ + a 2 
(the square of the horizontal wavenumber), the bulk con- 
centration Co c, and the other dimensionless parameters. 

In what follows, we characterize the stability of the 
nominally plane-front solution in terms of the growth ve- 
locity V } and bulk concentration c*, with all other pa- 
rameters taken as fixed. For each value of V h we seek 
one or more critical values of c* (denoted by c£), such 
that for Coc < c£, disturbances decay for all wave- 
numbers, while for c« > c*, disturbances grow for all 
wavenumbers in some range. In order to determine the 
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critical value(s) of c x , we First determine the neutral curve 
(c* vs the wavenumber a), separating those combina- 
tions of c x and a for which all temporal eigenvalues cr 
lie in the left half-plane (a stable basic state) from those 
for which at least one eigenvalue lies in the right half- 
plane (an unstable basic state). To determine the bulk 
concentration on the neutral curve for an arbitrarily cho- 
sen wavenumber a , we first guess a value of c x and com- 
pute all eigenvalues cr using Eq. [21]. If all eigenvalues 
have negative real parts, the value of c x is increased by 
doubling the previous value; otherwise, the new value 
of Cx is chosen as half the previous value. This process 
is continued until we determine two values of c* be- 
tween which at least one sign change is obtained in the 
real part of the least stable temporal eigenvalue. The 
concentration on the neutral curve is then determined using 
a bisection method. 

To compute the critical concentration c* at which in- 
stability first occurs ( i.e the minimum on the neutral 
curve), we arbitrarily choose a wavenumber and com- 
pute the corresponding c a on the neutral curve using the 
procedure described above. We then fix c x at the value 
computed at the previous step and compute the eigen- 
values a for a discrete set of wavenumbers in a chosen 
range. From this set, we select the wavenumber corre- 
sponding to the cr with largest real part. If this wave- 
number is at an endpoint of the chosen range, we extend 
the range to include the wavenumber corresponding to 
the largest Re(cr). We then select this wavenumber and 
determine the corresponding c x on the neutral curve. We 
continue this process until the relative change in c x is 
less than 10 -6 . 

V. RESULTS 

The solutions of Eqs. [13] through [19] depend on 15 
nondimensional parameters, as defined in Sections II and 
III. To determine the stability of a basic state with a 
nominally planar interface, numerical values of these pa- 
rameters need to be specified. (Solutions on a finite in- 
terval depend also on h.) Therefore, it is not possible to 
numerically explore the effects of more than a few com- 
binations of these parameters on the onset of instability. 

In this work, we have thus restricted ourselves to the 
Pb-Sn system which, due to the low melting points of 
both components, has been the subject of several ex- 
perimental studies. The parameter values (other than O*) 
are as used by Coriell et al . 121 at reference conditions 
corresponding to pure lead at its melting point. (The dif- 
fusivity of Sn in Pb corresponds to an infinitely dilute 
solution.) 

To test our code, we first considered the nonrotating 
case and compared our results to those of Coriell et al. 
Taking the liquid-side temperature gradient at the inter- 
face as G l = 200 K cm 1 , we computed neutral curves 
Cx vs a for various values of the solidification rate V,. 
The neutral curves were generally similar to those shown 
by Coriell et al. The only qualitative difference was that 
in our work, several new oscillatory neutral curves were 
found to branch from steady neutral curves found by 
Coriell et al. and in the present work. We believe that 
our detection of these additional oscillatory neutral curves 
(which in each case lie well above the critical value of 


Cx and are hence of no practical consequence) is due to 
use of a numerical technique which simultaneously com- 
putes a large number of temporal eigenvalues at each 
combination of c, and V,, as opposed to the shooting 
technique of Coriell et al. which individually computes 
the temporal eigenvalues by a one-point iteration scheme. 

We present our principal results in terms of stability 
boundaries in the V, - d plane. For ease of comparison 
to the work of Coriell et al., ll] results are presented in 
terms of dimensional variables. With G L fixed at 200 K 
cm ', Figure 1 shows stability boundaries for fl 0 = 0, 
100, 200, 300, and 500 rpm. For each value of fl 0 , the 
stability boundary consists of some portion of the mor- 
phological branch (c* decreasing with increasing V,) found 
by Coriell et al. (their Figure 1) joined to a convective 
branch. The critical value of c*. on the convective branch 
is an increasing function of fi 0 at any value of V,, clearly 
indicating the inhibitory effect of rotation on the onset 
of buoyancy-driven convection. We note that for V, = 

5 /i s , rotation at fi 0 = 500 rpm increases by slightly 
more than two orders of magnitude relative to the non- 
rotating case the critical Sn concentration above which 
the plane-front solution becomes unstable. We further 
note that the morphological branch is unaffected by ro- 
tation, whereas as ft 0 increases, the value of V, at which 
the onset of instability shifts from the convective branch 
to the morphological branch decreases from about 
40 p. s 1 in the nonrotating case to about 27 n s _1 for 
H 0 = 500 rpm. 

For each value of fl 0 , we see a local minimum near 
v i = 1ms with the minimum shifting to smaller growth 
velocities and becoming relatively more shallow as fl 0 
increases. We note that the maximum relative stabili- 
zation by rotation occurs near the local minimum, and 
that for Cl 0 = 500 rpm, the critical bulk concentration 
of Sn is increased more than a hundredfold. Although 
Coriell et al. m noted a local minimum in the stability 
boundary near V, = 1 ft s 1 in the nonrotating case for 
the largest gravitational acceleration considered, they of- 
fered no explanation for its existence. This minimum is 



Fig. 1 — V,-d stability boundaries for G, = 200 K cm"', with fi„ = 
0, 100, 200, 300, and 500 rpm. 
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a consequence of the fact that as V t — » 0, the concen- 
tration gradient G c = (k - ^PsC^Vj/ikp^D^ vanishes. 
As the temperature gradient is independent of c* and is 
stabilizing, the critical value of c* must ultimately in- 
crease as V; — ► 0. (Of course, the weight percent of the 
solute, Cx, cannot exceed 100.) As c* must initially de- 
crease with increasing V Jy there must be a local minimum 
on the convective branch before c * can increase to join 
up with the morphological branch of the stability bound- 
ary. This nonmonotonic dependence of c* on V f can be 
interpreted in terms of the existence of four critical val- 
ues of Vj for certain values of c*. (In addition to the three 
critical values shown in Figures 1 and 2, we note that 
for sufficiently large V h c * ultimately increases on the 
morphological branch/ 61 ) 

As discussed by Coriell et a// 21 for the nonrotating case, 
at each value of fl 0 , the onset of morphological insta- 
bility occurs via a short wavelength (large wavenumber) 
instability while convective instability sets in via distur- 
bances with relatively longer wavelengths. In the fol- 
lowing section, this point is discussed in the context of 
the mechanism by which rotation inhibits the onset of 
convection. 

For G l — 400 K cm -1 , Figure 2 shows stability bound- 
aries (c£ vs V/) analogous to those for G L = 200 K cm* 1 . 
Aside from a slight shift of the stability boundaries to 
higher values of c£, the results are qualitatively similar 
to those for G L = 200 K cm’ 1 . In particular, the onset 
of convection is suppressed but the morphological insta- 
bility is not influenced by rotation, and the critical value 
of c * passes through a minimum near V { = 1 /is’ 1 . 

VI. DISCUSSION 

The remarkable stabilization obtainable at low growth 
rates (more than a hundredfold increase in the critical 
value of c* at V) = 1 /is" 1 can be achieved by rotating 
the layer at 500 rpm for the two values of G L considered) 
is undoubtedly due to the well-known Taylor-Proudman 
mechanism described by Chandrasekhar/ 271 According 
to the Taylor-Proudman theorem, steady motion parallel 



Fig, 2 — Vf-ct stability boundaries for G L — 400 K cm \ with fi 0 = 
0, 100, 200, 300, and 500 rpm. 


to the axis of rotation in a uniformly rotating inviscid 
fluid is prohibited at any nonzero rotation rate. If this 
theorem were strictly applicable to a viscous fluid, the 
onset of steady convection would be prohibited, since 
the flow in convection cells must have a vertical com- 
ponent. Instead, in a viscous fluid, one sees an inhibition 
of the onset of steady convection, with the degree of 
inhibition (expressed here as an increase in c*) increas- 
ing with fl 0 . That the onset of oscillatory convection is 
hardly affected is due to the fact that the Taylor-Proudman 
theorem applies only to steady flows. It is also not sur- 
prising that the morphological instability is unaffected 
by rotation. The morphological instability occurs at very 
short wavelengths, so the motion is almost perpendicular 
to the solidification front ( i.e aligned with the axis of 
rotation). Hence, the Coriolis acceleration does not sen- 
sibly affect the morphological instability. 

For a horizontally unbounded binary fluid layer in which 
the density depends on temperature and one composition 
variable, Pearlstein’s linear stability analysis 1331 shows 
that Coriolis effects generally inhibit the onset of con- 
vection (by the Taylor-Proudman mechanism). Under 
some conditions, however, rotation can destabilize the 
layer, depending on the values of Pr, Sc, the dimension- 
less rotation rate (characterized by a Taylor number), and 
the dimensionless temperature or solute gradients (char- 
acterized by thermal and solutal Rayleigh numbers). For 
conditions under which destabilization (on a linear basis) 
occurs relative to the nonrotating case, instability sets in 
via an oscillatory mode, in which the natural frequency 
of oscillation of a buoyant fluid element is tuned (by 
rotation) in such a way that there is a local minimum in 
the critical value of Ra r as a function of the dimension- 
less rotation rate. This behavior was found r331 for Pr and 
Sc both less than unity, although there is no apparent 
reason why such destabilization cannot occur under other 
conditions when onset is via an oscillatory mode. 

Although Pearlstein [331 found in the rotating doubly 
diffusive case that for Pr < 1 < Sc (a condition satisfied 
in the present case, in which Pr = 0.023 and Sc = 81) 
there can exist as many as three critical values of the 
solute Rayleigh number for certain values of the Taylor, 
Prandtl, Schmidt, and thermal Rayleigh numbers, we have 
found no evidence of such multivalued stability bound- 
aries in the present calculations. As in the case investi- 
gated earlier/ 331 it is possible that such behavior occurs 
in relatively small regions of the parameter space (G L , 
V h etc.) and has gone undetected so far. As discussed 
in Section V, however, Figures 1 and 2 imply that for 
certain values of c^, there exist four critical values of V, 
(including the unshown portion of the morphological 
branch). 

Even though the foregoing analysis is restricted to a 
horizontally unbounded fluid layer, the work of Homsy 
and Hudson [291 and Biihler and Oertel 1471 suggests that its 
predictions will be qualitatively correct for finite aspect 
ratios (ratio of mold radius to height) if the parameter 
a 2 0 R Q /g (a Froude number, where R 0 is the mold radius) 
is sufficiently small. For the onset of thermal convection 
in rotating water or mercury layers heated from below, 
the excellent quantitative agreement between classical 
linear stability analysis for a horizontally unbounded 
layer [24,25 441 and experimental work for finite aspect 
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ratios 144 48 49] provides a clear demonstration of the 
potential of the Coriolis acceleration to suppress 
buoyancy-driven convection in a rotating fluid. 

Interpretation of the results of a stability analysis re- 
stricted to infinitesimal disturbances is obviously subject 
to the caveat that larger disturbances might grow, even 
though sufficiently small disturbances are predicted to 
decay. Indeed, it is known that for rotating fluid layers 
heated from below, the onset of thermal convection 
sometimes does occur f5 ° 51521 at lower Rayleigh numbers 
than predicted by linear theory. However, in that case, 
accounting for finite (i.e., noninfinitesimal) amplitude 
disturbances modifies the quantitative predictions of the 
theory; the basic qualitative prediction of stabilization by 
rotation remains unchanged. 

The relatively modest rotation rates required to sig- 
nificantly inhibit the onset of convection in the Pb-Sn 
system make the proposed method an interesting can- 
didate for a program of laboratory experiments. (We note 
here that the experimental work of Mtilleri 171 and Weber 
et al. m] focuses on the effect of rotation on the time- 
dependence of the supercritical flow, with no informa- 
tion given on the effect of rotation on the suppression of 
motion. Furthermore, although these authors have con- 
cluded that the Coriolis acceleration is key to the elim- 
ination of striations at high rotation rates, their 
experimental design complicates the separation of ef- 
fects of the Coriolis acceleration from the increased 
“pseudo-gravitational" effects associated with the cen- 
tripetal acceleration.) Experiments might be conducted 
using a completely filled cylindrical mold, thermally in- 
sulated on the vertical surface and mounted axi- 
symmetrically on a rotating horizontal turntable. If the 
liquid at the top of the rotating mold is in contact with 
a gas or vacuum, it will have a nearly paraboloidal free 
surface, on which the elevation above the point on the 
axis of rotation is A z = Vtlr 2 /(lg). Although this con- 
figuration is consistent with rigid-body rotation, it also 
leads to a nearly paraboloidal solid-melt interface. For 
high rotation rates or large mold radius, this will in turn 
lead to significant radial variations in the solidified alloy. 
Also, cooling at the radial boundary leads to a radial 
temperature gradient, which in turn leads to significant 
centrifugal effects. For these reasons, the results of the 
present analysis cannot be compared to the experimental 
work of Kou, f2n Kou et aL™ or Sample and Hellawell, r23] 
in which strong radial variations in macrosegregation are 
observed in the solid. Experiments of the type proposed 
above were initiated some time ago by Copley [531 for the 
crystallization of ammonium chloride from aqueous 
solution. 

Finally, we note that for many binary systems, plane- 
front solidification does not occur at practical growth ve- 
locities, and that the morphological instability results in 
dendritic solidification. In this case, rotation might also 
suppress buoyancy-driven convection in the melt and 
interdendritic liquid; this possibility is currently being 
investigated by us for the Pb-Sn system. 
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a magnitude of wavevector 

a x x-component of wavevector 
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*2 
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y 

e 

V 


©5 

K 
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y-component of wavevector 

disturbance concentration 

basic state concentration distribution 

amplitude of disturbance concentration 

bulk concentration 

critical bulk concentration 

diffusion coefficient 

magnitude of gravitational acceleration 

Gs/Gl 

(k ~ 1 )psC*V I /(kp L0 D L ) (nominal 
concentration gradient at interface) 
nominal liquid-side temperature gradient at 
interface 


(LVj + G L k L )/k s (nominal solid-side 
temperature gradient at interface) 
dimensionless height of computational 
domain in liquid and solid 
D L /V l 

unit vector in z-direction 
segregation coefficient 
thermal conductivity of liquid 
G L k L H/(LD L ) 

thermal conductivity of solid 
G L k s H/{LD L ) 

latent heat of fusion per unit volume 

D l /k l (Lewis number) 

m iGdG L 

slope of liquidus 

disturbance pressure 

v/k l (Prandtl number) 

position vector 

mold radius 

a TgG L H 4 /{k l v) (thermal Rayleigh number) 
<*cgc«H y /{D L v) (solutal Rayleigh number) 
v/D l (Schmidt number) 
time 

disturbance temperature in melt 
basic state temperature distribution in melt 
melting temperature of Pb 
f(G L H*) 

disturbance temperature in solid 

basic state temperature distribution in solid 

disturbance velocity 

basic state velocity 

nominal growth rate 

z-component of disturbance velocity 

amplitude of z-component of distui±>ance 

velocity 

cartesian coordinates 
dimensionless cartesian coordinates 
(2 Z + h)/h (vertical coordinate in solid) 

(2 Z - h)/h (vertical coordinate in liquid) 
solutal expansion coefficient 
thermal expansion coefficient 
amplitude of disturbed interface position 
(1 - k)/k 

Ps/Plo ~ 1 (fractional shrinkage) 
interface deflection relative to nominally 
planar condition 

amplitude of disturbance temperature in melt 
amplitude of disturbance temperature in solid 
*s/k l 

thermal diffusivity of liquid 
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k s thermal diffusivity of solid 

v kinematic viscosity 

P Ps/Plo 

p L local density of liquid 

p L 0 density of liquid at reference temperature and 
composition 

p s density of solid 

a temporal eigenvalue 

t dimensionless time 

'F capillary coefficient 

XI angular velocity 

fl 0 magnitude of (I 

ftj 2 Q, 0 H 2 /v (square root of the Taylor number 

Ta = 4SllH*/v 2 ) 

fi z amplitude of z-component of disturbance 

vorticity 

co z z-component of disturbance vorticity 
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Abstract 

A linear stability analysis is used to assess the potential of steady uniform rotation about the 
vertical axis to suppress the onset of convection during directional solidification of mercury 
cadmium telluride cooled from below. Since for sufficiently small x the density of Hg,_ x Cd x Te 
does not depend monotonically on temperature, there is a range of operating conditions (bulk 
CdTe mole fraction, nominal solidification rate, and liquid-side temperature gradient) for 
which, adjacent to the liquid-solid interface there is a sublayer in which the thermal 
stratification is destabilizing. This differs from the normal case (i.e., binary alloys 
characterized by a linear equation of state, such as Pb-Sn), in which the thermal stratification 
is stabilizing everywhere. Thus, in Hg 1 _ x Cd x Te melts there can exist a critical bulk CdTe mole 
fraction x x below which plane-front solidification is unstable at all solidification rates V p 
whereas in the normal case plane-front solidification at any x.. is linearly stable for 
sufficiently small V, Moreover, for Hg,_ x Cd x Te there can exist a critical solidification rate 
v , cr such that for Vj< V" plane-front solidification is unstable for all x B . Our results show 
that modest rates of uniform rotation can significantly suppress the onset of buoyancy-driven 
convection for the Hg 1 _ x Cd x Te system, and that rotation can substantially increase the critical 
solidification rate V^ r beyond which plane-front solidification is unstable for all bulk CdTe 
mole fractions. 
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1 . Introduction 

The crystal growth of mercury cadmium telluride (Hg,_ x Cd x Te) is of considerable interest 
due to its extensive use in the fabrication of infrared detectors and other electro-optical 
devices. This has led to a number of experimental [1,2] and theoretical [3-5] studies of the 
solidification of Hg,_ x Cd x Te from the melt. For a more extensive discussion, the reader is 
referred to the reviews by Micklethwaite [6] and Capper [7]. In those applications in which 
crystal size and defect density are critical, with the goal being to produce large, nearly perfect 
single crystals, it is frequently desired to achieve "plane-front” solidification in order to 
reduce segregation. If this could be accomplished, the process would be steady in a reference 
frame moving with the interface, and the only spatial variation would be in a direction normal 
to the interface. However, in real systems, several instabilities can cause departures from the 
nominally steady one-dimensional plane-front case. 

First, the liquid-solid interface may be subject to a so-called morphological instability, 
studied extensively since the early work of Mullins and Sekerka [8]. This instability deforms 
the nominally planar interface, ultimately leading to formation of a two-phase "mushy zone" of 
dendrites and interdendritic liquid. Departures from one-dimensionality and steadiness in the 
mushy zone can result in solute segregation in the solidified material. Second, depending on how 
the liquid density varies with temperature and composition, and according to whether the 
segregation coefficient is greater or less than unity (i.e., whether solute is rejected or 
preferentially incorporated at the interface), the melt may be subject to a buoyancy-driven 
convective instability. The resulting flow can contribute significantly to the redistribution of 
solute, and hence to segregation [9-1 2]. 

Convective and morphological instabilities in a binary liquid undergoing directional 
solidification were first studied by Coriell et al. [13] for Pb-Sn alloys. These authors showed 
that motion may occur due to either morphological or convective instabilities, and that the 
buoyancy force does not sensibly alter the criterion for onset of the former, which occurs at 
higher wavenumbers than the latter. Subsequent work has been reviewed by Glicksman et al. 
[14], and Sekerka and Coriell [15]. More general discussions of the effects of convection on 
plane-front and dendritic solidification have been given recently by Muller [12], Worster 
[16], Davis [17], Huppert [18], and Polezhaev [19]. 

One means by which the onset of buoyancy-driven convection in a density-stratified fluid 
layer can be inhibited is to subject the layer to a magnetic field aligned perpendicular or 
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parallel to the stratification [12]. For horizontal Pb-Sn layers solidified from below at 
several growth velocities, Coriell et al. [13] showed theoretically that the critical bulk Sn 
concentration above which instability occurs can be increased an order of magnitude by a 
vertical magnetic field of the order of one tesla. This technique requires the liquid to be an 
electrical conductor, and is applicable to metallic alloys, semiconductors, and aqueous solutions. 

Segregation might also be controlled by increasing or decreasing the magnitude of the 
gravitational acceleration or changing its direction. The theoretical studies of Coriell et al. 
[13], Alexander et al. [20], and Heinrich et al. [21,22] suggested that macrosegregation during 
directional solidification of binary alloys can be reduced by operating in a low-gravity 
environment. Muller [23] and Weber et al. [24] have recently discussed solidification under 
conditions where the magnitude of the body force is greater or less than that of normal gravity. 

Both gravity and an external magnetic field are body forces which act on the fluid. As 
opposed to contact forces such as pressure, viscous stress, and surface tension, which act on the 
surfaces of a fluid element, body forces act on the mass of a fluid element. Hence, their local 
strengths are proportional to the local fluid density. In addition to gravitational and magnetic 
body forces, there are other "pseudo-body forces" which manifest themselves as fictitious 
accelerations (centrifugal and Coriolis) when the reference frame to which the fluid motion is 
referred undergoes steady rotation relative to an inertial frame. (An additional fictitious 
acceleration associated with unsteady rotation relative to the inertial frame will not concern us 
here.) Although these accelerations do not correspond directly to forces (unlike the 
gravitational acceleration), they have the same mathematical form as accelerations associated 
with body forces, and can have dynamical consequences equally as profound as the gravitational 
and magnetic forces discussed above [25]. In light of this, and the fact that imposition of an 
external magnetic field or modification of the gravitational field is not always feasible, the 
possibility of using rotation to suppress the onset of convection in a liquid undergoing 
directional solidification is of interest. 

The idea of using rotation to affect fluid motion and solute distribution in the melt, and 
ultimately segregation in the solid, is not a new one. Beginning with the seminal work of 
Czochralski [26] on the use of differential rotation of a growing crystal relative to the crucible, 
there have been many attempts to use rotation to reduce compositional variations in crystal 
growth. Among these are the accelerated crucible rotation technique (ACRT) developed by 
Scheel and Schulz-DuBois [27] and Schulz-DuBois [28], in which axisymmetric azimuthal or 
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more complicated shear flows are generated by unsteady rotation of the crucible relative to the 
fluid. At about the same time, Copley [29] conducted preliminary experiments in which the 
solidification of aqueous ammonium chloride was conducted on a rotating turntable. Later, Kou 
[30], Kou et al. [31], and Poirier et al. [10] experimentally investigated the effects of rotation 
on macrosegregation in the solidification of Pb-Sn and Al-Cu alloys under conditions such that 
centrifugal effects were thought to be dominant. Sample and Hellawell [32,33] experimentally 
investigated the effects of mold precession on the solidification of Pb-Sn alloys and aqueous 
ammonium chloride solutions. More recently, Muller [23], Weber et al. [24], Rodot et al. 
[34], and Chen and Chen [35] have considered crystal growth in or on a centrifuge. 

Rotation can manifest its effects in a number of ways. The rotation of one or more solid 
boundaries can generate a shear flow (as in the A CRT technique) which directly affects the 
redistribution of solute through the advection term (u • Vx, where u is the velocity vector and 
x is the mole fraction) in the solute transport equation. Alternatively, if all solid boundaries 
rotate at the same steady angular velocity, then the fluid may be in solid-body rotation. In a 
rotating frame, the fluid is then motionless (u = 0) and there is no advective contribution to 
solute transport. However, the resulting fictitious accelerations can have important 
consequences for the stability of the solid-body rotation [25,36], either enhancing or 
suppressing the secondary flows responsible for solute redistribution. 

As discussed above, steady rotation at constant angular velocity Q about a fixed axis can 
influence the motion of a fluid by means of two "pseudo-body" forces. These correspond to the 
centripetal and Coriolis accelerations, for which the terms p^flxflxr and 2p L £Jxu, 
respectively, are added to the momentum equation. Here, r is the position vector measured 
from the axis of rotation, u is the local fluid velocity relative to the rotating frame, and p is 
the local density of the liquid. If the density variation is neglected, the centrifugal acceleration 
has no dynamical consequences. If the density variation cannot be neglected, then for a given 
geometry, the ratio of the magnitudes of these accelerations is 2u ± /(fl 0 r 0 ), where fi 0 , r , and 
u x are the rotation rate, characteristic length (e.g., mold radius), and component of fluid 
velocity relative to the moving frame and perpendicular to the axis of rotation, respectively. 
Thus, if u ± is independent of or increases sublinearly with fl 0 , the importance of the Coriolis 
acceleration relative to the centrifugal acceleration will increase with decreasing Cl Q . 

Oztekin and Pearlstein [37] have recently considered the effect of Coriolis acceleration on 
the stability of plane-front solidification of dilute Pb-Sn binary alloys, for which a linear 
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dependence of density on temperature and mass fraction is an excellent approximation to the 
equation of state. Their theoretical work shows that relatively modest rotation rates can 
significantly suppress the onset of buoyancy-driven convection in horizontally unbounded 
layers of that alloy undergoing directional solidification. 

The question naturally arises as to what effects rotation has on the onset of convection in the 
Hg,_ x Cd x Te system, in which the nonmonotonic dependence of liquid density on temperature 
[38] can give rise to, adjacent to the solid/liquid interface, a sublayer in which the thermal 
stratification is destabilizing. It is the purpose of the present work to investigate this question. 

We follow the approach of Oztekin and Pearlstein [37] and use a linear stability analysis to 
study the effects of the Coriolis acceleration on the onset of instability during solidification of 
H 9i-xCd x Te. The paper is organized as follows. In §2, we present the governing equations. The 
one-dimensional basic state and linear disturbance equations are given in §3. The latter are 
nondimensionalized and Fourier-decomposed in §4. The numerical solution technique is 
described in §5. Results for Hg, x Cd x Te are presented in §6, followed by a discussion in §7. 

2. Governing Equations and Formulation 

We adopt the solidification model used by Coriell et al. [13] and Oztekin and Pearlstein 
[37], in which the fluid properties, except density, are taken to be constant throughout the 
liquid. The variation of density with temperature and CdTe mole fraction will be accounted for 
in the buoyancy term and neglected everywhere else. The equations governing the fluid motion 
are the Oberbeck-Boussinesq equations in a reference frame translating with the nominally 
steady velocity ( i z V, ) of the moving interface and rotating at a constant angular velocity Q about 
the vertical axis 


V.u = 0, 


^ + u.Vu + 2Qxu + QxQxr = — Lyp-^k^lALj +V v 2 u, 


3t 


k Lo 

31 , 

-^ + u.VT l = k l v2 Tl , 

9x, , 

— - + u» Vx L =D L V 2 x L , 


Lo 


(la) 

(lb) 

(lc) 


- .~ L -- L ■ -L , (Id) 

where the dependent variables u, p, x L , and T L are respectively, velocity, pressure, CdTe mole 
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fraction, and temperature in the liquid, i z is the unit vector in the Z-direction, Vj is the 
nominal solidification rate (with dimensions of velocity), k l is the thermal diffusivity of the 
liquid, D l is the diffusivity of CdTe in the liquid, v is the kinematic viscosity, p u is the liquid 
density at a particular reference temperature and CdTe mole fraction, £2 = £2 0 i z is the constant 
angular velocity, and g is the magnitude of the earth's gravitational acceleration. We note that 
the fluid velocity u and interface velocity V s ( vide infra ) are referred to the moving frame, 
unlike the velocities employed in reference 13. Finally, we note that (Id) assumes that mass 
transport in the system can be accounted for using a binary model, in which CdTe diffuses 
through HgTe, and that the trace quantities of atomic Hg, Cd, and Te present in the liquid are 
negligible. 

We have used an equation of state based on the experimental data of Chandra and Holland 
[38] and Mokrovskii and Regel [39]. The variation of liquid density with temperature and CdTe 
mole fraction is represented by bivariate polynomials 


3 4-i 


P l (x,T) = II by x' (T - 1030 K) j 
i=o j=o 


for T<T max( X > 


(2a) 


and 


where 


3 1 

P l (x,T) = II dy x' (T - 1 030 K) j for T > T max (x) 

i=0 j-0 


(2b) 


3 

T max( x )= X f i Xi (3) 

i = 0 

is a polynomial fit to the largest temperature for which liquid density data are available in the 
literature. The linear dependence of density on temperature (2b) for T > T max (x) is chosen to 
match p L and BpJdJ at T max . Details are given in Oztekin and Pearlstein [40]. 

The energy equation in the solid is 


9T S 

at 


-V, 


2k 

az 


= k c V 2 T c 


( 4 ) 


where T $ is the temperature in the solid, and k s is the thermal diffusivity of the solid. 
The boundary conditions at the liquid-solid interface are 
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-(eV s +p* i z Vj) .n = u.n, 


(5a) 


(u+i z Vj).t = 0 (5b) 

where n and t are unit vectors normal and tangential to the interface, respectively, V $ is the 
local interface velocity, p* = P s /p Lo is the density ratio, p s is the (constant) solid density, and 
e = p*- 1 is the fractional shrinkage. The energy balance at the liquid-solid interface is 

l_( v s+ i z v i) • n = (- k L VT L+ k s VT s )» n , (5c) 

where the right-hand side represents the difference between the normal components of the heat 
flux vectors in the liquid and solid, and the left-hand side represents the product of the latent 
heat L(x oo ) and the normal component of the solidification velocity. Here k^x^) and k s (x oo ) are 
the thermal conductivities of the liquid and solid, respectively. The solute balance at the 
interface is 

Ps( x l" x s)( v s + «z v i) * n = -P L o D L Vx L* n * (5d) 


where x s is the CcTTe mole fraction in the solid, and solid-state diffusion has been neglected. We 
also require the temperature to be continuous 


W (Se) 

across the interface, the liquid and solid concentrations of CdTe at the interface to be related 
according to the binary phase diagram 

x s =kx L , (5f) 

where k is the segregation coefficient, and the temperature and CdTe mole fraction at the 


interface to be related by 


\“"^M + m L X L ^ 


f ± + y 
^ R 1 R 2 ) 


(5g) 


where T M is the melting temperature of HgTe, m L is the slope of the liquidus, R 1 and R 2 are the 
principal radii of curvature of the interface, and 'F is the capillary coefficient. The 
temperature- and composition-dependence of most of the thermophysical properties (other 
than p) is given in Appendix A. For the remaining properties, values measured at selected 
temperatures and compositions are taken as constant over the entire range. 


7 






v s =o, 


with 


u = -p*V, i z , 


3. Basic State and Linear Disturbance Equations 

In the reference frame described in §2, the steady one-dimensional basic state specified by 

(6a) 

(6b) 

K ' G ^[l - exp(-p* V, (6c) 

(6d) 

(6e) 


T, = T„ + 


m, x. 


k p*V, 


x L =x o 


l + i-^-exp(-p*V I Z/D L ) 

k 


in the liquid, and 


T s . T„+^=- + ^[1- <*P(-' V, Z/k s )], 


in the solid is, aside from the velocities, which we refer to a translating and rotating frame, the 
same as that used by Coriel, et a,. [13], Here, V„ 0. \, x L , and T s are, respectively, the 
basic state interface velocity, velocity, temperature and CdTe distributions in the liquid, and 
temperature distribution in the solid, and G L and G s are the liquid-s.de and solid-side 
temperature gradients, respectively. Our basic state differs from that of Cor, ell et a/. ,n t at 
v s = 0 and u is not proportional to the shrinkage e, because in our formulation the inter ace 

and fluid velocities are referred to the same (moving) inertial reference frame. 

As in stability analyses of other flows subject to uniform rotation [41-45], the I, near 
disturbance equations we use differ from those in the nonrotating case only by addition of the 
Coriolis acceleration to the disturbance momentum equation. (This approximation, discussed by 
Homsy and Hudson [48], neglects the consequences of density variation in the term 
corresponding to the centrifugal acceleration.) The equations governing small disturbances in 

the liquid are 

V.u=0, < 7a > 


W _ p* v + 2Q X u = - — Vp - f T r L i z + f x X' L i z + w 2 u' , 

dt p 1 az p L0 

_ p* v. +G L w'exp(-p* v, Z/k l ) =K L V 2 T L , 

3t 3Z 


(7b) 


(7c) 
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- P* V, + G X w' exp(- p* v, Z/D l ) =D l V 2 x' l 

at p 1 az x 


(7d) 


where u, V. T' L . *' L , and P are. respectively, the disturbances to the iiguid veiocity. ,ts 
Z-component, temperature, CdTe mole fraction, and pressure. Here, we e me 


T T" 3T 


f dp, (X,T) 


ax 


( 8 ) 


x l ,T l 


V T t 


and G x =(k-Dp*V,*_/(kD L ) is the CdTe gradient at the nominally planar interface. The 
disturbance energy equation in the solid is 

where T s is the temperature disturbance in the solid. 

Taking the curl of (7b) 


and the dot product of the result with i z , we obtain an equation 


_ P* v, ii - 2Q 0 ^ - vW 2 


( 10 ) 


at 


az 


az 


for the z-component of the vorticity, «' z . After twice 
component of the resulting equation is 

d(0 


taking the curt of (7b), the vertical 


. p * v, i-vV ♦ 2Q 0 %- =f T gv 2 n - ♦ W w 

at dz 6L 


(ID 


The disturbance boundary conditions 


linearized about the nominal interface position (Z 0) are 


u'(X,Y,0,t) = eV t g X Y — * 

v , (X,Y,0,t) = eV r &n( y ? -. 

_ % an(x,Y,t) 

w (X,Y,0,t) = _E— ’ 

dw'(X,Y,0,t) = _ eV V 2 tf(X,Y,t), 

az 


(12a) 
(12b) 
(12c) 
(1 2d) 


rwvvrf 8T,'(X,Y0,t) Jp^L.^ViKW). < 12e > 

L = _ kL _U_ + Ks-^ ^ *s J 
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( 12 f) 

(I2g) 

(12h) 


* \/2 X ;WX Y t) D, 8 xj (X, Y,0,t) 

i-^-T,(X,Y,t) + V I x' L (X,Y,0,t) "p* O-k) 3Z 

T'(X,Y.O,t) + G L n(X,Y,t>-T'(X,Y,0,t) *G s n(X.Y,t), 

T'(X,Y,0,t) + G L l\(X,Y,t)=m L G I n(X,Y,t)+m L x' L (X,Y,0,t) + T M 4'V^(X,Y,t), 

w here n is the interface deflection, v;(X,Y,t)= i, W* is the linearization of the disturbance 
interface velocity, and X and Y are horizontal coordinates. 

4 Nondimensionalization and Fourier Decomposition 

' we scale the velocity, time, length, temperature, CdTe mole fraction, and vortiaty with 
k /H, H 2 /v, H, G l H, x_. and * l /h 2 , respectively, where H = D L /v, is the ^aractenstic 
length [associated with the stabilizing gradient of preferentially incorporated CdTe, see < )], 

and write the horizontal and temporal dependence of the dimensionless disturbance quanbt.es as 
exp(ot ♦ ia x X* + ia v Y*), where a x and a, are the X- and Y-components of the wavevector X 
and Y* are the dimensionless horizontal coordinates, and o is the temporal eigenvalue. 
Substitution into (7) and (9-1 2) yields 

(13a) 


a(D 2 -a 2 )W-^-D(D 2 -a 2 )W = -F T a 2 0 L -F x a 2 x L +(D 2 -a z ) 2 W £1 0 DQ, 

Sc 

0 Q — DQ = Q*DW + (D 2 -a 2 )Q, 

Sc 

aPr0 L -p*LeD0 L = (D 2 -a 2 ) 0 L -Wexp(-p*LeZ ), 

■>* nv = fn 2 -a 2 W +^- 


(13b) 

(13c) 

(13d) 


oScx L -p* Dx l = (D 2 -a 2 )x L + 2 ^-Wexp(-p*Z < ) 
i the liquid, and 

a Pr0 s -LeD© s =K(D 2 -a 2 )0 s C ' 

, the solid. Here, Pr = v/* L is the Prandt, number, Scv/D is the Schmidt number, 

e = D /k is the Lewis number, F T = f T gG L H 4 /(x L v), F x -f x gx - H /(x L v), 0 0 

he square root of the Taylor number, k = * s /* u is the thermal diffusivity ratio, Z is the 

limensionless vertical coordinate, and T - <1-*)/k. The boundary conditions (12a-h) at the 
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interface become 


W(0) = -eapPr, 
DW(0) = ea 2 pi_e, 
Q(0)=0, 


aScp = -k L D9 L (0) + k s D0 s (O) + Le(k* p*-k* /k - Vk)P, 
(1-G)P=© s (O)-0 l (O), 


(ctSc + kp*)p = -kx L (0) Vdz l (0), 

YP 


(-a 2 T +m*-l)p = -2L_x (0) +0 L (O), 

Yp 


05a) 

05b) 

05c) 
0 5d) 

(1 5e) 
(1 5f) 

(15g) 


where W, X L , Q, and 0 L are, respectively, the amplitudes of the disturbances to the vertical 
velocity, concentration, vertical vorticity, and temperature in the liquid, 0 $ is the amplitude 
of the disturbance temperature in the solid, and P (a constant) is the amplitude of the disturbed 
interface position. Here we define dimensionless parameters k* = G L k L H/(LD L ), 
k s -G L k s H/(LD L ), m = m i_G x /G L , T M = T M 4y(G L H 2 ), and G = G s /G l . For the far-field 
boundary conditions, we follow Coriell et al. [1 3] and set all disturbances to zero 


W = DW = Q=x l = 0 l = 0 


©s = 0 


Z * 


Z " 


( 1 6a-e) 
0 6f) 


far from the interface. 


5. Numerical Solution 

Our objective is to find conditions under which infinitesimally small disturbances neither 
grow nor decay for a finite number of wavenumbers, and decay for all other wavenumbers. 
Disturbances which neither grow nor decay are said to be neutral. The neutral disturbances can 
be of two types, depending on the imaginary part of a. If the imaginary part of a is zero for a 
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neutral disturbance, the onset of instability will be via monotonicall y growing disturbances 

(steady onset). If the imaginary part of c is not zero, the neutral disturbance will oscillate in 
time (oscillatory onset). 

In the previous section, we formulated an eigenvalue problem for two systems of ordinary 
differential equations on two semi-infinite intervals (in the liquid and solid), coupled by 
boundary conditions at the deformable interface. For convenience, we follow Coriell et al. [13] 
and solve the problem on a finite interval [-h,h], where 2h is the dimensionless height of the 

computational domain. With regard to the far-field boundary conditions at Z*=±~, we set all 
disturbances to zero 


W = DW = Q=x L =e L =0 at Z*= h (17a-e) 

in the liquid and 

®s = 0 at Z*= -h (17f) 

in the solid. In this work, we have taken h to be at least 10, depending on the vertical structure 
of the basic state temperature and solute fields (6c, d). (For steady onset, we have checked a 
number of our results using more accurate asymptotic boundary conditions applied at Z*= ±h 
derived following Keller’s [46] procedure, and have found excellent agreement between the 
eigenvalues computed using the two sets of boundary conditions.) Since we use Chebyshev 
polynomials in our numerical solution, we scale the liquid and solid regions by z 2 = (2Z*-h) /h 
and z, = (2Z +h) / h, respectively, so that each region lies between -1 and +1 . The resulting 
system is then solved using a spectral Galerkin technique developed by Zebib [47], The problem 
is thus reduced to a matrix eigenvalue problem 

As + oBs = 0, (18) 

where o is the temporal eigenvalue, and the elements of the square matrices A and B depend on 
the square of the horizontal wavenumber a 2 - a 2 + a* , the bulk CdTe mole fraction x_ , and the 
other dimensionless parameters. 

In §6, we characterize the stability of the nominally plane-front solution in terms of the 
bulk CdTe mole fraction x m and solidification rate V p with all other parameters taken as fixed. 
As shown schematically in figure 1 , there is a critical value of V, (denoted by V cr ) such that 
for Vj > Vj , some disturbances (i.e., for some wavenumbers) grow for all values of x.,, while 
for Vj< Vj cr , disturbances of every wavenumber decay for some range of bulk CdTe mole 
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fractions. The upper and lower limits of that range are denoted by x‘ and x", respectively. To 
determine the critical values of x_ ( «; and x”) we first find the morphological and convective 
neutral curves ( x. versus a) separating those combinations of «_ and a for which all temporal 
eigenvalues a lie in the left half-plane (LHP) from those for which at least one eigenvalue lies 
In the right half-plane (RHP), corresponding to stable and unstable basic states, respectively. 

in order to determine the bulk CdTe mole fraction on the morphological and convective 
neutral curves for an arbitrarily chosen wavenumber a, we first compute all eigenvalues a of 

< 1 81 at each ° f N values of x - ( x - ■ 1 5 n < N) in the range b£*\0.2], where the lower bound 
is typically 1 0- 4 and the upper bound Is determined by the largest CdTe mole fraction for which 
we have data for the equation of state. We then attempt to determine a range of x. for which all 
temporal eigenvalues are in the LHP. If one of the original N values of x. selected is stable 
(i.e„ all temporal eigenvalues lie in the LHP), we then determine two intervals such that as x. 
increases, within one a transition from instability to stability occurs, and within the other a 
transition from stability to instability occurs. Through these two intervals pass the convective 
and morphological neutral curves, respectively. If none of the original N values of x_ is stable, 
we choose the value (say, x<?) for which the most unstable temporal eigenvalue has the 
smallest real part, and subdivide the interval [x^x^*')] unti | we either find a stab|e va|ue of 

x. (at which juncture we proceed to isolate the two intervals described above), or abandon the 

search when the real part of the least stable temporal eigenvalue and the difference between 

consecutive values of x. supports the expectation that all intermediate values of x. are 

unstable. Having found the intervals of x_ in which the transitions occur, we then compute the 

bulk CdTe mole fractions on the convective and morphological neutral curves using a bisection 
method. 

To compute the critical bulk CdTe mole fraction *1 (i.e., the minimum on the 
morphological neutral curve), we arbitrarily choose a wavenumber and compute the 
corresponding x„ on the morphological neutral curve using the procedure described above. We 
then fix x„ at the value computed in the previous step, and compute the eigenvalues a for a 
discrete set of wavenumbers in a chosen range. We next select the wavenumber from among this 
set at which the a with largest real part was obtained. If this wavenumber is at an endpoint of 
the chosen range, we extend the range until the wavenumber corresponding to the largest Re(a) 
is inside the range. We then select this wavenumber and determine the corresponding x„ on the 
neutral curve. We continue this process until the relative change in x m is less than HT 6 . We 
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compute x** following the same procedure. 


6. Results 

The solutions of (13-16) depend on a number of parameters, as defined in §2-4 
(Solutions on a finite interval depend also on h.) Restriction of the analysis to the pseudobinary 
system Hg,_ x Cd x Te reduces the parameters to the solidification rate V,, liquid-side 
temperature gradient G L , bulk CdTe mole fraction x_, and rotation rate 0 The 
thermophysical properties of the liquid and solid (other than liquid density) are taken as 
constants and evaluated using the functional forms given in Appendix A at the bulk CdTe mole 
fraction and corresponding liquidus or solidus temperatures. 

We present our results in figures 2 and 3 in terms of stability boundaries in the V -x 
plane for G L = 25 and 50 K cm -1 and several rotation rates Q 0 . The stability boundaries 
consist of some portion of the morphological branch joined to an oscillatory convective branch. 
The solid curve in each figure represents the minimum on the neutral curve corresponding to 

morphological instability (x^). The dashed curve denotes the maximum on the neutral curve 

associated with the onset of buoyancy-driven convection $<**). 

For G L = 50 K cm \ figure 2 shows the stability boundaries for £2 C = 0, 25, 50, and 100 
rpm. Above the morphological branch, disturbances grow for some range of wavenumber a and 
plane-front solidification is unstable. Similarly, below the convective branch appropriate to 
each rotation rate shown, disturbances grow for some range of a. However, between the 
morphological branch and the convective branch associated with each rotation rate (i.e„ for 

*" ' **' X -*' dlsturban ces decay for all wavenumbers and plane-front solidification is stable 
Note that the morphological branch is unaffected by rotation, whereas as n„ increases, the 
convective branch is shifted downward. The critical bulk CdTe mole fraction x" on the 
convective branch is a decreasing function of at any solidification rate, clearly Indicating 
the inhibitory effect of rotation on the onset of buoyancy-driven convection. Beyond the critical 
solidification rate (denoted by V") at which the convective and morphological branches 
intersect, there is no stable range of bulk CdTe mole fraction x_. Thus, plane-front 
solidification is unstable at all bulk CdTe mole fractions for sufficiently high solidification 
rates. The critical solidification rate at which the morphological and convective branches 
intersect occurs at higher x. with increasing rotation rate Q„. The critical value V" is 
increased by more than a factor of ten at 0 o = 100 rpm relative to the nonrotating case. We' also 
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note that the convective branch has a relatively shallow local minimum near V, = 1 nsec’ 1 , the 
location of which depends only weakly on Q 0 . Thus, for a given rotation rate, operation at the 
solidification rate corresponding to this local minimum allows plane-front solidification to be 
conducted stably at the lowest value ofx^. 

For G l = 25 K cm -1 and fl 0 = 0 and 100 rpm, figure 3 shows stability boundaries in the 

V I ~ X "° P ' ane analogous t0 those for G l = 50 K cm' 1 . As expected on the basis of the results of 
Conell et al. [13] for Pb-Sn without rotation, the morphological branch is shifted to the left, 
reducing the range of solidification rates and bulk mole fractions for which plane-front 
solidification is stable On the other hand, the convective branch is shifted downwards relative 
to the G l = 50 K cm' 1 case, corresponding to an increase in the range of stable operating 
conditions. We note that decreasing the temperature gradient has very little effect on the onset 
of convection in a nonrotating layer, but reduces the range of stable bulk CdTe mole fractions by 
more than twofold for a layer rotating at 100 rpm. The critical value of V t is higher by more 
than a factor of five at fi 0 = 100 rpm relative to the nonrotating case, although the degree of 
stabilization is less than the factor of ten predicted for G L = 25 K cm" 1 . Although the reduction 
in x„ has very little effect on the critical value of V, at which the morphological and convective 
stability boundaries intersect in the nonrotating case, for Q q = 100 rpm the critical value of V, 
decreases by a factor of three when the temperature gradient is reduced. 

As discussed for the Pb-Sn case [37], at each value of the onset of morphological 
instability occurs via a short wavelength (large wavenumber) instability, while convection 
sets in via disturbances with relatively longer wavelengths. In the following section this point 
is discussed in the context of the mechanism by which rotation inhibits the onset of convection. 

7. Discussion 

During directional solidification of binary alloys cooled from below and characterized by a 
linear equation of state, with rejection of a light solute at the liquid-solid interface (e.g., Pb- 
Sn), the vertical temperature and solute gradients are stabilizing and destabilizing, 
respectively, throughout the liquid layer. In the Pb-Sn system, the onset of buoyancy-driven 
convection occurs via monotonically growing disturbances, as shown by Coriell et al. [13], 
However, for the pseudobinary Hg^CdJe system, the CdTe gradient is everywhere 
stabilizing, and adjacent to the liquid-solid interface there can exist a sublayer in which the 
thermal stratification is destabilizing. In this system, the onset of convection occurs via 
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oscillator disturbances. We note that, as shown in Figures 2 and 3, the system is unstable at 
the indicated values of G L for al, V, when x. is sufficiently small, corresponding to the absence 
of a stably stratified thin liquid layer adjacent to the interface for small x . 

For Pb-Sn alloys we have shown [37] that significant suppress^ of the onset of 

buoyancy-driven convection occurs at modes, rotation rates, and is due to the well-known 

ayor-Proudman mechanism. According to the Taylor-Proudman theorem [36], asady motion 

para e to t e axis of rotation in a uniformly rotating inviscid fluid is prohibited If this 

theorem applied to viscous fluids, steady convection would be prohibited, since cellular 

convective flow must have a vertical component. Instead, in viscous fluids the onset of steady 

convection is inhibited, with the degree of inhibition (expressed as an increase x , 

increasing with rotation rate. Although the Taylor-Proudman theorem strictly applies only'to 

Steady flow, the Coriolis acceleration also suppresses the oscillatory onset of buoyancy-driven 

convection during directional solidification of Hg^Te as shown in §6. It is no, surprising 

that the morphological instability is unaffected by rotation since i, occurs at very short 

wavelengths, w„h the motion nearly perpendicular to the solidification front (i.e„ aligned with 
the rotation axis). 

For a horizontally unbounded fluid layer in which density depends linearly on temperature 
and composition and no phase change occurs, a linear analysis [45] shows that Coriolis effects 
generally inhibit the onset of convection (by the Taylor-Proudman mechanism). However for 
some combinations of Pr, Sc, and the dimensionless rotation rate (characterized by a Taylor 
number Ta) and solute and temperature gradients (characterized by solutal and thermal 
Rayleigh numbers Ra s and Ra T ), rotation destabilizes the layer (on a linear basis, relative to 
the nonrotating case. This occurs only when convection sets in via oscillations and the natural 
requency of a buoyant fluid element is tuned (by rotation) so that the critical Ra T has a local 
minimum as a function of Ta. With linear gradients and no phase change, this destabilization is 
found for Pr and Sc both less than unity [45], although there is no apparent reason why it 
cannot occur under other conditions when onset is via an oscillatory mode. We have found no 
such destabilization by rotation for the solidification of Hg,_ x Cd x Te. 

Even though our analysis is limited to horizontally unbounded layers, work for single- 
component fluids not undergoing phase change [48,49] suggests that for fixed finite aspect 
ratios (ratio of mold radius to heigh,) our predictions will be qualitatively correct if R L 
(a Froude number, where R„ is the mold radius) is sufficiently small. For the onset of , “hemal 
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convection in rotating water or mercury layers heated from below, the excellent quantitative 
agreement between linear stability theory for a horizontally unbounded layer [41-43] and 
experimental work for finite aspect ratios [43,50,51] clearly demonstrates the potential of the 
Coriolis acceleration to suppress buoyancy-driven convection in a rotating fluid. 

The results of a linear stability analysis are subject to the caveat that even when 
infinitesimal disturbances are predicted to decay, larger disturbances might grow. Indeed, for 
rotating fluid layers heated from below, the onset of thermal convection can occur at lower 
Rayleigh numbers than predicted by linear theory [52-54]. However, in that case accounting 

for disturbances of noninfinitesimal amplitude modifies only the degree of stabilization 
predicted. 

The modest rotation rates that can significantly inhibit the onset of convection in the 
H 9|-x Cd x Te and pb ‘Sn systems make uniform rotation about the vertical a good candidate for 
experimental verification. Experiments should be conducted with completely filled cylindrical 
molds, insulated on the vertical surface, and mounted on the axis of a rotating horizontal 
turntable. If these precautions are not taken, the Coriolis-related stabilization may be 
overwhelmed by other effects. For example, if the liquid's top surface is in contact with a gas or 
vacuum, it and the liquid-solid interface will be nearly paraboloidal, which may lead to radial 
segregation in the solid for high rotation rates or large mold radii. Cooling the vertical surface 
leads to a radial temperature gradient and centrifugal effects. For these reasons, our results 
cannot be compared to earlier experimental work [30-33]. 

We note that the experimental work of Muller [23] and Weber et al. [24] focuses on the 
effect of rotation on the time-dependence of the supercritical flow, with no information given on 
suppression of the onset of motion. Furthermore, although these authors concluded that the 
Coriolis acceleration is key to eliminating striations at high rotation rates, their experimental 
design complicates the separation of Coriolis effects from the "pseudo-gravitational" effects 
associated with high centrifugal accelerations. 

Finally, Antar [55] has recently presented a linear analysis of the onset of convection in a 
fluid layer cooled from below, using an approximate equation of state for Hg,_ x Cd x Te. Although 
this work purports to consider "convective instabilities in the melt for solidifying mercury 
cadmium telluride", it differs from ours in that it includes no phenomena associated with 
solidification (existence of a moving or deformable interface where phase change occurs, 
nonlinear basic state solute stratification due to rejection or preferential incorporation at the 
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interface, latent heet effects, etc, Comparison to its resuits is further precluded since the 
..near basic state therma, stratihcation and the definition of the temperature difference AT 
across the fluid layer are inconsistent. Specifically, the bottom temperature (see Antar's eqn 
( » is T„-aT, which is inconsistent with the definitions of T, (temperature where a p / 3 T is 
zero, determined by the equation of state) and at, Pu , differently, Antar defines AT as the 
temperature difference across the layer, (e.g., aT = |t -T I . | T ( d ) - TfOih „■ 

temperature profile, however, yields T(d, -T 0 -(d 0 -d)AT/d o and T^It -aT, from which it 
follows that |T(d) - T(0)| = d AT /d 0 , which is inconsistent wth the dehnition of AT This error 
vitiates the analysis, which uses a therma, Rayleigh number based on a temperature difference 
inconsistent with the basic state. That this inconsistency is not simply an apparent one due to a 

blT/sTt riT ' S C ° nf ' rmed by re,erence to an ear lier paper [56] which uses the same 
asic state and disturbance equations. A different (and very unusual) definition of at (related 

to the overall temperature difference by a consent dependent on the top or bottom temperature; 

see igure of [56]) was used to maintain consistency. That definition of at was the only one 

compatible w„h the basic state, and hence with the anaiysis and results. Unfortunately, in 

re erence 55, AT ,s defined as the overall temperature difference, which is wrong. 
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Appendix A 

The thermophysical properties of the liquid and solid (other than the density of the liquid) 
are taken as constants evaluated using the functional forms shown below at the bulk mole 
fraction and corresponding liquidus or solidus temperatures T |jq and T ). 

We have taken the density of solid Hg 1 _ x Cd x Te as 

_L = _L_ 1-c 


_ an d P 
CcTTe *HgTe 


r- k 1/ (A1) 

S CdTe K HgTe 

where C - xM cdTe /M is the CdTe mass fraction, the molecular weights of CdTe and HgTe are 
W 240 9/m0 ' and 328.2 g/mol, M=xM OTe+0 - x) u and p 
are the solid densities of pure CdTe and HgTe, respectively, represented by 

p CdTe = W (A2a) 

and 

P HgTe~ **0 + V’ (A2b) 

Where the coefficients X, and ♦, are determined by least-squares fits to the experimental data of 
Glazov et al. [57] and Mokrovskii and Regel [39], respectively. The values of k, and ♦, are 


= -8. 095 x 1 0 5 g cm -3 K -1 
V -1-230 x 10 -4 gcm -3 K -1 


(A3) 


= 5.820 g cm -3 
< f> 0 = 8.201 g cm -3 

The viscosity is represented over the entire range by 

v = 1 x 1 0“ 2 exp^o-KtyT) cm 2 sec-\ 
where * 0 = -3.401 and 3445 K are computed by least-squares fit to the experimental datl 
of Glazov et al. [58] for pure CdTe. Due to the lack of experimental viscosity data for HgTe we 
use (A3) for the viscosity of Hg^CdJe, independent of composition. 

The thermal diffusivity of liquid and solid Hg^CdJe are approximated by 

k l = (^o + EiX 1 / 2 +E 2 x)ln(T/K) - (U 0 + U 1 x 1 / 2 +(J 2 x), 


where the coefficients 

E o 


j=0 


(A4a) 

(A4b) 


E n = 0.104 cm 2 sec -1 


Ej= -0. 1 46 cm 2 sec 
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9=^ 


E 2 = 0.118 cm 2 sec _1 

U,= -0.941 cm 2 sec -1 

^oo = “5.945 x10 -3 cm 2 sec -1 

S 2 o = “ 0114 cm 2 sec -1 

S 1 1 = 1. 646 x 1 0" 4 cm 2 sec -1 

5q2 = “7.829x10 -8 cm 2 sec -1 K -2 


U 0 = 0.668 cm 2 sec _1 
U 2 = 0.783 cm 2 sec- 1 
S 10 = 1-698 xl0~ 2 cm 2 sec _1 
5 01 = -7.148xl0~ 3 cm 2 sec~ 1 K~ 1 
S 21 = 4.475x1 0 -4 cm 2 sec _1 K _1 
5 12 = 1.969x1 0 -7 cm 2 sec -1 K~ 2 


S 22 - “4. 120x10 7 cm 2 sec ’K -2 

are determined by least-squares fit to the experimental data of Holland and Taylor [57], 
The specific heat of liquid Hg lx Cd x Te is represented by 


2 2-i 


c pj_ (T, x) = £ £ q jj x j (T - 943 K) j . 

i=0 j=0 


(A5 


The coefficients q/ound by least-squares fit to the calculated specific heat of Su [60] are 

q 00 = 65.08 J K -1 mol -1 q 01 = -0.101 J K" 2 mol' 1 

q 02 = 3.961 x 10- 4 JK- 3 mol" 1 q 1Q= ,.648 JK" 1 mol"’ 

q n = 7.004 x 10- 3 J K^mor 1 q 20 = -5.044 x 1 0" 3 J K" 1 mol" 1 , 

where we have used the relation 1 mol - 2 g-atom for the pseudobinary system Ha Cd Te t 
convert Su's data (in cal K" 1 g-atom" 1 ) to J K" 1 mol" 1 . * 

The specific heat of solid Hg^CdJe is taken as 

C P,S = X ) C p,S,HgTe + XC p,S t C<fTe> (A6j 

where the specific heats of pure HgTe and CdTe vary with temperature [61 ] according to 


and 


p,S,HgTe 


“p.S.CcJTe 


52.09 + 9.08 x 1 0" 3 (T / K)J J K^mol" 1 
40.0 + 3.3 x 10~ 2 (T/KjJjK-’mor 1 . 


(A7a 


(A7b; 


The latent heat of fusion is approximated by 


L(x) - ( 1 x)L HgTe + x Lccjje , 


where L HgTe = 3.6 x 1 0 4 J / mol and L HgTe = 4 x 1 0 4 J / mol 


(A8) 

are used for the pure components 
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[61] The Werfacia' energy is estimated by comparing Hg,. x Cd x Te and InSb (in reference 62 
Ge and InSb are compared to estimate the interfacial energy of InSb) 





f 7 \ 


L 


a 2 

a H gi . x Cd x Te- 

a 2 


0 

L 


°'Hg^ x Cd x Te 

^ A 


J lnSb- 


(A9) 


where the latent heat of fusion per uni, volume for Hg,_ x Cd x Te Is taken to va,y linearly with 
CdTe mole fraction (A8), a 0 . 6.4797 A for InSb. and a, . 6.465 A is used for Hg,_ Cd Te 
in ependent of x [58], Here the surface tension is taken as ff ha> . 8. 5 x 1 O' 6 J cm' 2 , and the 
latent heat of fusion per unit volume of InSb is L = 1 20 x 10 3 J cm' 3 Th» „ 

UUXIUJcr "- ' he capillary coefficient 

is then computed from 


^ CT Hg Vx Cdj e -j^ . 

where T, jq is the liquidus temperature for a given bulk mole fraction of CdTe in the liquid. 
We determine the liquidus and solidus temperatures 


(A 10) 



T| jq = (943 + 681 x - 372 x 2 ) K , 

(Alla) 

the liquidus slope 

T S oi = (943 + 202 x - 1 50 x 2 + 324 x 3 ) K , 

(A 1 1 b) 


m L = (681 - 745 x)K, 

(A 12) 

and the segregation coefficient 


k = x(0.30 - 2.24 x lO^K-’T +2.67 x 10- 6 K- 2 T 2 ,) 

SOI ' 

k = 3.74 


for x> 0.1, (A13a) 

for x < 0. 1 , (A 1 3b) 


as functions of CdTe mole fraction by least-squares fits to data of Szofran and Lehoczky (63). 
Finally, the solute diffusion coefficient is taken as D L . 5.5 x HT 5 cm 2 /sec, the value most 


commonly used in the literature [3], 
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Figure Captions 

Schematic depiction of the stability boundary. 

V,-x. stability boundaries for G L = 50 K cnr’ with £2,, = 0, 25, 50, and 100 rpm. 
v,-x. stability boundaries for G L - 25 K cm-' with Q„ - 0 and 100 rpm. 
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Abstract 

The convective and morphological instabilities In a horizontally unbounded binary liquid 
undergoing directional solidification by cooling from below are studied by means of linear 
stability analysis. The possibility of using steady uniform rotation about a vertical axis to 

suppress the onset of buoyancy-driven convection is considered for both plane-front and 
dendritic solidification. 

For dilute Pb-Sn alloys, our results clearly show that the onset of convection in a 
horizontally unbounded layer undergoing plane-front solidification can be suppressed 
significantly at modest rotation rates. Specifically, plane-front solidification is linearly 
stable at higher Sn concentrations In a rotating configuration than in a nonrotating one. The 
predicted inhibitory effects of rotation on convection are discussed in terms of previous 
experimental and theoretical studies of the effect of rotation on the onset of buoyancy-driven 

convection in single-component fluids heated from below and in binary fluids subject to 
thermal and solutal stratification. 

We also consider the stability of one-dimensional plane-front solidification of the 
pseudobinary Hg^Cd/Te system (of considerable importance due to its wide application in the 
fabrication of electro-optic detectors) in which the liquid density does not depend 
monotonically on temperature (e.g., has a local maximum) for some range of the bulk solute 
composition. In contrast to the normal case where the density depends monotonically on 
temperature and composition (e.g., the Pb-Sn alloys considered by Coriell et al. (1980)), 
for certain combinations of the operating parameters (solidification rate, nominal liquid- 
side vertical temperature gradient, and bulk solute concentration) there can exist a critical 
value of the bulk mole fraction (CJ below which plane-front solidification is unstable at all 
dimensionless solidification rates y, whereas in the normal case plane-front solidification at 
any C„ is linearly stable for all sufficiently small dimensionless solidification rates. 
Moreover, when the density varies nonmonotonically with temperature, there can exist a 
critical value of the dimensionless solidification rate y c such that for y> y c plane-front 
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solidification is unstable for all C m . In this case, for y <y c there is a finite range of C„ for 
which plane-front solidification Is stable. This latter result differs from the normal case, 
for which at all dimensionless solidification rates, plane-front solidification is stable for all 
values of C„ lying below some critical value. The stability boundaries and neutral curves, 
differing qualitatively from those for the normal case, are discussed in terms of the existence 
adjacent to the liquid-solid interface of a sublayer in which the thermal stratification is 
destabilizing. For Hg^CdJe. our results show that uniform rotation at modest rotation rates 
can significantly suppress the onset of buoyancy-driven convection. 

For a binary liquid undergoing solidification by cooling from below, we assess the 
stability of the one-dimensional dendritic solution. The mushy zone, consisting of liquid and 
solid phases, is modeled as a porous medium with anisotropic permeability. The local 
porosity, as well as the location of the boundaries separating the solid from the mushy zone 
and the mushy zone from the liquid, are taken to be dynamical variables. The basic state, 
computed using a thermodynamically self-consistent nonlinear model of solidification, exists 
for only some combinations of the operating parameters (solidification rate, nominal vertical 
temperature gradient, and bulk solute concentration). The dendritic solution also exists 
under conditions for which plane-front solidification is linearly stable with respect to 
morphological disturbances. This point is discussed in the light of previous work on 
nonlinear morphological instability. A division of the parameter space according to the 
existence and stability of solutions corresponding to plane-front and dendritic solidification 
is presented and discussed for the Pb-Sn system. Uniform rotation is shown to be less 
stabilizing than in the plane-front case. 
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CHAPTER 1 

r 

Introduction 

1.1. Motivation and Previous Work 

Many solid materials of practical interest are prepared by directional solidification of 
binary and multicomponent liquids. Examples include high-strength, high-temperature 
alloys for use in turbine blades, electronic and photonic materials, and protein crystals for 
use in x-ray crystallography. 

In those applications in which the goal is to produce large, nearly perfect single crystals, 
it is frequently desired to achieve "plane-front" solidification. In such a situation, the 
liquid-solid interface remains planar, the temperature, concentration, and pressure fields 
are one-dimensional with their gradients normal to the interface, and the fluid motion (due 
solely to shrinkage) is in the direction of the interface motion and independent of position. 
The process is steady in a reference frame moving with the interface. It is easily shown that a 
necessary condition for the existence of this steady one-dimensional state is that the density 
gradient be parallel or anti-parallel to the gravitational acceleration. Even when this 
necessary condition is satisfied (e.g., by cooling the liquid from above or below, rather than 
from the side), however, several instabilities can cause departures from the nominally 
steady one-dimensional plane-front case. 

First, the liquid-solid interface is subject to a morphological instability (Mullins & 
Sekerka 1964) deforming the planar interface and ultimately leading to formation of a two- 
phase “mushy zone" of dendrites and interdendritic liquid. In turn, departures from one- 
dimensionality and steadiness in the mushy zone result in nonuniform solute distribution in 
the solidified material. Second, the density of a binary liquid depends on both temperature and 
composition. When a liquid is solidified by cooling from below, preferential rejection or 
incorporation of solute at the interface is potentially destabilizing if the solute-enriched (or 
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depleted) liquid at the interface is less dense than the warmer overlying bulk liquid. This 
solute gradient can overcome the stabilizing temperature gradient, leading to buoyancy- 
driven convection in the liquid, thus providing another mechanism, besides molecular 
diffusion, for redistributing solute from the interface into the bulk liquid. 

Convective and morphological instabilities in a binary liquid undergoing directional 
solidification were first studied by Coriell et ah (1980). They showed that the buoyancy 
force does not sensibly alter the criterion for onset of the morphological instability, which 
occurs at shorter wavelengths than does the convective mode. Subsequent work has been 
reviewed by Glicksman et ah (1986) and Sekerka & Coriell (1987); recent work includes 
papers by Polezhaev (1988), Muller (1988), Young & Davis (1989), Huppert (1990), 
Davis (1990), and Worster (1991). 

Buoyancy-driven convection during directional solidification of binary alloys has been 
shown to be the dominant factor in the formation of "freckles", “channel segregates , and 
other macrosegregation defects deleterious to the mechanical properties of directionally 
solidified alloys. The formation and characterization of freckles in nickel-based superalloys 
were first studied experimentally by Giamei & Kear (1970). Following their work, Poirier 
et ah (1981) investigated macrosegregation in electroslag ingots, showed that convection in 
the melt results in freckling in the solidified material, and suggested that rotation might 
reduce freckling. Sarazin & Hellawell (1988) have experimentally demonstrated the 
formation of freckles in Pb-Sn, Pb-Sb, and Pb-Sn-Sb ingots. Specifically, compositional 
convection, in the form of discrete plumes emanating from chimneys in a mushy zone of 
dendritic crystals, has been observed by Copley et ah (1970) and Chen & Chen (1988) in 
laboratory experiments using transparent aqueous ammonium chloride solutions. The 
resulting convection leads to freckles strikingly similar to those found in ingots cast of 
metallic alloys. Ridder et ah (1981) studied the effects of fluid flow on macrosegregation in 
nominally axisymmetric ingots and showed that melt convection results in macrosegregation 
in the solidified materials. In a theoretical study of a binary alloy solidifying radially 
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inward, Maples & Poirier (1984) concluded that macrosegregation results from mushy-zone 
convection driven by nonuniform temperature and solute distributions. A recent review of 
the effects of buoyancy-driven convection on macrosegregation in binary and pseudobinary 
systems has been given by Muller (1988). In particular, much attention is currently 
focused on the questions of when and how the freckles and channel segregates and other 
macrosegregation defects are formed, with the main goal being to grow materials that are free 
of such macrosegregation defects resulting from convection during solidification. 

One means by which the onset of convection can be inhibited in a density-stratified fluid 
layer is to subject the layer to a magnetic field aligned parallel or perpendicular to the 
stratification (MUIIer 1988). For horizontal Pb-Sn layers solidified from below at several 
growth velocities, Coriell et al. (1980) showed theoretically that the critical Sn 
concentration above which instability occurs can be increased an order of magnitude by 
applying a vertical magnetic field of the order of one tesla. This technique requires that the 
liquid be an electrical conductor, and so is applicable to metallic alloys, semiconductors, and 
aqueous solutions. 

Macrosegregation might also be controlled by increasing or decreasing the magnitude of 
the gravitational acceleration or changing its direction. Alexander et al. (1989) and Heinrich 
et al. (I989a-b) theoretically studied the effect of reduced gravity on macrosegregation in 
directionally solidified alloys. These authors suggested that macrosegregation in alloys can be 
reduced by solidification in a low-gravity environment. MOIIer (1990) and Weber et al. 
(1990) have recently discussed solidification under conditions where the magnitude of the 
body force is greater or less than that of normal gravity. 

Both gravity and an external magnetic field are body forces which act on the liquid. As 
opposed to contact forces such as pressure, viscous stress, and surface tension, which act on 
the surfaces of a fluid element, body forces act on the mass of a fluid element. Hence, their 
local strengths are proportional to the local fluid density. In addition to gravitational and 
magnetic body forces, there are other "pseudo-body forces" which manifest themselves as 
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fictitious accelerations (centripetal and Coriolis) when the reference frame to which the 
fluid motion is referred undergoes steady rotation relative to an inertial frame. (An 
additional fictitious acceleration, with which we will not be concerned, manifests itself if the 
rotation of the noninertial frame relative to the inertial frame Is unsteady.) Although these 
accelerations do not correspond directly to forces (as in the case of the gravitational 
acceleration), they have the same mathematical form as accelerations associated with body 
forces, and can have dynamical consequences equally as profound as the gravitational and 
magnetic forces discussed above (Greenspan 1968). In light of this, and the fact that 
modification of the gravitational field or imposition of an external magnetic field may not 
always be feasible, the possibility of using rotation to suppress the onset of convection in a 
liquid undergoing directional solidification is of interest. 

Several effects of rotation on solidification have been discussed by Schulz-DuBois 
(1972). The experiments of Kou (1978), Kou et al. (1978), Sample & Hellawell (1984), 
Muller (1990), and Weber et al. (1990) have shown that rotation can significantly reduce 
the degree of macrosegregation in binary alloys directionally solidified under plane-front or 
dendritic conditions. Kou and Kou et al. studied the effect of steady rotation about a vertical 
axis on freckle formation in Sn-Pb alloys. Sample & Hellawell considered solidification of 
the transparent alloy NH 4 CI-H 2 0 in a crucible rotating about an axis inclined between 0 and 
30 degrees with respect to the vertical. MOller (1990) and Weber et al. have studied the 
effect of rotation on the growth of Te-doped InSb crystals, and have shown that striations can 
be suppressed at sufficiently high rotation rates. These studies show that the degree of 
macrosegregation (freckles and striations) can be significantly reduced by rotation. 

There are two means by which steady rotation about a fixed axis can influence the motion 
of a fluid. In terms of a reference frame rotating with constant angular velocity Q about an 
axis, these correspond to the centripetal and Coriolis accelerations, for which the terms 
p L £lx£lxr and 2p L Qxu, respectively, are added to the momentum equation. Here, r is the 
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position vector measured from the axis of rotation, u is the local fluid velocity relative to the 
noninertial reference frame, and p L is the local density of the liquid. 

Previous theoretical studies of the effects of rotation on the onset of buoyancy-driven 
convection have been restricted to cases where no solidification occurs. For a horizontally 
unbounded layer of a single-component fluid, Chandrasekhar (1953), Chandrasekhar & 
Elbert (1955), and Niiler & Bisshopp (1965) have shown that steady uniform rotation 
about a vertical axis can significantly inhibit the onset of convection, with the Coriolis- 
related Taylor-Proudman mechanism (Chandrasekhar 1961) playing the dominant role. The 
effects of centripetal, Coriolis, and gravitational accelerations on convection in horizontally 
confined rotating fluids in cylindrical containers of various aspect ratios were considered in 
an early series of papers by Homsy & Hudson (1969, 1971a-b, 1972). More recently, for 
a single-component fluid Weber et at. (1990) have computed buoyancy-driven flows 
equivalent to those that can be driven by a temperature gradient maintained between the ends 
of an otherwise insulated right circular cylinder rotating at constant angular velocity about 
an axis perpendicular to and intersecting the cylinder axis but not passing through the 
cylinder. Three-dimensional computations, in which the variation of the magnitude of the 
centripetal acceleration OxQxr along the cylinder axis was neglected (a good approximation 
when the cylinder length is small compared to the shortest distance between the axis of 
rotation and the cylinder), were performed with the term accounting for the Coriolis 
acceleration 2flxu either included or omitted. Weber et al. found excellent agreement 
between experiment and computation when the Coriolis acceleration was included. 

For a binary fluid, Pearlstein (1981) has shown that the Coriolis acceleration can 
either stabilize or destabilize a horizontally unbounded layer, depending on the values of the 
Prandtl and Schmidt numbers, the dimensionless rotation rate (expressed in terms of a 
Taylor number), and the dimensionless temperature or solute gradient (expressed in terms 
of thermal and solutal Rayleigh numbers). Other work concerning the effect of rotation on 
doubly-diffusive convection in binary fluids (with no phase change) has been reported by 
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Sengupta S Gupta (1971), Masuda (1978), Antoranz & Velarde (1978, 1979), Schmitt & 
Umbert (1979), Riahi (1983). Worthem eta/. (1983). and Bhattacharjee (1988a-c). 

1.2. Overview of the Dissertation 

The purpose of the present work is to investigate the onset of morphological and 
buoyancy-driven convective instability in a binary liquid undergoing directional 
solidification and to identify the mechanism by which Coriolis effects affect the onset of 
buoyancy-driven convection. We consider the stability of plane-front solidification of a 
binary alloy (Pb-Sn) and the pseudobinary mercury cadmium telluride system, and the 

dendritic solidification of Pb-Sn. 

The dissertation is organized as follows. In Chapter 2. the convective and morphological 
instabilities in a binary liquid undergoing plane-front solidification by cooling from below 
are studied. The effects of rotation on the onset of buoyancy-driven convective and 
morphological instabilities are considered using a linear stability analysis. Results for dilute 
Pb-Sn alloys are presented and discussed. 

In Chapter 3, we use a linear stability analysis to study the onset of instability in binary 
Squids which exhibit a density maximum in the interior, and solidify with a nominally planar 
interlace. We consider the stability of plane-front solidification of the pseudobinary 
mercury cadmium telluride system (Hg t .,CdJe, where x is the bulk mole fraction of CdTe) 
from the melt. The effect of a steady uniform rotation about a vertical axis has also been 

investigated, as for the Pb-Sn system. 

In Chapter 4, we use a linear stability analysis to study the buoyancy-driven convection 
in the liquid and mushy zone and to investigate the effects of rotation on convective instability 
for alloys solidified dendritically by cooling from below. The mushy zone, consisting of 
dendrites saturated by interdendritic liquid, is modelled as a porous medium with anisotropic 
permeability. The local porosity, as well as the locations of the mushy-zone/solid and 
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liquid/mushy-zone interfaces are taken as dynamical variables, and so are determined as part 
of the solution. The analysis is illustrated by results for the Pb-Sn system. 

The notation is generally consistent throughout the dissertation excepting the usage of C 
and its subscripted variants as composition variables, which are defined when first used in 
each Chapter. The numerical techniques employed in each Chapter are generally similar, but 
differ in important details, and hence are explained in full In each case. 


CHAPTER 2 


Coriolis Effects on the Stability of Plane-front 
Solidification of Dilute Pb-Sn Binary Alloys 

2.1. Introduction 

During directional solidification of alloys, it is frequently desired to produce large single 
crystals with very low densities of macrosegregation defects and other imperfections. In 
principle, this can sometimes be achieved by -plane-front- solidification, in which the melt- 
solid interface remains perfectly planar. In such a case, the solidification process would be 
steady in a reference frame moving with the interface, and the only spatial variation would be 
in a direction normal to the interface. However, in real systems, several instabilities can 
cause departures from the nominally steady and one-dimensional plane-front case. 

Buoyancy-driven convection occurring during directional solidification of binary alloys 
leads to the formation of freckles and other macrosegregation defects deleterious to the 
mechanical properties of directionally solidified alloys. In particular, much attention is 
currently focused on the questions of when and how the freckles, channel segregates, and other 
macrosegregation defects are formed, with the main goal being to grow materials that are free 
of such macrosegregation defects resulting from convection during solidification. 

The possibility of using steady uniform rotation about a vertical axis to suppress the 
onset of buoyancy-driven convection during solidification of a binary alloy is considered 
using a linear stability analysis. The predicted inhibitory effects of rotation on convection 
are discussed in terms of previous experimental and theoretical studies of the effect of 
rotation on the onset of buoyancy-driven convection in single-component fluids heated from 
below. The plane-front solidification of Pb-Sn binary melts, for which the density depends 
monotonically on temperature and composition, is considered. 
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In this Chapter, we use linear stability analysis to study the effect of the Coriolis 
acceleration on convective and morphological instability for Pb-Sn alloys which solidify with 
a nominally planar interface. This work serves to identify the mechanism by which Coriolis 
effects affect the onset of convection in solidifying binary alloys. The analysis is illustrated 
by results for the Pb-Sn system. 

2.2. Basic State and Linear Disturbance Equations 

We adopt the model of solidification used by Coriell et al. (1980), in which the 
Oberbeck-Boussinesq equations govern motion in the liquid. In a reference frame translating 
with the nominally steady velocity (Vj) of the moving interface and rotating with constant 
angular velocity fl, as shown in figure 2.1 , the basic state 


U =(o,0,-f§-Vj), 
Plo 

(2.1) 


(2.2) 

^C.[t + JJSex P (-^Lz)]. 

(2.3) 

T S - T M + ~^ S + [l - exp( - z)] , 

(2.4) 


is the same used by Coriell et al., where u, T L , C^ andTg are, respectively, the basic state 
velocity (referred to the moving frame), temperature and solute concentration in the liquid, 
and temperature in the solid, k l and k s are the thermal diffusivities of the liquid and solid, 
respectively, D L is the diffusion coefficient, p s is the solid density. p L0 is the liquid density of 
pure lead at its melting point, C„ is the bulk concentration, G L is the liquid-side temperature 
gradient at the planar interface, Vj is the nominal growth rate, k is the segregation 
coefficient, m L is the slope of the liquidus, T M is the melting temperature of Pb, e« ps/plo- 1 
is the fractional shrinkage, Gs * (LVj + G L k L )/k s is the solid-side temperature gradient at the 
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planar interface, k L and k s are the thermal conductivities of the liquid and solid, 
respectively, and L is the latent heat of fusion per unit volume. (All concentrations used in this 
Chapter are weight percentages.) We note that the basic-state temperature distribution in 

the solid (2.4) is valid only near the interface. 

As in stability analyses of other flows subject to uniform rotation (Chandrasekhar 
1953; Chandrasekhar & Elbert 1955; Niiler & Bisshopp 1965; Pearlstein 1981; Nakagawa 
& Frenzen 1955), the linear disturbance equations we use differ from those for the 
nonrotating case only by addition of the Coriolis acceleration to the disturbance momentum 
equation. (This approach neglects the consequences of density variation in the term 
corresponding to the centripetal acceleration, and Is discussed in §2.6.) The equations 
governing small disturbances in the liquid are 

V*u=0, ( 2 - 5 ) 

3u _ v. — + 2Qxu = ttj- g TJ Z - t*c 9 CJ Z + v ^ 2u * (2.6) 

dt Plo oz Plo 

_ -H§- V, + G L w exp( — z ) = k l V 2 T l , (2.7) 

at Plo dz L Plo k l 

39l_ PS. Vl ^+G c wexp(--^-z) = D L V 2 C L , (2.8) 

at Plo «*z Plo d l 

where the dependent variables u, w, p, T L , and C L are, respectively, the disturbance values of 
the velocity, its z-component, pressure, temperature, and concentration in the liquid. Here. 
i z is the unit vector in the z-direction, g is the magnitude of the earth's gravitational 
acceleration, v is the viscosity, a T andac are the thermal and solutal expansion coefficients, 
respectively, G c = (k-DpsC.V^kpLoDO is the concentration gradient at the planar 
interface, Cl = Ci Q h is the angular velocity, and Q 0 *s a constant. The disturbance energy 



equation in the solid is 


(2.9) 
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where T s is the disturbance temperature in the solid. 

Taking the curl of (2.6) and the dot product of the result with l 2 , we obtain an equation 


Plo dz dz 


d). 


( 2 . 10 ) 


for the 2-component of the vorticity, © 2 . After twice taking the curl of (2.6), the vertical 
component of the resulting equation is 


^ V ^'S V 4 V2w + 2fi °| t=OT9V ^ + ac 9 v l C L + ' ,v4 w- 


( 2 . 11 ) 


The disturbance boundary conditions at the interface are 

u(x,y,0,t) = eVj-i frOW) , 
dx 

v(x,y,0,t) = eV, ^t x >y.t) 

3y 

w(x,y,0,t) = — c gT l ( x fY .t) t 
dt 

^>0,1 ) .-eV 1 vjrt ( x,y,t), 


(2.12a) 
(2.12b) 
(2.12c) 
(2.1 2d) 


n(x . y . 1) , (2 . 12e) 


V. 2 p S C, 


dz 

C 


+ V, C L (x,y,0,t) + ^ _ -PlqPl 3C L (x,y,0,t) 

d lPlo k at p s (1-k) dz ’ 11 ^ 


T L (x,y,0,t) + G l i)(x,y,t) = T s (x t y f 0,t) + G s n(x,y,t). 


(2.1 2g) 


T^x.y.O.t) + G l tj (x.y ,t) = m L G c Ti(x,y,t) + m L C L (x,y,0,t) + T M V vjr\ (x,y,t), (2. 1 2h) 
where i\ is the interface position and ¥ is the capillary coefficient. 
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2.3. Nondimenslonalization and Fourier Decomposition 

We scale the velocity, time, length, temperature, concentration, and vorticity with 
k l /H, H 2 /v, H, G l H, C„, and k l /H 2 , respectively, where H - D L /V, is the characteristic 
length (associated with the destabilizing gradient of rejected solute, see (2.3)), and write the 
horizontal and temporal dependence of the dimensionless disturbance quantities as exp(cn + 
ia X + ia v Y), where a x and a y are the x- and y-components of the wavevector. and a is the 
temporal eigenvalue. Substitution into (2.7)-(2.11) yields 

a(D 2 - a 2 ) W - D(D 2 - a 2 )W = - Ra T a 2 0 l - Ras L® ® Z c l + (° 2 “ ® 2 ) w_ ^o D£1 z >( 2 - 13 ) 

Sc 

oQ £ D q z DW + (D 2 -a 2 )fl z , (2- 14 ) 

Sc 

Pr CT 0 L - P Le D0 l = (D 2 - a 2 )e L - W exp(- pLe z) . (2.15) 

Sc a X L - p DX l =(D 2 -a 2 )X L + ^Wexp(-pZ) (2.16) 

in the melt, and 

Prc0 s -Le De s = K(D 2 -a 2 )0 s (2.17) 

in the solid. Here. Pr = v/k l is the Prandtl number, Sc = v/D L is the Schmidt number, 
Le = D|_/K|_ is the Lewis number, Ra T = ajgG L H 4 /(K L v) and Ra s - a c 9C«,H 3 /(D L v) are the 
thermal and solutal Rayleigh numbers, respectively, n* 0 = 2Q 0 H 2 /vis the square root of the 
Taylor number, p = p s /pLO fe the density ratio, k - k s /k l is the jhermal diffusivity ratio, 
and m (i_k)/k. The boundary conditions (2.12a-h) at the interface become 

W(0) = - e Pr o p, (2.18a) 

DW(0) = e Le a 2 p, (2.18b) 


fl z (0) = 0, 


(2.18c) 
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Sccrp = - k.L D©l( 0) + ks D©s(0) + Le (k[p - k^/ic - 1 /k) p, 
(1 - G) p = © s (0) - © L (0), 

(Sc ct + k p) p = - k Z L (0) - ^ DX l (0), 

(- a 2 Tm + m* - 1 ) p = ^ X L (0) + © L (0), 


(2.1 8d) 

(2.1 8e) 
(2.1 8f) 

(2.1 8g) 


where W, X|_, and © L are, respectively, the amplitudes of the disturbances to the vertical 
velocity, concentration, vertical vorticity, and temperature in the melt, © s is the amplitude 
of the disturbance temperature in the solid, and p (a constant) is the amplitude of the 
disturbed interface position. Here we define dimensionless parameters k^ ■ G L k L H/LD L , 
k| = G L k s H/LD L , m*= m L Gc/G L , = T M 4VG L H 2 , and G - G S /G L . For the far-field boundary 
conditions, we follow Coriell et at. (1980) and set all disturbances to zero 


W= DW = n z = X L = © L = 0 as Z -> oo. 


(2.19a) 


© s = 0 


as Z -> - oo, 


(2.19b) 


far from the interface. 


2.4. Numerical Solution 

Our objective is to find conditions under which infinitesimally small disturbances 
neither grow nor decay for a finite number of wavenumbers, and decay for all other 
wavenumbers. Disturbances which neither grow nor decay are said to be neutral. The 
neutral disturbances can be of two types, depending on the imaginary part of a. If the 
imaginary part of o is zero for a neutral disturbance, the onset of instability will be via 
monoton ically growing disturbances (steady onset). If the imaginary part, of a is not zero, the 
neutral disturbance will oscillate in time (oscillatory onset). 
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In the previous section, we formulated an eigenvalue problem for two systems of 
ordinary differential equations on two semi-infinite intervals (in the melt and in the solid), 
coupled by boundary conditions at the deformable interface. For convenience, we follow 
Coriell et al. (1980) and solve the problem on a finite interval [-h,h], where 2h is the 
dimensionless height of the computational domain. With regard to the far-field boundary 
conditions at Z = ± we set all disturbances to zero 

W= DW = £l z = X L =e L =0 at Z = h (2.20a) 

in the melt and 

6 S = 0 at Z = -h (2.20b) 

in the solid. We have used h - 10. (For steady onset, we have checked a number of our 
results using more accurate asymptotic boundary conditions applied at Z - ± h derived 
following Keller's (1976) procedure, and have found excellent agreement between the 
eigenvalues computed using the two sets of boundary conditions.) Since we use Chebyshev 
polynomials in our numerical solution, we scale the vertical coordinates in the liquid and 
solid regions by Z£ •= ( 2 Z - h)/h and z-j * ( 2Z + h)/h, respectively, so that each region lies 
between -1 and +1. The resulting system is then solved using a spectral Galerkin technique 
developed by Zebib (1987). The problem is thus reduced to a matrix eigenvalue problem 

As + oBs = 0, (2.21 ) 

where a is the temporal eigenvalue, and the elements of the square matrices A and B depend on 
a 2 _ a 2 + 3 y, the square of the horizontal wavenumber, the bulk concentration C.., and the 
other dimensionless parameters. The details are given In Appendix A. 

In what follows, we characterize the stability of the nominally plane-front solution in 
terms of the growth velocity V, and bulk concentration C., with all other parameters taken as 
fixed. For each value of Vj, we seek one or more critical values of C„ (denoted by C„ ) such 
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that for C„< Cl, disturbances decay for all wavenumbers, while for C_> Cl, disturbances 
grow for all wavenumbers in some range. In order to determine the critical value(s) of C„, 
we first determine the neutral curve ( C„ versus the wavenumber a) separating those 
combinations of C. and a for which all temporal eigenvalues a lie in the left half-plane (a 
stable basic state) from those for which at least one eigenvalue lies in the right half-plane 
(an unstable basic state). To determine the bulk concentration on the neutral curve for an 
arbitrarily chosen wavenumber a, we first guess a value of C„ and compute all eigenvalues a 
using (2.21). If all eigenvalues have negative real parts, the value of C„ is increased by 
doubling the previous value; otherwise the new value of C„ is chosen as half the previous 
value. This process is continued until we determine two values of C„ between which at least 
one sign change is obtained in the real part of the least stable temporal eigenvalue. The 
concentration on the neutral curve is then determined using a bisection method. 

To compute the critical concentration Cl at which instability first occurs (i.e., the 
minimum on the neutral curve), we arbitrarily choose a wavenumber and compute the 
corresponding C«, on the neutral curve using the procedure described above. We then fix C„ 
at the value computed at the previous step and compute the eigenvalues a for a discrete set of 
wavenumbers in a chosen range. From this set we select the wavenumber corresponding to 
the a with largest real part. If this wavenumber is at an endpoint of the chosen range, we 
extend the range to include the wavenumber corresponding to the largest Re(a). We then 
select this wavenumber and determine the corresponding C„ on the neutral curve. We 
continue this process until the relative change in C„ is less than 1 0" 6 . 

2.5, Results 

The solutions of (2.1 3)-(2.1 9) depend on fifteen dimensionless parameters, as defined 
in §§2.2 and 2.3. To determine the stability of a basic state with a nominally planar 
interface, numerical values of these parameters need to be specified. (Solutions on a finite 
interval depend also on h.) It is therefore not possible to numerically explore the effects of 
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more than a few combinations of these parameters on the onset of instability. In this work, 
we have thus restricted ourselves to the Pb-Sn system which, due to the low melting points of 
both components, has been the subject of several experimental studies. The parameter values 
(other than n* 0 ) are as used by Coriell et al. (1980) at reference conditions corresponding to 
pure lead at its melting point. (The diffusivity of Sn in Pb corresponds to an infinitely dilute 

solution.) 

To test our code, we first considered the nonrotating case, and compared our results to 
those of Coriell et al. Taking the liquid-side temperature gradient at the interface as 
<\ = 200 K cm- 1 , we computed neutral curves (C_ versus a) for various values of the 
solidification rate V,. The neutral curves were generally similar to those shown by Coriell et 
al. The only qualitative difference was that, in our work, several new oscillatory neutral 
curves were found to branch from steady neutral curves found by Coriell et al. and in the 
present work. We believe that our detection of these additional oscillatory neutral curves 
(which in each case lie well above the critical value of C„ and are hence of no practical 
consequence) is due to use of a numerical technique which simultaneously computes a large 
number of temporal eigenvalues at each combination of C„ and V,, as opposed to the shooting 
technique of Coriell et al., which individually computes the temporal eigenvalues by a one- 
point iteration scheme. The solid and dashes curves represent the steady and oscillatory onset 
of convection in each of the neutral curves. 

For G L = 200 K cm" 1 . n 0 - 0 rpm, and V, - 5 and 80 p sec- 1 , figure 2.2 shows neutral 
curves similar to those obtained by Coriell et al. For Vj-Spsec" 1 , steady onset of 
buoyancy-driven convective instability is found, and the critical concentration below which 
the plane-front solidification is stable for all wavenumbers occurs on the convective branch, 
as shown in figure 2.2(a). However, for ^ - 80 p sec- 1 , morphological and oscillatory 
convective neutral curves are obtained, and the critical concentration occurs on the 

shown in figure 2.2(b). Note that the critical 


morphological instability curve, as 


wavenumber for the morphological instability mode is larger than that for the buoyancy- 
driven convective mode. 

For this temperature gradient, a transition between these two extreme cases is observed 
for intermediate values of the solidification rate. Figure 2.3 shows a sequence of neutral 
curves for increasing V,, beginning with figure 2.3(a) for V, « 5 p sec -1 (figure 2.2a). 
When the solidification rate is increased to Vj * 12 5 psec -1 , the morphological neutral 
curve appears at much higher wavenumbers. For Vj - 15 p sec -1 , one of the buoyancy-driven 
convective branches pinches off and disappears (figure 2.3c). Another branch pinches off by 
Vj = 20 p sec -1 (figure 2.3d). As the solidification rate is further increased to 
V, - 30psec -1 , only one steady buoyancy-driven convective branch is left, and the critical 
concentration has shifted to larger values (figure 2.3e). The steady convective branch 
pinches off and becomes smaller with increasing solidification rate as shown in figure 
2.3(f-h) for V, - 35, 37.5, and 39.375 p sec -1 , respectively. Finally, for Vj * 40 psec -1 
(figure 2.3i) the steady convective branch disappears and the minimum value of C„ 
(corresponding to the critical condition) occurs on the morphological instability branch. For 
comparison, the extreme case V, - 80 p sec -1 (figure 2.2b) is shown again in figure 2.30). 

For G l « 200 K cm -1 and Vj - 30 p sec -1 , figure 2.4 shows the neutral curves (C„-a) 
for different rotation rates. Figure 2.4(a) shows the neutral curve for the nonrotating case 
(identical to figure 2.3e). If the system is rotated at Q 0 - 90 rpm, the critical concentration 
increases to higher values, as shown in figure 2.4(b). As the rotation rate is increased 
further to flo * "180 rpm, the steady convective branch pinches off, and becomes smaller for 
fl 0 - 270 rpm as illustrated in figure 2.4(c-d). Finally, the steady convective branch 
disappears, and the critical concentration occurs on the morphological branch for n 0 « 360 
rpm. The transition is similar to that shown in figure 2.3. 

We present our principal results in terms of stability boundaries in the Vj-Cl plane. 
For ease of comparison to the work of Coriell et at. (1980), results are presented in terms of 
dimensional variables. With Gj_ fixed at 200 K cm -1 , figure 2.5 shows stability boundaries 


for i2 0 - 0, 100, 200, 300, and 500 rpm. For each value of Q 0 . the stability boundary 
consists of some portion of the morphological branch (Cl decreasing with increasing Vj) found 
by Coriell et al. (their figure 1), joined to a convective branch. The critical value of C.on 
the convective branch is an increasing function of Qo 3t any value of Vj, clearly indicating the 
inhibitory effect of rotation on the onset of buoyancy-driven convection. We note that for 
Vj - 5 n sec -1 , rotation at £2 0 - 500 rpm Increases by slightly more than two orders of 
magnitude relative to the nonrotating case the critical Sn concentration above which the 
plane-front solution becomes unstable. We further note that the morphological branch is 
unaffected by rotation, whereas as Q 0 increases, the value of V, at which the onset of 
instability shifts from the convective branch to the morphological branch decreases from 
about 40 p sec -1 in the nonrotating case to about 27 p sec -1 for Q 0 - 500 rpm. 

For each value of n 0 . we see a local minimum near Vj - 1 psec -1 , with the minimum 
shifting to smaller growth velocities and becoming relatively more shallow as fl 0 increases. 
We note that the maximum relative stabilization by rotation occurs near the local minimum, 
and that for n 0 = 500 rpm, the critical bulk concentration of Sn is increased more than a 
hundredfold. Although Coriell et al. (1980) noted a local minimum in the stability boundary 
near Vj - 1 p sec -1 in the nonrotating case for the largest gravitational acceleration 
considered, they offered no explanation for its existence. This minimum is a consequence of 
the fact that as Vj -» 0, the concentration gradient G c = (k-1)p s C. o V 1 /(kp L0 D) vanishes. As the 
temperature gradient is independent of and is stabilizing, the critical value of must 
ultimately increase as Vj -> 0. (Of course, the weight percentage of the solute, C„, cannot 
exceed 100). As Cl must initially decrease with increasing V,. there must be a local 
minimum on the convective branch before Cl can increase to join up with the morphological 
branch of the stability boundary. This nonmonotonic dependence of Cl on Vj can be 
interpreted in terms of the existence of four critical values of Vj for certain values of C„ . 
(In addition to the three shown in figures 2.5 and 2.6, we note that for sufficiently large Vj , 
Cl ultimately increases on the morphological branch (Davis 1990).) 


As discussed for the nonrotating case by Coriell et al. (1980), at each value of Q 0 the 
onset of morphological instability occurs via a short wavelength (large wavenumber) 
Instability, while convective instability sets In via disturbances with relatively longer 
wavelengths. In the following section this point is discussed in the context of the mechanism 
by which rotation inhibits the onset of convection. 

For G l = 400 K cnr 1 , figure 2.6 shows stability boundaries (C* versus Vj), analogous to 
those for G L = 200 K cm -1 . Aside from a slight shift of the stability boundaries to higher 
values of Cl, the results are qualitatively similar to those for G L = 200 K cm -1 . In 
particular, the onset of buoyancy-driven convection is suppressed but the morphological 
instability is not influenced by rotation, and the critical value of C„ passes through a 
minimum near V,- 1 nsec -1 . 

2.6. Discussion 

The remarkable stabilization obtainable at low growth rates (more than a hundredfold 
increase in the critical value of C„ at Vj = 1 \i sec -1 can be achieved by rotating the layer at 
500 rpm for the two values of G L considered) is undoubtedly due to the well-known Taylor- 
Proudman mechanism, described by Chandrasekhar (1961). According to the Taylor- 
Proudman theorem, steady motion parallel to the axis of rotation in a uniformly rotating 
Inviscid fluid is prohibited at any nonzero rotation rate. If this theorem were strictly 
applicable to a viscous fluid, the onset of steady convection would be prohibited, since the flow 
in convection cells must have a vertical component. Instead, In a viscous fluid, one sees an 
inhibition of the onset of steady convection, with the degree of inhibition (expressed here as an 
increase in CJ increasing with £i 0 . That the onset of oscillatory convection is hardly affected 
is due to the fact that the Taylor-Proudman theorem applies only to steady flows. It is also not 
surprising that the morphological instability is unaffected by rotation. The morphological 
instability occurs at very short wavelengths, so the motion is almost perpendicular to the 
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solidification front (i.e., aligned with the axis of rotation). Hence, the Coriolis acceleration 

does not sensibly affect the morphological instability. 

For a horizontally unbounded binary fluid layer in which the density depends on 
temperature and one composition variable. Pearlstein’s linear stability analysis (Pearlstein 
1981) shows that Coriolis effects generally inhibit the onset of convection (by the Taylor- 
Proudman mechanism). Under some conditions, however, rotation can destabilize the layer, 
depending on the values of Pr, Sc, the dimensionless rotation rate (characterized by a Taylor 
number), and the dimensionless temperature or solute gradients (characterized by thermal 
and solutal Rayleigh numbers). For conditions under which destabilization (on a linear basis) 
occurs relative to the nonrotating case, instability sets in via an oscillatory mode, in which the 
natural frequency of oscillation of a buoyant fluid element is tuned (by rotation) in such a way 
that there is a local minimum in the critical value of Ra T as a function of the dimensionless 
rotation rate. This behavior was found (Pearlstein 1981) for Pr and Sc both less than unity, 
although there is no apparent reason why such destabilization cannot occur under other 

conditions when onset is via an oscillatory mode. 

Although Pearlstein (1981) found in the rotating doubly-diffusive case that 
for P r <1 <Sc (a condition satisfied in the present case, in which Pr - 0.023 and Sc - 81) 
there can exist as many as three critical values of the solute Rayleigh number for certain 
values of the Taylor, Prandtl, Schmidt, and thermal Rayleigh numbers, we have found no 
evidence of such multivalued stability boundaries in the present calculations. As in the case 
investigated earlier (Pearlstein 1981), it is possible that such behavior occurs in relatively 
small regions of the parameter space (G u , V„ etc.) and has gone undetected so far. As discussed 
in § 2.5, however, figures 2.5 and 2.6 imply that for certain values of C_ there exist four 
critical values of V, (including the unshown portion of the morphological branch). 

Even though the foregoing analysis is restricted to a horizontally unbounded fluid layer, 
the work of Homsy & Hudson (1971a) and Buhler & Oertel (1982) suggests that its 
predictions will be qualitatively correct for finite aspect ratios (ratio of mold radius to 
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height) if the parameter flj R 0 /g (a Froude number, where R 0 is the mold radius) is 
sufficiently small. For the onset of thermal convection in rotating water or mercury layers 
heated from below, the excellent quantitative agreement between classical linear stability 
analysis for a horizontally unbounded layer (Chandrasekhar 1953; Chandrasekhar & Elbert 
1955; Nakagawa & Frenzen 1955) and experimental work for finite aspect ratios (Nakagawa 
& Frenzen 1955; Fultz & Nakagawa 1955; Goroff 1960) provides a clear demonstration of 

the potential of the Coriolis acceleration to suppress buoyancy-driven convection in a 
rotating fluid. 

Interpretation of the results of a stability analysis restricted to infinitesimal 
disturbances is obviously subject to the caveat that larger disturbances might grow, even 
though sufficiently small disturbances are predicted to decay. Indeed, it is known that for 
rotating fluid layers heated from below, the onset of thermal convection sometimes does occur 
(Veronis 1959, 1966, 1968) at lower Rayleigh numbers than predicted by linear theory. 
However, in that case, accounting for finite (i.e., non-infinitesimal) amplitude disturbances 
modifies the quantitative predictions of the theory; the basic qualitative prediction of 
stabilization by rotation remains unchanged. 

The relatively modest rotation rates required to significantly inhibit the onset of 
convection in the Pb-Sn system make the proposed method an interesting candidate for a 
program of laboratory experiments. (We note here that the experimental work of Muller 
(1990) and Weber et at. (1990) focuses on the effect of rotation on the time-dependence of 
the supercritical flow, with no information given on the effect of rotation on the suppression 
of motion. Furthermore, although these authors have concluded that the Coriolis acceleration 
is key to the elimination of striations at high rotation rates, their experimental design 
complicates the separation of effects of the Coriolis acceleration from the increased “pseudo- 
gravitational effects associated with the centripetal acceleration.) Experiments might be 
conducted using a completely filled cylindrical mold, thermally insulated on the vertical 
surface, and mounted axisymmetrically on a rotating horizontal turntable. If the liquid at the 
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top of the rotating mold is in contact with a gas or vacuum, it will have a nearly paraboloidal 
free surface, on which the elevation above the point on the axis of rotation is 
Az-rf 0 r 2 /(2g). Although this configuration is consistent with rigid-body rotation, it also 
leads to a nearly paraboloidal solid-melt interface. For high rotation rates or large mold 
radius, this will in turn lead to significant radial variations in the solidified alloy. Also, 
cooling at the radial boundary leads to a radial temperature gradient, which in turn leads to 
significant centrifugal effects. For these reasons, the results of the present analysis cannot 
be compared to the experimental work of Kou (1978), Kou et al. (1978), or Sample & 
Hellawell (1984), in which strong radial variations in macrosegregation are observed in the 
solid. Experiments of the type proposed above were initiated some time ago by Copley 
(1976) for the crystallization of ammonium chloride from aqueous solution. 

Finally, we note that for many binary systems, plane-front solidification does not occur 
at practical growth velocities, and that the morphological instability results in dendritic 
solidification. In this case, rotation might also suppress buoyancy-driven convection in the 
melt and interdendritic liquid; this possibility is considered for the Pb-Sn system in 

Chapter 4 . 
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CHAPTER 3 

Stability of Plane-Front Solidification of a Binary Liquid for 
which the Density Depends Nonmonotonically on 
Temperature, including Coriolis Effects 

3.1. Introduction 

In the previous work on binary liquids (e.g., the study of lead-tin alloys by Coriell et at. 

1980) solidified by cooling from below (as in Bridgman growth), consideration was 
restricted to the case in which the liquid density (p) depended linearly on temperature (T) 
and solute mole fraction (C), i.e., 


P(T.C) = p 0 [l - a(T-T 0 ) - P(C-C 0 )J (3.1) 

where T 0 and C 0 are reference values. In general, the vertical variation of density is given by 


m dp 9T , 9p 3C 
dz dTdz dCdz 


(3.2) 


so that a local density extremum can occur even if the density depends monotonically on each 
of the stratifying agencies (as in (3.1)), and the profiles of the stratifying agencies vary 
monotonically with the vertical coordinate z . For example, when T increases linearly with 
increasing elevation above the liquid-solid interface (an excellent approximation for many 
systems; cf. §3.3) and solute is rejected from or preferentially Incorporated into the solid, 
the steady one-dimensional temperature and solute profiles near the interface can be 
approximated by 


and 


T(z)«T 0 + G l z 
C( z) • C_[l + !=£ e" Vl$/Dl -] 


(3.3a) 


(3.3b) 
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where G L > 0 is the liquid-side temperature gradient at the interface, V, is the solidification 
rate. D L is the solute diffusion coefficient, and k is the segregation coefficient. We see that 
dC/dz < 0 for k < 1 (rejection of solute) and dC :/dz > 0 for k > 1 (preferential incorporation 

of solute). From (3.2) and (3.3), we have 


oGl + 

Po dz 


PViC, 1-k e ~VjZ/D L ^ 
D|_ k 


(3.4) 


so that if the inequality 


PViC, 1-k r 1 
aG L D L k 


is satisfied, the density will assume a maximum value at 


z* 



cG(_D|_k ~ l 
MCJI-k)] 


(3.5) 


(3.6) 


in the liquid, notwithstanding the monotonic dependence of p on T and C. and of T and C on z. 

The density can also depend nonmonotonically on z if p(T,C) varies nonmonotonically 
with one of the stratifying agencies, say T. Among the binary systems in which density does 
not depend monotonically on temperature are many dilute aqueous solutions, for which the 
local density maximum is associated with the 3.98°C density maximum of pure water (at one 
atmosphere), and dilute solutions of cadmium telluride in mercury telluride. As noted in 
§§3.5 and 3.6, the density can also have more than one local extremum, even though the 

temperature and solute distributions depend monotonically on z. 

Although fluid layers having local density maxima have been the subjects of previous 
studies of the onset of buoyancy-driven convection in water (Veronis 1963; Merker et at. 
1979; Normand & Azouni 1992), directional solidification of a binary liquid differs in 
several important ways. First, one must account for the deformable moving interface. 
Second, in a fluid layer with a linear vertical temperature gradient, nonmonotonic dependence 
of density on temperature corresponds to a local density maximum in a single-component 
fluid, whereas (3.2) shows that in the binary case with density depending on temperature and 
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composition, the elevation where 3p/3T - 0 generally does not correspond to a local density 
maximum. Moreover, it is well-known that in a binary fluid, onset of buoyancy-driven 
convection depends on the individual stratifying agencies (temperature and solute in the 
present case), rather than on the density gradient. (Turner 1973, Chapter 8). Although this 
point has been recognized in analyses of convection in dilute aqueous solutions of NaCI (Foster 
1972; Gebhart & Mollendorf 1978; Qureshi & Gebhart 1986), the stability of plane-front 
solidification of a liquid with a nonmonotonic dependence of density on temperature has not 
been previously considered. 

The effects <5f rotation on convective and morphological instabilities have been reviewed 
by Oztekin & Pearlstein (1992) and in Chapter 2 herein. In the linear stablity analysis for 
Pb-Sn alloys, the results of Chapter 2 show that the Coriolis acceleration suppresses 
buoyancy-driven convection in the melt, and plane-front solidification becomes stable in a 
larger range of solidification rate and concentration for any given temperature gradient. 
Here, we investigate the effect of rotation on the onset of instability for Hg^CdxTe 
pseudobinary alloys. 

We use a linear stability analysis to study the onset of instability in binary liquids which 
exhibit a density maximum in the interior, and solidify with a nominally planar interface. 
We also assess the potential of uniform rotation to suppress the onset of buoyancy-driven 
convection. We consider the stability of plane-front solidification of the pseudobinary alloy 
mercury cadmium telluride (Hg^CdjJe, where x is the bulk mole fraction of CdTe, herein 
denoted by CJ from the melt. The growth of Hg^CdxTe crystals is of considerable practical 
interest because of the uses this material finds in the fabrication of infrared detectors and 
other electro-optical devices. For these applications, crystal size and defect density are 
critical, with the goal being to produce large single crystals with very high degrees of 
uniformity and very low densities of macrosegregation defects and other imperfections. This 
has led to a number of experimental (Galazka et al. 1981; Capper et al. 1986) and theoretical 
(Bourret et al. 1985; Kim & Brown 1989; Apanovich & Ljumkis 1991) studies of 
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solidification of Hg^CdJe from the melt. For a more detailed discussion, the reader is 
referred to the reviews by Micklethwaite (1981) and Capper (1989). 

This Chapter is organized as follows. In §3.2, we present the governing equations and an 
appropriate nondimensionalization. The one-dimensional basic state and linear disturbance 
equations are given in §3.3. The numerical solution technique is described in §3.4. Results 
for the Hg^xCdJe pseudobinary system are presented in §3.5a and §3.5b for the nonrotating 
and rotating cases, respectively, followed in §3.6 by a general discussion of the solidification 
of binary liquids in which the density depends nonmonotonically on temperature. 

3.2. Governing Equations and Formulation 
3.2. 1. Equation of State 

For the solidification of mercury cadmium telluride, we have used an equation of state 
based on the experimental data of Chandra & Holland (1983) and Mokrovskii & Regel 
(1952). Since the density of pure HgTe liquid does not depend monotonically on temperature, 
the density of Hg^CdxTe will also depend nonmonotonically on temperature for sufficiently 
small CdTe mole fractions. The variation of liquid density with temperature and composition 
is represented by bivariate polynomials of the form 

p,(T,C)-i X b.C'(T-1030K)l lor T<T ma „(C) (3.7a) 

i-0 j-0 

and 

pi (T,C) =£ X d ij C , (T-1030K)J for T > T max (C) (3.7b) 

l-o j=o 

where 

3 

ijd - X 'i c 1 

i-0 


(3.8) 
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Is a polynomial fit to the largest temperature (as a function of C) for which density data are 
available in the literature. The linear dependence of density on temperature (3.7b) for 
T >T max( c ) *s chosen to match p L and 9p L /3T at T max . Details are given in Appendix B. 

3.2.2. Dimensionless Governing Equations 

We adopt the model of directional solidification used by Coriell et ah (1980), In which 
the thermophysical properties in the liquid are taken to be constant, except for the density. 
The spatial and temporal variation of liquid density according to (3.7a,b) will be accounted 
for in the buoyancy term and neglected everywhere else. 

We scale the velocity, time, length, pressure, solute mole fraction, and the difference 
between the local and interface temperatures with V,. H/V Jf H, p L0 vV,/H, C„, and G L H , 
respectively, where H ■ D|_/Vj is the characteristic length associated with the solute gradient 
at the interface, Vj is the nominal solidification rate (with dimensions of velocity), k l is the 
thermal diffusivity of the liquid, D L is the binary diffusivity of CdTe in the liquid, v is the 
kinematic viscosity, C„ is the bulk value of the solute mole fraction, p L0 is the liquid density 
at the bulk value of the solute mole fraction and the corresponding liquidus temperature, and 
G l is the nominal liquid-side temperature gradient at the interface. The dimensionless 
equations governing the fluid motion are the Oberbeck-Boussinesq equations in a reference 
frame translating with the nominally steady velocity (i 2 V,) of the moving interface 



V . u = 0 

• 

(3.9a) 


u . Vu + 2D* x u L + Cl* x (£1* x r)^ =- 

Vp-M[k2li, 2+ V 2 U ( 

(3.9b) 

Plot Sc 


u ( 3 i + u.vr L ; 

) - ? 2 Tl . 

(3.9c) 


$ — VC l = V 2 C l . 


(3.9d) 
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where the dimensionless variables u, p, T L , and C L are. respectively, velocity, pressure, 
temperature, and solute mole fraction in the liquid, Sc - v/D L is the Schmidt number. 
Le» D L /K L is the Lewis number, y - V,/(gD L ) 1/3 is the dimensionless solidification rate, g 
is the magnitude of the gravitational acceleration, and Cl « ClH/Vi - Cl (D L /g ) /y is the 
dimensionless angular velocity. We note that in (3.9a-d), the velocity is referred to the 
moving frame, unlike the mixed formulations employed by Coriell et al. (1980) and 
subseqent authors in which the velocity in the laboratory frame appears in equations written 

in the moving frame. 

The energy equation in the solid is 

*. v «T Sl (3.10) 

V 9t dz / 

where T s is the dimensionless temperature in the solid, k*- k s /k l ,s the thermal diffusivity 

ratio, and k s is the thermal diffusivity of the solid. 

Dimensionless boundary conditions at the interface are derived from conservation of 

mass 

-e V s • n - p* l 2 • n = u • n, (3.11a) 

and the no-slip assumption 

(u + l 2 ) • t = 0 , (3.11b) 

where n and t are unit vectors normal and tangential to the interface, respectively. V s is the 
local interface velocity scaled by V„ p* * p s /p L o is the density ratio, Ps is the (constant) solid 
density, and e = p*-1 is the fractional shrinkage. The dimensionless energy balance at the 

interface is 


£l* (V s + I 2 ). n = -kLVT L . n + A VT S • n , 


(3.11c) 


where the right-hand side corresponds to the difference between the normal components of 
the heat flux vectors in the liquid and solid, and the left-hand side corresponds to the product 
of the latent heat (L ■ L(C^)/Lo, where L(C.J is the latent heat per unit volume and 
Lq-L(O)) and the normal component of the solidification velocity. Here 
r-G L kLo/[L 0 (gD L ) 1/3 ] is the dimensionless liquid-side temperature gradient at the 
interface, A - k s (C.J/k L0 is the thermal conductivity ratio, k£- k L (C„)/k L0 , k L0 is the 
thermal conductivity of pure HgTe liquid at its melting temperature, and k L (C„) and k s (C„) 

are the thermal conductivities of the liquid and solid, respectively. The dimensionless solute 
balance at the interface is 


p*(C L - C s ) (V s + i z ) • n=-VC L » n, (3.1 Id) 

where C s is the solute mole fraction in the solid. We also require the temperature to be 
continuous 

t l= t s. (3. lie) 

across the interface, the liquid and solid solute mole fractions at the interface to be related 
according to the binary phase diagram 

Cs=kC Lf (3. Ilf) 

where k is the segregation coefficient, and the temperature and solute mole fraction on the 
liquid side of the interface to be related by 

(3.1,8, 

where m£ - rn L C«k L0 /(D L L 0 ), ¥*- T M k L0 ^g 1/3 /(DL /3 L 0 ), T M is the melting temperature of 
HgTe, m L is the slope of the liquidus, R, and R 2 are the principal radii of curvature of the 
interface, and 'P is the capillary coefficient. The temperature- and composition-dependence 
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of most of the thermophysical properties (other than p L ) is given for Hg^Cd/Te in Appendix 
C. For the remaining properties, values measured at selected temperatures and compositions 
are taken as constant over the entire range. 

3.3. Basic State and Linear Disturbance Equations 

In the reference frame described in §3.2.2, the steady one-dimensional basic state 

specified by 


v s = 0 , 

(3.12a) 

u = - p* i 2 , 

(3.12b) 

Tl “ p-Le L 1 - exp < " z) ] ' 

(3.12c) 

Cl= 1 + Z exp(- p* z) , 

(3.1 2d) 


in the liquid, and 

- LVr) r _ exp(-tf z)l , <3.l2e) 

in the solid is, with one exception, the same as that used by Coriell et al. (1980). Here V s , 
u,T l , C L ,andT s are, respectively, the dimensionless basic state interface velocity, velocity, 
temperature, and solute distributions in the liquid, and temperature distribution in the solid, 
3 ^ s « (i-k)/k. Our basic state differs from that of Coriell et al. in that V s - 0 and u is not 
proportional to the shrinkage e, because in our formulation the Interface and fluid velocities 
are referred to the same (moving) inertial reference frame. Note that for p* = 1 and 
Le- 0(1 0 -3 ), (3.12c) can be approximated by a linear function of z (e.g., (3.3a)) from the 
interface out to fairly large values of z. 

To determine the conditions under which a disturbance grows, we write the interface 
velocity, liquid velocity, temperature, solute mole fraction, and pressure as 
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V s (x,y,T) =0 + Vs(x,y,T) . < 3 .13a) 

u(x,y,z,T) = u + u'(x,y,z,x) , (3.13b) 

T L (x,y,z,x) -TJz) + T[(x,y,z,x) . (3.13c) 

C L (x,y,z,x) -C L (z) I + CJx.y.z.x) , (3.13d) 

p(x,y,z,x) =p(z) + p'(x,y,z,x) , (3.1 3e) 

and the solid temperature as 


T s (x.y,z,x) =Ts(z)+ Ts(x,y,z,x) , 


(3.1 3 f ) 


where V s , u , T L , C L , and p' are, respectively, disturbances to the dimensionless interface 
velocity, liquid velocity, temperature, solute distribution, and pressure, T s is the 
dimensionless temperature disturbance in the solid, and p is the dimensionless basic state 


pressure distribution. Substituting (3.13a-f) into the governing equations and boundary 
conditions, subtracting the basic state equations, and retaining only linear terms, we obtain 
dimensionless disturbance equations and boundary conditions. The equations governing small 
disturbances in the liquid are 


V • u' = 0, 

/ r 

u [ * " p * it + w ' exp (~ P * U 2 )] = 

- P* ^ - Hp* w' exp(- p* z) = V 2 C l , 


(3.14a) 

(3.14b) 

(3.14c) 

(3.1 4d) 


where w is the dimensionless z-component of the disturbance velocity. Here, we define 
f T - f T LoD L /(p L0 k LO ). F c - C„f c /p L0 , f T - ap L (T L ,C L )/aT, f c - 3p L (f L ,C L )/3C, fl* = Q* i z , 
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and fljj = flo^i/Q 2 ) 1 ^? 2 ' s dimensionless rotation rate. We note that although the left- 
hand sides of the disturbance equations are similar to the corresponding left-hand sides of the 
disturbance equations of Coriell et al. and others, the second term in each of (3.14b-d) arises 
from a convective term of the type u .V( ), whereas in previous work these terms arose 
from the mixed nature of the formulation. 

The disturbance energy equation in the solid is 




(3.15) 


Taking the curl of (3.14b) and the dot product of the result with I 2 , we obtain an equation 


J_/9^1_p*^el-2n 0 — ^ = v2 »' 

Sc V dx p 3z 0 3z / 


(3.16a) 


for the z-component of the vorticity, © L . Taking the curl of (3.14b) twice and the dot 
product of the result with i z . we obtain an equation 

-L viTlx-S- vft+vV. <3.i6b> 

Sc \ di dz oZ f y *Sc T Sc 

for the z-component of the velocity, where V? « a 2 /3x 2 + 3 2 /3y 2 is the horizontal portion of 
the Laplacian operator. 

The disturbance boundary conditions linearized about the nominal interface position 


( z - 0) are 


u'(x,y,0,x) = e- 


v'(x,y,0,T) = e J 


w'(x.y.O.x) = -e- 


-c VXTt(x,y,x). 
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1 L » 9n(x,y,x) _ _ ■♦ aT L (x,y,0 > T) 
r L 3x ■ Kl dz + 


kp* n(x,y,x) + k C^x.y.O.x) + = _ J_ ^(x.y.O.x) 

ax So* ?z 


TJx.y.O.x) + Ti(x,y,x) = Ts(x,y,0,x) + + L*j?) xi(x,y,x), 

~ [ T L.(x.y.O,x) +Ti(x,y,x)] = [Ci(x,y t O,x) -p*S n(x,y,x)] + V^(x,y,x). 


(3.1 7e) 
(3.1 7f) 
(3.1 7g) 
(3.1 7h) 


9 

where V s (x,y,x) - I 2 aTj/ax is the linearization of the disturbance interface velocity, and tj is 
the local dimensionless interface deflection. 

We write the horizontal and temporal dependence of the dimensionless disturbance 
quantities as exp(cjT + ia x x + ia y y), where a x and a y are the x- and y-components of the 
horizontal wavevector, and a is the temporal eigenvalue. Substitution into (3.14c,d), 
(3.15), and (3.16) yields 


^ [ a(D 2 - a 2 )W - p*D(D 2 - a 2 )W + 2£2 0 * Da>] = 


■^ aSeL '^ a!xL+(D2 - a2)2w ' (3 - 18a) 

~ (aco-p*Dco + 2flo DW)= (D 2 -a 2 )©, (3.18b) 

Le(o0 L -p* D6 L )= (D 2 -a 2 )e L - Le Wexp(-p*Lez) , (3.18c) 

a X L - p* D3 C l = (D 2 - a 2 )z L + Sp* W exp(- p*z) (3 . 1 8d) 


in the liquid, and 


Le (a e s - De s ) = k* (D 2 - a 2 )e s 


(3.19) 
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in the solid; the boundary conditions (3.17a-h) at the interface become 


W(0) = - e o p. 

(3.20a) 

DW(0) = e a 2 p. 

(3.20b) 

e 

/«“s 

o 

II 

o 

(3.20c) 

T L* a p = - k£ D6l(0) + A D0 S (O) + Le (kj>* - P. 

(3.20d) 

(a + k p*)p = - k X L (0) - ^ DX L (0), 

(3.20e) 

0 L (O)+ p = e s (0)+^ (kj;+ L’^p, 

(3 . 2 Of ) 

[o L (0) + p] = [x L (0> - p’H p] - ^xa 2 p . 

(3.20g) 


where W, a>, X L , and e L are, respectively, the amplitudes of the disturbances to the vertical 
velocity, vertical vorticity, solute distribution, and temperature in the liquid. 6 S is the 
amplitude of the disturbance temperature in the solid, and p (a constant) is the amplitude of 
the disturbed interface position. For the boundary conditions far from the interface, we have 
followed Coriell et al. (1980) and set all disturbances to zero 

W = DW = co = X L = 9|_ = 0 asz->«o, (3.21 a-e) 

e s = 0 as z ->-«». (3.21 f) 

3.4. Numerical Solution 

Our objective is to find conditions for which infinitesimally small disturbances decay 
(Re (a) < 0) for all but a finite number of critical wavenumbers (typically one), and are 
neutral (Re(a) - 0) for the critical wavenumber(s). These conditions separate basic states 
that are linearly stable from those that are not. 
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In the previous section, we formulated an eigenvalue problem for systems of ordinary 
differential equations on two semi-infinite intervals (in the liquid and solid), coupled by 
boundary conditions at the deformable interface. For convenience, we follow Coriell at al. 
(1980) and solve the problem on a finite computational domain [-h,h]. With regard to the 
far-field boundary conditions at z = ± «>, we set all disturbances to zero at z = h 

W = DW = co = X L = e L = 0 at z = h , (3.22a-e) 

in the liquid and at z = -h 

e s = 0 at z = -h . (3 . 2 2 f ) 

in the solid, where we have taken h to be at least 10, depending on the vertical structure of 
the basic state temperature and solute fields given by (3.1 2c, d). (For steady onset, we have 
checked a number of our results using more accurate asymptotic boundary conditions applied 
at z “ ± h derived following Keller's (1976) procedure, and have found excellent agreement 
between eigenvalues computed using the two sets of boundary conditions.) Since we use 
Chebyshev polynomials in our numerical solution, we scale the vertical coordinate in the 
liquid and solid regions by z 2 - ( 2 z - h)/h and z, - ( 2 z + h)/h, respectively, so that each 

region lies between -1 and +1. The resulting system is solved using a spectral Galerkin 
technique developed by Zebib (1987). 

We approximate the highest derivatives of the amplitude of the disturbance velocity, 
solute distribution, and temperature in the liquid, and temperature in the solid, by truncated 
sums of Chebyshev polynomials of the form 

J 

w< W - XWs). (3.23a) 

j=o 

J 

“ (2) (Z2> - £ M jTj(Z 2 ), 
j=0 



(3.23b) 


36 


x[ 2) (Z 2 ) = X P j T i (Z 2 ) ' 
j=0 

J 

e[ 2) (z z ) = ^t,(2 2 ), 

j=0 
j 

e^ 2) (zi) = XZjTj^). 
i=o 

where Tj is the j-th Chebyshev polynomial and the coefficients Kj. M jt P jt Qj, and Zj are to be 
found. Representations of lower order derivatives can be found by integrating (3.23a-e) and 
using standard properties of Chebyshev polynomials. The procedure described by Zebib 
(1987) reduces the problem to a matrix eigenvalue problem 

As + oBs = 0, (3.24) 

where a is the temporal eigenvalue, and the elements of the square matrices A and B depend on 
the square of the horizontal wavenumber a 2 = a£ + a 2 , the bulk mole fraction C„, and the other 
dimensionless parameters. A more detailed description is given In Appendix A. 

In what follows, we characterize the stability of the nominally plane-front solution in 
terms of the bulk mole fraction and dimensionless solidification rate y, with all other 
parameters taken as fixed. As shown in §3.5, there is a critical value y c such that for y > y c , 
there is no range of stable C„ (i.e., for any C_ a disturbance at some wavenumber will grow), 
while for y < y c , disturbances of every wavenumber decay for some range of the bulk mole 
fraction. (See figure 3.1 for a schematic representation.) The upper and lower limits of this 
range are denoted by Cl and Cl* , respectively, and depend on y. These critical values of C.„ 
correspond to the extrema on the morphological and convective neutral curves (C„ versus a) 
separating those combinations of C_ and a for which all temporal eigenvalues o lie in the left 
half-plane (LHP) from those for which at least one eigenvalue lies in the right half-plane, in 
turn corresponding to stable and unstable basic states, respectively. 


(3.23c) 

(3.23d) 

(3.23e) 
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In order to determine the bulk mole fractions on the morphological and convective 
neutral curves separating stable and unstable basic states for an arbitrarily chosen 
wavenumber a, we first compute all eigenvalues a of (3.24) at each of N values of C„ (C^, 1 
£ n £ N) in the range [C™ in , 0.2], where the lower bound is typically 10 -4 and the upper 
bound is determined by the largest bulk mole fraction of CdTe for which we have data for the 
equation of state. We then attempt to determine a range of C„ for which all temporal 
eigenvalues are in the LHP. If one of the original values of C„ selected is stable (i.e., all 
temporal eigenvalues lie in the LHP), we then determine two intervals such that as C„ 
increases, in one a transition from instability to stability occurs, and in the other a 
transition from stability to instability occurs. Through these two intervals pass the 
convective and morphological neutral curves, respectively, as discussed in §5. If none of the 
original N values of C„ is stable, we choose the value (say, C^) for which the most unstable 
temporal eigenvalue has the smallest real part, and subdivide the interval [C2* 1) , di +1 *J until 
we either find a stable value of C„ (at which juncture we proceed to isolate the two intervals 
described above), or abandon the search when the real part of the least stable temporal 
eigenvalue and the difference between consecutive values of C„ supports the expectation that 
all intermediate values of C„, are unstable. Having found the intervals of C„ in which the 
transitions occur, we then compute the bulk mole fractions on the convective and 
morphological neutral curves using a bisection method. 

To compute the critical bulk mole fraction (i.e., the minimum on the morphological 
neutral curve), we arbitrarily choose a wavenumber and compute the corresponding C„ on 
the morphological neutral curve using the procedure described above. We then fix C„ at the 
value computed in the previous step, and compute the eigenvalues a for a discrete set of 
wavenumbers in a chosen range. We next select the wavenumber from among this set at 
which the o with largest real part was obtained. If this wavenumber is at an endpoint of the 
chosen range, we extend the range until the wavenumber corresponding to the largest Re(a) is 
inside the range. We then select this wavenumber and determine the corresponding C.. on the 
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neutral curve. We continue this process until the relative change In C„ is less than 10 .We 
compute C” following the same procedure. Spectral and domain convergence results are 
shown in Tables 1 and 2. 

3.5. Results 

3.5.1. Nonrotating case 

We begin by considering the basic state profiles of temperature, solute distribution, and 
density in the liquid. Equation (3.12c) shows that (for p*=1) the dimensionless basic state 
temperature distribution decays exponentially away from the interface on a lengthscale 1/Le, 
where 4 x 10 -3 < Le < 10 -2 for Hg.,. x Cd x Te, and so is essentially linear on the scale of the basic 
state concentration profile (3.1 2d). As discussed in §3.1, the liquid density of Hg^CdxTe 
does not depend monotonically on temperature for sufficiently small C. Hence, qualitatively 
different vertical density stratifications are obtained for fixed values of y and r, depending on 
C„. The nature of the density stratification and the thermal and solutal contributions thereto 
have profound consequences for the onset of convection. 

For y - 1.32 x 10 -4 and r - 8.2 x 10 -4 (corresponding to the dimensional solidification 
rate and temperature gradient Vj = 0.5 p/sec and Gl ■ 25 K/cm, respectively), figure 3.2(a) 
shows basic state density profiles for different values of C„. For C. - 0.1, the density varies 
monotonically with elevation in the liquid. For C„ ■ 0.045, the density variation is still 
monotonic, but there is a considerable reduction in the magnitude of the density gradient near 
the interface. For smaller values of (0.025 and 0.01), the density assumes a local 
maximum in the liquid. Figure 3.2(b) shows the vertical density stratifications for different 
values of C„ for y = 2.64 x 10" 4 and r * 8.2 x 10" 4 . The basic state density profiles are 
generally similar to those for y* 1.32x lO^andr- 8.2 x 10 - *, being monotonic for larger 
values of C„, and nonmonotonic for smaller values. We note, however, that for C„ - 0.025, 
the basic state density stratification has both a local minimum and a local maximum. (In 
figure 3.2, the approximate coincidence of the pl/Plo curves is only apparent, as can be seen 
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from figure 3.3, in which for y - 2.64 x 10 -4 and r « 8.2 x KT 4 , the derivative of p L /p L0 
with respect to C.. vanishes at a slightly different z near 1.7 for each value of C„ shown.) 

The solutions of (3.18)-(3.22) depend on sixteen dimensionless parameters, as defined 
in §§3.2 and 3.3. To determine the stability of a basic state with a nominally planar 
interface, numerical values of these parameters need to be specified. (The approximate 
solutions on the finite interval also depend on h.) It is therefore not feasible to numerically 
explore the effects of more than a few combinations of these parameters on the onset of 
instability. In this work, we present neutral curves (C_ versus a) for fixed values of y and 
r, amplitudes of the disturbances as a function of the vertical coordinate for selected values of 
C„, y, and r, and stability boundaries (critical values of the bulk mole fraction C„ as a 
function of y) for selected values of r. The thermophysical properties are taken as constants, 
with most being evaluated at C„ using the functional forms shown in Appendices B and C at the 
melting temperature for liquid properties, and at the freezing temperature for solid 
properties. We have used g * 9.80 m/sec 2 . 

Figure 3.4(a-c) shows neutral curves for r« 1.64 x10~ 3 (G L - 50 K/cm) and three 
values of y. The solid neutral curve at the top of each figure corresponds to morphological 
instability. The minimum on that curve (denoted by C*) determines one point on the stability 
boundary. The dashed neutral curve in each figure corresponds to an oscillatory buoyancy- 
driven convective instability mode. Just above that curve, all disturbance decay, whereas 
just below, oscillatory disturbances grow. The maximum on that curve determines another 
point (C**) on the stability boundary. In the unstable region below the oscillatory convective 
neutral curve, solid neutral curves correspond to steady convective modes. Since we have 
found no liquid density data for C„ > 0.2, computations were not performed for bulk mole 
fractions in excess of 0.2 (above the dashed-dot line in each of figures 3.4(a-c). 

Figure 3.4(a) shows the neutral curves for y « 2.64 x 10 -4 . The buoyancy-driven 
instability modes occur in a wavenumber range of approximately 0.01 to 1 0. The extrema of 
the morphological and oscillatory convective neutral curves occur at « 0.0989 and 
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C* - 0.0369, respectively. For C** < C„ < Cl , infinitesimal disturbances decay for all 
wavenumbers. Note that the critical bulk mole fraction Cl (below which there is no 
morphological instability) occurs at a larger wavenumber than C‘* (above which there is no 
convective instability). When y is increased to 5.28 x 10 -4 (Vj« 2 p/sec), C^, decreases to 
0.0464 and Cl* increases to 0.042, so that the range of C„ in which disturbances decay for all 
wavenumbers is smaller, as shown in figure 3.4(b). The wavenumber range in which 
buoyancy-driven convection occurs is shifted slightly to the left. For y - 7.92 x 10 -4 
(Vj - 3 p/sec), the extremum on the morphological curve lies below the extremum on the 
convective curve, as shown in figure 3.4(c), and there is no range of C„ for which 
disturbances decay for all wavenumbers. For r ■ 1.64 x 10~ 3 , we find a critical value of the 
dimensionless solidification rate y c - 5.61 x 10' 4 (V, - 2.12 p/sec) above which plane- 
front solidification is unstable for all values of C m . 

The basic state density, its gradient and the thermal and solutai contributions thereto, 
and the amplitudes of disturbances to temperature, solute distribution, and vertical velocity, 
are shown as functions of z in figure 3.5(a-d) for y - 2.64 x 10 -4 , r- 1.64 x 10 -3 , and four 
values of C., at a - 2.85. Figures 3.5a(i), 3.5b(i), 3.5c(i), and 3.5d(i) show the basic state 
density profiles for increasing values of C„. We note that below C„ - 0.035, the density 
depends nonmonotonically on z. We point out that for the linear equation of state (3.1), the 
stabilizing or destabilizing character of the thermal and solutai stratifications can be 
determined by inspection of the basic state temperature and concentration profiles. For more 
complicated equations of state, it is necessary to individually examine the two terms on the 
right-hand side of (3.2). Thus, in figures 3.5a(ii), 3.5b(ii), 3.5c(ii), and 3.5d(ii), the 
solid, dashed, and dotted curves are the dimensionless basic state density gradient, and 
contributions to it from the solute and temperature gradients, respectively. Note that the 
solute gradient is stabilizing [(3p L /dC)(9C L /3z) < 0) throughout the liquid layer for all 
values of the bulk mole fraction C„. However, the temperature gradient is destabilizing 
l(3p L /3T)(3T L /3z) > 0] in a sublayer adjacent to the liquid-solid interface, and stabilizing 
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[(3p L /3T)(3T L /3z) < 0] above this sublayer, as shown in figures 3.5a(ii), 3.5b(ii), 
3.5c(ii), and 3.5d(ii). As C„, increases, the height of this sublayer decreases, the magnitude 
of the destabilizing thermal contribution (3p L /3T)(3T L /3z) to the density gradient decreases, 
and the magnitude of the stabilizing solute contribution (3p L /dC)(dC L /3z) increases. 

Figures 3.5a(iii-v), 3.5b(iii-v), 3.5c(iii-v), and 3.5d(iii-v) show the z-dependence 
of the real and imaginary parts (denoted by solid and doited curves respectively) of the 
disturbances to the temperature, solute distribution, and vertical velocity associated with the 
most unstable oscillatory convective instability mode for different values of C„. For 
C.. = 0.01 (a value that for a = 2.85 lies well below the oscillatory neutral curve in figure 
3.4a), figure 3.5a(i) shows that the basic state density profile density assumes a maximum 
near z - 2.7, and that between the liquid-solid interface and this maximum, the liquid layer 
is hydrostatically unstable. Also, in this lower part of the liquid the temperature and solute 
gradients are destabilizing and stabilizing, respectively. As is clearly seen in figure 
3.5a(iii-v), the disturbances are confined largely within the sublayer where the 
temperature gradient is destabilizing, although they penetrate weakly into the liquid above 
the density maximum, in which both the temperature and solute gradients are stabilizing. 

For C„ ■ 0.035 and 0.0369 figures 3.5b(i) and 3.5c(i) show that although there is no 
region in which the overall density gradient is hydrostatically unstable, there is a sublayer 
adjacent to the interface in which the temperature gradient is destabilizing. The disturbances 
associated with the oscillatory mode of convection are confined to this sublayer, as shown in 
figures 3.5b(iii-v) and 3.5c(iii-v). In this sublayer the temperature gradient is 
destabilizing and the solute gradient is stabilizing. As discussed in §3.6, these conditions in 
the sublayer thus correspond to the diffusive regime of doubly-diffusive convection, with 
oscillatory onset of buoyancy-driven convection. The vertical structure of the disturbances 
is very similar to that shown in figure 3.5a(iii-v). In both cases, the disturbances are 
confined largely within a sublayer in which the thermal stratification is destabilizing, and 
penetrate weakly into the liquid above the sublayer. The disturbances for C„ « 0.01, 
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however, show more structure in the sublayer than those for C„ - 0.035 or 0.0369. This is 
because the effect of the destabilizing temperature gradient is stronger than for C. - 0.035 
and 0.0369, as discussed above. For C„ - 0.07, however, the structure of the eigenfunctions 
is completely different, as shown in figure 3.5d(iii-v). In this case, the eigenfunctions are 
associated with the least stable morphological instability mode. Although there still exists a 
small sublayer in which the temperature gradient is destabilizing, the disturbances are not 
confined to that region. The disturbances decay strongly with Increasing z. 

We present our principal results in terms of stability boundaries in the y-C„ plane for 
four values of r in figure 3.6(a-d). The stability boundary consists of some portion of the 
morphological branch joined to a convective branch. The dashed and solid curves in each 
figure correspond to minima on the morphological neutral curve (Cl) and maxima on the 
neutral curve associated with buoyancy-driven convection (Cl*), respectively. Figure 
3.6(a) shows the stability boundary for r - 8.2 x 10" 4 . Above the morphological branch, 
disturbances grow in some wavenumber range and the plane-front solution is unstable. 
Similarly, disturbances grow below the convective branch for some range of a. On the other 
hand, below the morphological branch and above the convective branch, disturbances of all 
wavenumbers decay and plane-front solidification is stable. Beyond the critical solidification 
rate y c ■ 2.70 x 1 0 -4 (Vj - 1 .02 p/sec) at which the convective and morphological branches 
intersect, there is no stable range of bulk mole fraction C„. For r= 1.64 xIO -3 , figure 
3.6(b) shows the stable region in the plane analogous to that for r = 8.2 x 10 -4 . The 
morphological branch has shifted to the right and the convective branch has shifted slightly 
downward. Their intersection occurs at a higher y(y c - 5.61 x 10“ 4 , Vj - 2.12 p/sec), and 
the stable region is larger. Increasing r still more leads to further enlargement of the stable 
region and larger values of y c , as shown in figures 3.6(c) and 3.6(d) for r = 3.28 x 10~ 3 and 
6.56 x 10 -3 (G l « 100 and 200 K/cm), respectively. 
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3.5.2. Rotating case 

For r = 1 .64 x 1 CT 3 (G l - 50 K/cm), figure 3.7 shows stability boundaries for Q 0 « 0, 
25, 50, and 100 rpm. Above the morphological branch, disturbances grow for some range of 
wavenumber a and plane-front solidification is unstable. Similarly, below the convective 
branch appropriate to each rotation rate shown, disturbances grow for some range of a. 
However, between the morphological branch and the convective branch associated with each 
rotation rate (i.e., for C" < C.. < C^.), disturbances decay for all wavenumbers and plane- 
front solidification is stable. Note that the morphological branch is unaffected by rotation, 
whereas as Ci 0 increases, the convective branch is shifted downward. The critical CdTe bulk 
mole fraction C" on the convective branch is a decreasing function of CIq at any dimensionless 
solidification rate, clearly indicating the inhibitory effect of rotation on the onset of 
buoyancy-driven convection. Beyond the critical solidification rate (denoted by y c ) at which 
the convective and morphological branches intersect, there is no stable range of CdTe bulk 
mole fraction C„. Thus, plane-front solidification is unstable at all CdTe bulk mole fractions 
for sufficiently high solidification rates. The critical solidification rate y c at which the 
morphological and convective branches intersect occurs at higher y with increasing rotation 
rate. The critical value of y c is increased by more than a factor of ten at G 0 - 100 rpm 
relative to the nonrotating case. We also note that the convective branch has a relatively 
shallow local minimum near y * 2.64 x 10~*, the location of which is only weakly dependent 
on il 0 . Thus, for a given rotation rate, operation at the solidification rate corresponding to 
this local minimum allows plane-front solidification to be stably conducted at the lowest 
value of C„. 

For r = 8.2 x 10 (Gl « 25 K/cm) and Qq ■ 0 and 100 rpm, figure 3.8 shows stability 
boundaries in the y-C_ plane analogous to those for r = 1.64 x 10" 3 . As expected on the basis 
of the results of Coriell et at. (1980) for Pb-Sn without rotation, the morphological branch 
is shifted to the left, reducing the range of solidification rates and bulk mole fractions for 
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which plane-front solidification is stable. On the other hand, the convective branch is shifted 
downwards relative to the r = 1.64 x 10" 3 case, corresponding to an increase in the range of 
stable operating conditions. We note that decreasing the temperature gradient has very little 
effect on the onset of convection in a nonrotating layer, but reduces the range of stable CdTe 
bulk mole fractions by more than twofold for a layer rotating at 100 rpm. The critical 
solidification rate y c is more than a factor of five higher at Q 0 ■ 100 r P m than in the 
nonrotating case, although the stabilization is less than the factor of ten predicted at 
r = 1 .64 x 10 -3 . Although the reduction in r has very little effect on the critical value of y at 
which the morphological and convective stability boundaries intersect in the nonrotating case, 
for fi 0 * 100 rpm the critical value of y decreases by a factor of three when the temperature 
gradient is reduced. 

As discussed for the Pb-Sn case (Chapter 2), at each value of fto i^e onset of 
morphological instability occurs via a short wavelength (large wavenumber) instability, 
while convective instability sets in via disturbances with relatively longer wavelengths. In 
the following section this point is discussed in the context of the mechanism by which rotation 
Inhibits the onset of convection. 

3.6. Discussion 

The qualitative differences between the neutral curves and stability boundaries presented 
in §3.5 and those characteristic of the normal case in which density varies linearly with 
temperature and solute mole fraction (e.g., the Pb-Sn results of Chapter 2, cf. Coriell at al. 
1980, Oztekin & Pearlstein 1992) lead us to consider how the equation of state affects or 
alters the mechanism by which the onset of motion occurs during directional solidification of 
a binary liquid co ole d from below. 

We begin by noting that the mechanism responsible for the morphological instability is 
insensitive to the net density stratification in the liquid, and is insensitive to variations of 
temperature and composition outside a relatively thin layer adjacent to the interface. This is 
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evident from a comparison of the neutral curves and stability boundaries for the present case 
to those for Pb-Sn. 

On the other hand, the onset of buoyancy-driven convection depends strongly on the 
nature of the stabilizing or destabilizing temperature and solute gradients in the liquid, as 
well as on the diffusivities of the stratifying agencies (Stern 1960; Baines & Gill 1969). As 
discussed in §3.1 , the density of a binary liquid solidified by cooling from below can depend 
nonmonoton ically on the vertical coordinate even if the density depends monotonically on 
temperature and composition. For example, if the equation of state is of the form (3.1) and 
a >0, p >0, and k < 1 (as for Pb-Sn), then (3.5) provides a necessary criterion for the one- 
dimensional basic state to have a density maximum (3.6) within the liquid. However, for any 
system (e.g., Pb-Sn) characterized by the linear equation of state (3.1), even if the density 
has a maximum in the liquid, the temperature and composition distributions (3.3a) and 
(3.3b) will be stabilizing and destabilizing, respectively. Following the standard "parcel- 
argument for a linearly stratified doubly-diffusive fluid (Turner 1973, page 251), we note 
that if a fluid element is displaced downward, it will lose its excess heat more rapidly than its 
composition will adjust to the relatively solute-enriched surroundings, because the thermal 
diffusivity is much larger than the solute diffusivity. Hence, the displaced parcel will be 
dilute relative to the surroundings, and if p > 0 and the resulting buoyancy force is sufficient 
to maintain the downward motion, the one-dimensional basic state will be unstable. In the 
context of doubly-diffusive flows, this configuration is in the "fingering" regime, and the 
onset of convection is said to be "steady" or "monotonic". 

For binary liquids such as Hg-,. x Cd x Te, for which the dependence of density on 
temperature or solute mole fraction is not monotonic, solidification by cooling from below for 
some combinations of the bulk mole fraction, liquid-side temperature gradient, and 
solidification rate can lead to, as discussed in §3.1, a configuration in which there exists 
adjacent to the interface a sublayer in which the positive temperature and solute gradients 
are destabilizing and stabilizing, respectively. If a liquid parcel in this sublayer is displaced 
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downward, it will loss its excess heat more rapidly than its excess CdTe mole fraction, again 
because the thermal diffusivity is much Jarger than the solute diffusivity. Its CdTe mole 
fraction will thus exceed that of the surroundings, and so the parcel will be less dense than 
the surrounding liquid. The resulting buoyancy force will tend to move the fluid parcel 
vertically upward. Hence, the initial displacement engenders a restoring force, which can 
result in overshoot of the parcel's initial (equilibrium) position. Such a configuration is in 
the "diffusive" regime, and the motion is said to be "overstable", with the temporally growing 
oscillatory parcel displacements leading to the oscillatory onset of convection. 

The difference between the mechanisms by which the onset of motion occurs in 
configurations in which the temperature gradient is everywhere stabilizing, and in those in 
which temperature is destabilizing in a sublayer adjacent to the interface, is manifested in 
qualitative differences between the neutral curves, and hence the stability boundaries, for the 
Pb-Sn and Hg.,_ x Cd x Te systems. 

The most notable feature of the neutral curves in the present case is that the onset of 
buoyancy-driven convective instability is associated with an oscillatory mode, unlike the 
normal (e.g., Pb-Sn) case. Moreover, for each wavenumber for which it exists, the 
corresponding oscillatory neutral curve bounds the range of stable C„ from below, rather 
than from above as in the normal case. As the extremal values of C. on the morphological and 
oscillatory convective neutral curves approach, the Tange of bulk mole fractions for which 
plane-front solidification is stable disappears. This contrasts to the normal case, in which 
plane-front solidification is always stable for sufficiently dilute solutions (i.e., C„ lying 
below the minima of the morphological and steady convective neutral curves). 

These features of the neutral curves have profound consequences for the stability 
boundaries in the y-C. plane. For Hg 1 . x Cd x Te, plane-front solidification can be unstable at all 
solidification rates if the bulk mole fraction is chosen sufficiently low, as shown in §3.5. 
This is because the effect of a destabilizing temperature gradient in the sublayer adjacent to 
the interface becomes stronger as the bulk mole fraction decreases (due to a more prominent 
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density maximum as the composition approaches pure HgTe), and the stabilizing effect of the 
solute gradient diminishes (because it is proportional to CJ. On the other hand, in Pb-Sn 
directionally solidified by cooling from below, plane-front solidification is stable for any 
solidification rate and liquid-side temperature gradient if the solute bulk mole fraction Is 
sufficiently low (Coriell et al. 1980; Oztekin & Pearlstein 1992). 

For Hg^CdxTe, plane-front solidification can be unstable at all values of the bulk mole 
fraction (less than C„ « 0.2) if the solidification rate exceeds a critical value y c . The 
explanation for this result is made clear by reference to figure 3.6(a-d). For values of C„ 
below the dashed curves, plane-front solidification is unstable with respect to the onset of 
overstable buoyancy-driven convection, because the stabilizing influence of the bottom- 
heavy solute stratification is insufficiently strong to overcome the destabilizing temperature 
gradient in the sublayer adjacent to the interface. Above the solid curve, the configuration is 
unstable with respect to the morphological instability. As y approaches y c from below, the 
stable region is ’pinched" from below and above until at y c it finally disappears. This 
situation differs from that for Pb-Sn, for which plane-front solidification is stable in a 

region bounded above by the morphological and convective portions of the stability 
boundaries. 

In a single-component fluid with an unstably stratified layer overlying or underlying a 
stably stratified region, convective motion is typically localized in the unstably stratified 
layer, but may penetrate into the adjacent stable region. In his early analysis of the onset of 
thermal convection in a water layer with the temperature maintained at 0°C at the bottom and 
in excess of 4°C at the top, Veronis (1963) showed that convection occurs in the unstably 
stratified region and penetrates into the stably stratified region as well. He determined the 
extent of penetration for different ratios of the stably stratified layer thickness to the total 
layer thickness. He also found that convection in the stably stratified region is viscously 
coupled to the more vigorous motion in the unstable layer. Walton (1982) and Zangrando & 
Bertram (1985) have considered a doubly-diffusive fluid layer with a uniform vertical 
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temperature gradient and a nonuniform vertical solute gradient. Both of these studies show 
the disturbances to be localized about the neighborhood of the critical depth at which the 
salinity gradient reaches its minimum value. In our work, when the liquid-solid interface is 
adjacent to a sublayer in which the temperature gradient is destabilizing, the disturbances 
associated with the oscillatory onset of buoyancy-driven convection are largely confined to 
that thermally unstably stratified sublayer. 

As discussed above, the oscillatory onset of buoyancy-driven convection and the 
localization of the disturbances are associated with the occurrence of a sublayer in which the 
temperature gradient is destabilizing. Hence, the results should be expected to apply 
qualitatively to the solidification of other binary or multicomponent liquids in which 
(3p L /aT)(aT L /3z) changes sign within the layer. 

We recall from §3.5 (see figure 3.2b for C„ ■ 0.025) that for Hgi. x Cd x Te, basic state 
density profiles more complicated than monotonic and unimodal can occur. For several cases 
(both unstable and neutrally stable) in which the basic state density profile exhibits both a 
local maximum and a local minimum, examination of the eigenfunctions of the disturbances to 
the temperature, solute, and vertical velocity shows that the disturbances are localized 
within the sublayer in which the thermal stratification is destabilizing. The eigenfunctions 
are qualitatively similar to those for the case in which a single extremum (a maximum) 
occurs, which is not surprising since in each case disturbance localization is a consequence of 
confinement of the destabilizing temperature gradient to a sublayer adjacent to the liquid- 
solid interface, and even when the basic state density profile has two local extrema, the 
thermal contribution to the right-hand side of (3.2) has only one sign change. 

Our results demonstrate that nonmonotonic variation of density with temperature can 
have dramatic qualitative effects on the onset of Instability in an unbounded horizontal fluid 
layer undergoing plane-front solidification. The existence of a sublayer in which the thermal 
stratification is destabilizing should also have important consequences for the convection 
which occurs in vertical Bridgman growth in ampoules of finite radius, as well as in zone 
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melting and other processes used to grow Hg l0( Cd x Te crystals from the melt. To date, 
however, these effects have not been observed in simulations of buoyancy-driven convection 
in Hg^CdxTe for these geometries, which have used equations of state that do not properly 
account for the nonmonotonic dependence of density on temperature. The axisymmetric 
simulations of vertical Bridgman growth by Kim & Brown (1989) employ a linear equation 
of state of the form (3.1), even though the CdTe mole fraction in the liquid near the interface 
Is about 0.075 in each computation and the density is known to depend nonmonotonically on 
temperature for CdTe mole fractions up to at least 0.1. The computations of Apanovich & 
Ljumkis (1991) for zone melting of Hg t . x Cd x Te use an unusual equation of state in which the 
density of Hg 1 . x Cd x Te depends quadratically on temperature for x *0.13, and is independent of 
temperature for x - 0.13. 

For Pb-Sn alloys, we have shown in Chapter 2 (cf. Oztekin and Pearlstein 1992) that 
the onset of convection can be suppressed significantly at modest rotation rates. This is 
undoubtedly due to the well-known Taylor-Proudman mechanism, described by 
Chandrasekhar (1961). According to the Taylor-Proudman theorem, steady motion parallel 
to the axis of rotation in a uniformly rotating inviscid fluid is prohibited at any nonzero 
rotation rate. If this theorem were strictly applicable to a viscous fluid, the onset of steady 
convection would be prohibited, since the flow in convection cells must have a vertical 
component. Instead, in a viscous fluid, one sees an inhibition of the onset of steady convection, 
with the degree of inhibition (expressed here as an Increase in C J increasing with fl 0 . 
Although, the Taylor-Proudman theorem applies only to steady flows, we find that the 
rotation still suppress the oscillatory onset of buoyancy-driven convection during directional 
solidification of Hg^Cd^e, as presented in §3.6. It is not surprising that the morphological 
instability is unaffected by rotation. The morphological instability occurs at very short 
wavelengths, so the motion is almost perpendicular to the solidification front (i.e., aligned 
with the axis of rotation). Hence, the Coriolis acceleration does not sensibly affect the 
morphological instability. 
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Finally, Antar (1991) has presented a linear analysis of the onset of convection in a 
horizontal fluid layer cooled from below, using an approximate equation of state for 
Hg! x Cd x Te. Although he purports to consider "convective Instabilities in the melt for 
solidifying mercury cadmium telluride" his analysis differs from ours in that no phenomena 
associated with solidification (existence of a moving or deformable interface at which phase 
change occurs, nonlinear basic state solute stratification due to rejection or preferential 
incorporation at the interface, latent heat effects at the interface, etc.) are included. 
Moreover, comparison to Antar's results is precluded because his linear vertical basic state 
thermal stratification is inconsistent with his definition of the temperature difference across 
the fluid layer. Specifically, the basic state temperature distribution (Antar's equation (1)) 
requires the bottom temperature to be T 0 -AT, which is inconsistent with his definitions of 
T o (temperature where 3p/3T is zero, determined by Antar’s equation of state) and AT 
(temperature difference between top and bottom of the layer). This error vitiates his 
analysis, which uses a thermal Rayleigh number defined in terms of a temperature difference 
inconsistent with the basic state. That this inconsistency is not simply an apparent one due to 
a typographical error is confirmed by reference to an earlier paper (Antar 1987) employing 
the same basic state and disturbance equations, in which a different (and highly unusual) 
definition of AT (related to the actual overall temperature difference by a constant dependent 
on the temperature at the top or bottom wall; see figure 1 of Antar 1987) was used to achieve 
consistency. That definition of AT was the only one consistent with the basic state, and hence 
with the analysis and results. Unfortunately, Antar's 1991 paper defines AT as the overall 
temperature difference, and so is wrong. 
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CHAPTER 4 

Stability of Dendritic Solidification of a Binary Liquid, 

including Coriolis Effects 

4.1. Introduction 

Although plane-front solidification (in which the liquid and solid are separated by a 
well-defined interface that remains planar and translates at a constant velocity) is desirable 
in many applications, solidification of many binary alloys under conditions of practical 

interest occurs dendritically. Dendritic solidification occurs In a so-called "mushy" zone 

% 

between the solid and melt in which dendrites, consisting of long vertical stems with 
secondary and tertiary branches, are immersed in interdendritic fluid. This phenomenon is 
illustrated in photographs by Huppert & Worster (1985) and Chen & Chen (1991) of the 
mushy zone in a dendritically solidifying aqueous ammonium chloride solution. In such a 
case, the dendritic region might be modeled as a single domain, in which the macroscopic 
structure (i.e., features that scale with the mushy-zone thickness) are considered to be more 
important than the microscopic morphology (i.e., features dependent on a detailed stochastic 
description of dendrite location and shape). Thus, in many previous theoretical investigations 
of its evolution and behavior, the mushy zone has been considered as a porous medium with 
anisotropic permeability. A review of plane-front and dendritic solidification has been given 
by Glicksman et al. (1 986). 

The gradients of temperature and solute which necessarily occur in a solidifying binary 
or multicomponent liquid can destabilize the nominally motionless basic state, leading to the 
onset of buoyancy-driven convection. Analysis of the convective motion in solidifying alloys 
is complicated by the fact that temperature and composition in the mushy zone are related by 
the phase diagram, and by the fact that phase change leads to motion and deformation of the 
boundaries between the melt, mushy zone, and solid. We will refer to these boundaries as the 
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liquid/mushy-zone and mushy-zone/solid interfaces. Characteristic of a moving boundary 
problem, their locations are not known in advance, and must be computed as part of the 
solution. 

For binary liquids with k < 1 (e.g.. Pb-Sn alloys) for which the equation of state is 
adequately represented by the linear relation 

P = Po[1-«t( t - t o)-“c( c - c o)] - ^ 4 - 1) 

cooling from below leads to dendritic solidification with a stabilizing temperature gradient 
and a destabilizing composition gradient in both the mushy zone and liquid. The vertical 
composition gradient in the mushy zone is almost constant, while in the liquid there is a thin 
compositional boundary layer adjacent to the liquid/mushy-zone interface. These 
stratifications can lead to two different convective modes, as discussed by Worster (1992) in 
his analysis of the stability of dendritically solidified binary alloys. One mode, referred to as 
the boundary layer mode, associated with the thin compositional boundary layer in the liquid 
adjacent to the liquid/mushy-zone interface, is largely confined within this layer, and has a 
relatively small critical wavelength (on the order of its thickness). The boundary layer mode 
penetrates only weakly into the mushy zone, which is thus essentially stagnant. On the other 
hand, the mushy layer mode is driven by the adverse compositional stratification in the 
mushy zone and has a larger critical wavelength that scales with the mushy-zone thickness. 
These two convective mode types are analogous to those found by Chen & Chen (1988) in 
superposed porous and liquid layers heated from below. ~ 7t - " 

The onset of buoyancy-driven convection in dendritically solidified binary liquids has 
been investigated by Fowler (1985), Nandapurkar ©/ a/. (1989), and Worster (1992) by 
means of linear stability analysis. In all three studies, the mushy zone was modeled as a 
porous medium, in which the kinetics of the phase change are sufficiently fast so that 
nucleation and other nonequilibrium effects can be neglected. Thus, the temperature and 
composition of the mushy zone were taken to be related by the liquidus of the phase diagram. 
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Fowler (1985) considered solid volume fraction of solid in the porous layer to be a 
dynamical variable, but did not admit perturbations to the iiquid/mushy-zone or mushy- 
zone/solid boundaries. He used Darcy's law as an approximate description of the conservation 
of momentum in the porous layer, and considered the case in which the solute diffusivity 
vanishes, so that the thickness of the solutal boundary layer adjacent to the Iiquid/mushy- 
zone interface is zero. Thus, the onset of convection must occur via the mushy layer mode. 
Although Fowler's general formulation allowed for the solid volume fraction to be treated as a 
dynamical variable, his analysis was restricted to a case in which that quantity was 
asymptotically zero. Thus, no interaction was allowed between buoyancy-driven convection 
and solidification. 

Nandapurkar et al. (1989) used the momentum equation developed by Ganesan & Poirier 
(1990) from systematic averaging of the Navler-Stokes equation over a "representative 
volume" large compared to the microscopic scale (primary dendrite arm spacing) and small 
compared to the characteristic length scale (i.e., thickness) of the mushy zone. In their 
analysis of the one-dimensional basic state, the volume fraction of liquid was computed from a 
thermodynamically self-consistent formulation. However, the volume fraction of liquid was 
not allowed to be disturbed in the stability analysis, which is thermodynamically 
inconsistent. Moreover, the stability analysis of these authors did not allow the mushy- 
zone/liquid interface to move (relative to an inertial frame) or deform. These assumptions 
significantly suppressed potential interactions between convection and solidification. 

Like Fowler, Worster (1992) uses Darcy's law in the porous medium and takes the 
porosity to be a dynamical variable, but admits disturbances to the location and shape of the 
Iiquid/mushy-zone boundary, although not to the mushy-zone/solid interface. He concluded 
that onset via the boundary layer mode is relatively independent of the structure of the mushy 
zone, while onset via the mushy layer mode is strongly dependent on the structure of the 
porous medium, and in particular on the porosity profile in the mushy zone. 
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We use a linear stability analysis to study the onset of buoyancy-driven convection in the 
horizontally unbounded liquid and underlying mushy zone of a binary alloy undergoing 
dendritic solidification. The mushy zone is modeled as a porous medium with anisotropic 
permeability. The local porosity, as well as the locations of the boundaries between the solid 
and mushy zone and between the mushy zone and liquid, are taken as dynamical variables, to 
be determined as part of the solution. This work serves to identify the mechanism by which 
Coriolis effects affect the onset of buoyancy-driven convection in dendritically solidifying 
binary alloys. The analysis is illustrated by results for the Pb-Sn system. 

This Chapter is organized as follows. We present the governing equations in §4.2 and an 
appropriate nondimensionalization in §4.3. The one-dimensional basic state and linear 
disturbance equations, both derived from the dimensionless equations, are given in §4.4. The 
numerical solution technique is described in §4.5. Results for Pb-Sn binary alloys are 
presented in §4.6 and discussed in §4.7. 

4.2. Formulation 

4.2. 1. Governing Equations 

A schematic of the solidifying system is shown in figure 4.1. A binary liquid of infinite 
horizontal extent in which the bulk mass fraction of solute is C„ solidifies at a constant 
nominal rate V, due to cooling from below, and is rotated at a constant angular velocity ft 
about a vertical axis parallel to the gravity vector g - -gi z . A mushy zone of nominal height 
H m (to be determined as part of the solution) lies between the melt and solid. For 
convenience, we adopt the practice of previous investigators (e.g., Coriell et al. 1980; Fowler 
1985; Worster 1992) and write the governing equations in a reference frame translating 
with the nominally steady velocity (i 2 Vj) of the moving interface and rotating at a constant 
angular velocity ft about l 2 , but refer the velocities to the laboratory frame. 

The fluid density, p f , is taken as a constant (p 0 , the liquid density at its melting point for 
C„ * 0) except in the buoyancy terms in the liquid and mushy-zone momentum equations, 
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where it is approximated by (4.1). The density of the solid Is also taken as p 0 , so that the 
effects of solidification shrinkage are neglected. 

In the liquid, the motion is governed by the Oberbeck-Boussinesq equations 

V • U[_ = 0 , (4.2a) 
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l 2 + v V 2 u L , (4.2b) 


in a reference frame translating with the nominally steady velocity (l 2 V,) of the moving 
interface and rotating at a constant angular velocity £2 about i 2 , where u L and p L are the 
velocity and pressure in the liquid, and v is the kinematic viscosity. The equations for 
conservation of energy and species in the liquid are 



(4.2c) 
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vc l =d l v 2 c l . 


(4. 2d) 


where T L , C L , k l , and D L are the temperature, solute concentration, thermal diffusivity, and 
the solute diffusivity in the liquid, respectively. Since the reference state corresponds to the 
pure solvent at its melting point, the equation of state in the liquid can be written as 

Pf = Po [l -“ t(T l -T m ) -ocClJ , (4.2e) 

where a T and a c are the thermal and solutal expansion coefficients, respectively, and T M is 
the melting temperature of pure solvent. 

The mushy zone is modelled as a porous medium saturated by interdendritic fluid. Its 
permeability is taken to be orthotropic, with horizontal isotropy. For conservation of mass 
and momentum, we adopt the model of Ganesan & Poirier (1990) developed by systematic 
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averaging of the Navier-Stokes equations over a "representative volume" large compared to 
the microscopic scale (primary dendrite arm spacing) and small compared to the thickness of 
the mushy zone. In a reference frame translating with the nominally steady velocity (1^) of 
the moving interface and rotating at constant angular velocity £2 about i 2 , the equations for 
conservation of mass and momentum In the mushy zone are 

V.u m =0, (4.3a) 

♦ • ^) + 2 “ x Um + 40 x ( “ xr 1 = 

$£i9j _ v <(£(♦) u m + vV 2 u m , (4.3b) 

Po Po 

where u m and p m are the superficial velocity and pressure in the mushy zone, <t> is the volume 
fraction of liquid (henceforth referred to as the volume fraction), and P(4») is an orthotropic 
tensor with horizontal isotropy whose horizontal and vertical components are 
p h ($) - 1/K h (<}>) and P v (<t>) * 1/K V (<|>), respectively, with K h (<|>) and K v (<t>) the components 
of the orthotropic permeability tensor in the mushy zone. The functional dependence of K h and 
K v on $ is discussed in Appendix D. 

The mushy zone is considered to be a continuum, in which the porosity variation is 
continuous. The thermophysical properties of the mushy zone (specific heat, thermal 
conductivity, permeability, and solute diffusivity) are taken to be functions of the temporally 
and spatially varying volume fraction <|>. The models we use for conservation of energy and 
solute in the mushy zone have been employed by several investigators (Worster 1986; 
Nandapurkar et al. 1989). These conservation equations are also formulated on a scale large 
compared to the spacing between primary dendrite arms but small compared to the thickness 
of the mushy zone. The temperature and interdendritic solute concentration are assumed to be 
uniform in the representative volume. Transport of heat and solute in the mushy zone is by 
diffusion and convection. Moreover, as interdendritic fluid solidifies, latent heat is released 
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within the mushy zone and solute is rejected into the interdendritic fluid. Therefore, the 
conservation equations for energy and species must include terms proportional to the rate of 
change of volume fraction. They can be written as 

+c P.L“m- ^ m = V.k nl VT m -L|4 + V,L^. (4.3c) 

* lt“ V I* 1? + ‘ VC ™ ’ V * D ™ VC "> * V,(.-k)C m & . (4.3d) 

where T m is the local temperature in the mushy zone, C m is the solute concentration in the 
interdendritic fluid, c p m and c pL are the specific heats of the mushy zone and liquid, 
respectively, k m and D m are the thermal conductivity and solute diffusivity in the mushy 
zone, k is the segregation coefficient, and L is the latent heat of fusion per unit volume. The 
thermophysical properties of the mushy zone are taken as volume-fraction-weighted 
averages of the fluid and solid phases within the mushy zone 


= <I>c P( l + - 40 c p,s * 

(4.3e) 

* 4> k L + (1 - *) k s • 

(4.3f) 

■ 4>D L , 

(4.30) 


where c p>s and k s are the specific heat and thermal conductivity of the solid. In the mushy 
zone, mass diffusion in the solid phase is orders of magnitude slower than in the 
interdendritic fluid, and so is neglected. The solidification kinetics in the mushy zone are 
taken to be sufficiently fast that nucleation and other nonequilibrium effects can be neglected. 
Hence, the temperature and composition within the mushy zone are related by the liquidus 
relationship, a linear approximation to which is given by 


"^m — T m + n\ C m , 


( 4 . 3 h ) 
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where m L is the slope of the liquidus. Equations 4.3(c,d) are written as two evolution 
equations for the algebraically related temperature and interdendritic solute concentration In 
the mushy zone. Thus 4.3(c,d) should be interpreted as evolution equations for 4> (time 
derivatives of which appear on the right-hand sides), and either T m or C m . The equation of 
state of the interdendritic fluid in the mushy zone is 

Pf = P0 f 1 “ a T( T m- T M)-OC C mj ■ (4.3i) 

In the solid the energy equation is taken as 

^-V.^v’Ts, (4.4) 

where T§ and k s are the temperature and thermal diffusivity of the solid, respectively. 

4.2.2. Boundary conditions 

The boundary condition applied far from the liquid/mushy-zone interface as z — » « is 

C L -» C„ , (4.5) 

where C» is the bulk mass fraction of solute. At the liquid/mushy-zone interface, 
Tj L (x,y,t)-z- 0, the volume fraction, velocity, normal stress, and shear stress are 
continuous 

4,-1 , (4.6a) 

u L -u m , (4.6b) 

£ (Prn - Pl) - vn • j[vu m+ (vu,*) 7 ] .n - [vu L+ (vu,) 7 ] . n} - 0 . (4.6c) 

t. {[pu„,+ (Vu m ) 7 J .n - [vu l +(Vu l ) t ] .nj. 0 , (4.6d) 
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where n and t are normal and tangential unit vectors, respectively. We also require 
temperature and composition to be continuous 


T L -T m , 


Ol_ ~ 0 m i 


(4.6e) 


(4.6f) 


at the liquid/mushy-zone interface, and that the energy and solute balances across that 
interface 


VT l . n - VT m . n , 

VC L . n - VC m . n , 

be satisfied. 

O n TlsfW) - z - 0, the general frame-invariant mass balance 

( p 2“ pi ) (Piui- P2 u 2 ) . n 
across the interface between two regions reduces to 

u m • n - 0, 


(4.6g) 


(4.6h) 


(4.7a) 


(4.7b) 


because we have assumed that the solid and fluid densities are equal. The no-slip condition is 


u m • t - 0. 


(4.7c) 


The composition and temperature at the mushy-zone/solid interface satisfy the 
thermodynamic conditions 

C m - minfCe, C„/k) , (4.7d) 


T m- T s. (4.7e) 

where C E is the solute concentration on the eutectic isotherm. Conservation of energy across 
the mushy-zone/solid interface requires 
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-L«J»V S . n = k m VT m . n-k s VT S . n , (4.7t) 

where V s is the velocity of the mushy-zone/solid interface referred to the laboratory frame. 

4.3. Nondimensionalization 

4.3. 1. Dimensionless Governing Equations 
We scale the length, time, temperature in the liquid and solid, and the velocity, vorticity, 
pressure, and solute mass fraction in the liquid with H L , H L /V,. G L H L , V,. V,/H L . p 0 V,v/H L , and 
C„, respectively, where H L - D L /V, is the characteristic length associated with the solute 
gradient at the liquid/mushy-zone Interface. The length, time, temperature, velocity, 
vorticity, pressure, and solute mass fraction in the mushy zone are scaled with H m , H m /Vj, 
G L H m , V„ Vj/H m , PoV,v/H m , and C«, respectively, where H m is the mushy-zone thickness. The 

dimensionless governing equations in the liquid are 


V • u L = 0 , 


lP“L u. • Vu. + 2£1* x u L + fl* x (£2* x r) 

sclatL dz L L 


-vp L - 


pf 


POT 


Sc 


l z+ V"U L 




Ul 


VT, 


) 





+ m_ • vCl = v 2 c l 


(4.8a) 


(4.8b) 


(4.8c) 

(4.8d) 


where the dimensionless variables u L , pj_ , ^.andCLare, respectively, velocity, pressure, 
temperature, and solute mass fraction in the liquid, Sc - v/D L and Le - D l /k l are the 
Schmidt and Lewis numbers, respectively, y - V,/(gD L ) 1/3 is the dimensionless 
solidification rate, and H L /V,-Q (D L /g 2 ) 1/3 /Y 2 is the dimensionless angular 

velocity. 


The dimensionless equations In the mushy zone are 


Sc 


V m- U m =0, 

(4.9a) 

* A 

P* L J + "m * V «"(^ L ) + x «m + 5 2 4>n‘ 


- - ~rt 'im + V»u m . 

POT SC 

(4.9b) 


(4.9c) 

- V mC m J = V m . ♦V m C m -5(1-k)C m 

(4.9d) 

T _ t M , m*A 
‘m - r + p- . 

(4.9e) 


where V m - (d/dx m , d/dy m , 3/azJ, the dimensionless variables u m , p m , f m , and C m are, 
respectively, the velocity, pressure, temperature, and solute mass fraction in the mushy 
zone, r = §_k L /[L(gD L ) 1/3 ] is the dimensionless liquid-side temperature gradient at the 
liquid/mushy-zone interface, § - Y g 1 ' 3 H m /D L 2 ' 3 is the ratio of the mushy-zone thickness to 
the characteristic length associated with the solute gradient in the liquid, k^ = $ + (i - 
is the thermal conductivity of the mushy zone scaled by the thermal conductivity of the liquid, 
k* - k s /k L is the ratio of the thermal conductivity in the solid to that in the fluid, 
c p,m = 4> + ( 1 -<fr)c p is the specific heat in the mushy zone scaled by the specific heat of the 
fluid, c p = c P(S /c PiL is the ratio of the specific heat in the solid to that in the fluid, p*» H m 2 P, 
m* = m L C.k L /[(LH m (gD L ) 1/3 )], and T m * - T M k L /[(LH m (gD L ) 1/3 )J. 

The dimensionless energy equation In the solid is 
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where f s Is the dimensionless temperature in the solid, and «•- k s /k l Is the thermal 
diffusivity ratio. 

4.3.2. Dimensionless Boundsry Conditions 
The dimensionless boundary condition applied tar from the Hquid/mushy-aone interface 


asz — > °° is 



Cl ->1 . 

(4.11) 

At the liquid/mushy-zone 

interface hL^.Vm.^l-^-O. the dimensionless 

boundary 

conditions are 

0 - 1 . 

(4.12a) 


u L «Um - 

(4.12b) 

(Pm - ^ PL ) = n « {[ 

v m G m + (v m 0 m ) T ] • n - * [ V “L + ( Vul ) T ] * n } ’ 

(4.12c) 

t • 

+ (VmUo,) 7 ] • n - 1 [vu L + (VU L ) T ] • nj - 0 , 

(4.1 2d) 


T L -4T m , 

(4.1 2e) 


Cl* Cm , 

(4.1 2 f ) 


VT|_ • n « V m T m • n , 

(4.1 2g) 


^ VC L • n ■ V m C m • n . 

(4.1 2h) 

The boundary conditions at the mushy-zone/solid interface iis( x m.ym. x m) ^ 0 are 


u m . n - 0, 

(4.13a) 


u m . t - °, 

(4.13b) 
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C m - minfCg/C., 1/k) , (4.13c) 

A Js 

S T m = T S - (4.13d) 

“ r* V s * n= ^m • n - k VTg • n , (4 . 1 3e) 

where V s is the velocity of the mu shy-zone/solid interface scaled by Vj. 

4.4. Basic State and Linear Disturbance Equations 

4.4. 1. Basic State 

The governing equations and boundary conditions (4.8)-(4.13) in the reference frame 
translating upward with the nominal interface velocity V, and rotating at constant angular 
velocity ft about the vertical axis admit a one-dimensional steady basic state solution. The 
temperature and composition in each region, along with the porosity of the mushy zone, 
depend only on the vertical coordinate z. The motionless basic state in the laboratory frame 
corresponds to a nonuniform superficial velocity distribution in the mushy zone in the 
moving frame. The mushy-zone/solid and liquid/mushy-zone interfaces are horizontal, 
stationary, and located at z - 0 and H m , respectively, with the thickness H m being determined 
as part of the solution. The dimensionless equations governing the steady temperature, 
composition, and volume fraction distributions in this motionless basic state are 


*? +Le 



(4.14a) 


in the liquid, 




= 0 . 


d 

*m 



+ L -e$Cp^(<|>) 



r *m * 


(4.14b) 


(4.14c) 
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^°+W1 - k > C m 
CknA *m 

II 

o 

(4.1 4d) 





■ 

f = lM_ + nL*c m , 

* m p p rn 


(4.14e) 





■ 

in the mushy zone, and 






(4.1 4f) 

* 

In the solid. Here T L . C L . T m . C m . *,andT s are the dimensionless basic state distributions ol 

■ 

temperature and concentration in the liquid, temperature. 

concentration, and volume fraction 


in the mushy zone, and temperature in the solid, respectively. 


— 

The dimensionless boundary conditions applicable to the basic state are 



Cl-1 


(4.15a) 

i 

as z -> ~ , 



= 

♦(i) = i . 


(4.15b) 

i 

t L a)-^T m (D . 


(4.15c) 

i 

£ 

ii 

Ol 

3^ 


(4.1 5d) 

i 

, dT L ($) dT m ( 1 ) 
1_ dz - dz* ’ 


(4.1 5e-f) 

B 

,<£ L ($) dC m ( 1 ) 

5 * - d^ ' 


(4.1 5g) 


at the liquid/mushy-zone interface, and 



- 

k T m (0) -Ts(0). 


(4.1 5h) 

— 


0^,(0) = min(CE/C„, 1/k), 


(4.1 5 i) 
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at the mushy-zone/solid interface. We note that (4.15a-J) constitute ten boundary 
conditions for a ninth-order differential equation system. However, the mushy-zone 
thickness H m is unknown, so that the problem is not overdetermined. Dimensionless basic 
state temperature and concentration distributions in the liquid and the temperature 
distribution in the solid can be determined by integrating (4.14a,b,f) subject to (4.15a), 
(4.1 5e), and a linear combination of (4.15f,g) to get 

T L (z) = -1 +±|i -exp[-Le(z-^)] J , (4.16a) 

CJ 2 ) - 1 - rjjn «*[-<*-«] . (4.16b) 

T s (z)=A 1 + A 2 exp^-J=pj, (4.16c) 

where A-) and A 2 are integration constants. The remaining nonlinear ordinary differential 
equations (4.14c, d) subject to the thermodynamic constraint (4.14e) are integrated subject 
to (4.15b,d,g,i) using a shooting technique to obtain T m ,C m , and$ along with the mushy- 
zone height H m . The integration constants A 1 and A 2 in (4.16c) can be found using (4.15h,j). 
Our basic state composition and volume fraction distributions in the mushy zone are in 
excellent agreement with those of Flemjngs (1974). 

4.4.2. Disturbance Equations and Boundary Conditions 
To determine the conditions under which disturbances grow or decay, we write the liquid 
velocity, temperature, composition, and pressure as 


u L (x.y,z,T) = 0 + ujx.y.z.t) , 


(4.17a) 
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T L (x.y,z,x) =T l (z) + TJx.y.z.x) . 

(4.17b) 

C^x.y.z.x) =C L (z) + C' L (x,y f z,x) , 

(4.17c) 

p L (x,y,z,x) =p L (z) + PL(x.y,z,x) , 

(4.1 7d) 


the velocity, temperature, composition, pressure, and volume fraction in the mushy zone as 


Um( x m>ym' Z m' T m) * 0 + u m( x m>ym* z m' t m) * 

A • . 

T m( x m'ym' z nv T m) = T m ( z m) + T m( x m*ym- z m* T m) • 

Cm( x m , ym» z m*' t in) — ^m( z m)"^ ^fn( x m’yrn ,z fn> x ni) * 

_ * 

Pm( x m’ym» z m ,x m) = Pm( z m) + Pm( x m»ym» z m ,T m) » 
<|)(x m ,y m ,z m ,x m ) = ♦(z m ) + ♦( x m ,y m .z m ,t m ) . 

the solid temperature as 

T s (x,y,z,x) =T s (z) +Ts(x,y f z,x) , 

the locations of the liquid/mushy-zone and mushy-zone/solid Interfaces as 

, nL( x fn*ym' T m) = + T lJ x m , ym> T m) • 

Tts( x m-ynv x m) = 0 + • 

and the local velocity of the mushy-zone/solid interface as 

V s (x,y,x) = 0 + Vs(x.y.x) . 


(4.1 7e) 
(4.1 7f) 
(4.17fl) 
(4.1 7h) 
(4.171) 


(4.1 7 j ) 


(4.17k) 

(4.171) 


(4.1 7m) 


where u^, T^, C* L . p^. and u' m , T,; . c; , p* . are the dimensionless velocity, temperature, 
concentration, and pressure disturbances in the liquid and mushy zone, T s is the 
dimensionless temperature disturbance in the solid, 4 is the disturbance volume fraction, 
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Hl and ii s are the dimensionless disturbances to the liquid/mushy-zone and mushy- 
zone/soiid interface positions, and Vg(x,y,t) = Is the dimensionless disturbance 

to the mushy-zone/solid interface velocity. 

Substituting (4.17a-m) into (4.8)-(4.13), subtracting the basic state equations, and 
retaining only linear terms, we obtain dimensionless linear disturbance equations and 
boundary conditions. The equations governing small disturbances in the liquid are 



V . U|_ = 0 , 

(4.18a) 

'IF 
1 1 

- ^ 2Q ' * “l] - - VPL* ^ Tj, * C L 

(4.18b) 



(4.18c) 

3C l dCi dCi ' 2a ; 

-*--* + -* 


(4.1 8d) 


i 

where w L is the disturbance to the z-component of the liquid velocity, F T - a T LD L /p 0 k L , 
F c = acC«/p 0 , £2* = n 0 * I z is the angular velocity, and Qq = ft 0 (D L /g i 2 ) 1/3 /y 2 is the 
dimensionless rotation rate. Taking the curl of (4.18b) and the dot product of the result with 
i z , we obtain an equation 


Sc 




- 2a 



= V 2 <Dl 


(4.19) 


for the z-component of the disturbance vorticity in the liquid, © L ! After twice taking the curl 
of (4.18b), the vertical component of the resulting equation is 





V 2 T' + 



v xC' L + V 4 w'. 


(4.20) 
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where vj - 3 2 /3x 2 + 3 2 /3y 2 is the horizontal portion of the Laplacian operator. The 
equations governing small disturbances in the mushy zone are 


V 


m * 



(4.21a) 



(4.21b) 



aT m._ c * m (<j>) ^ GL - ( 1-c*) ^♦ + £lw 
3x m prn ^' 3z m P dZm 



V m - k^(4>)V m T m 


+ (1 - k ’> i 



(4.21c) 






az, 


m 



(4.21 d) 





(4.21 e) 


where w^, is the disturbance to the z-component of the fluid velocity in the mushy zone. 
Taking the curl of (4.21b) and the dot product of the result with l 2 , we obtain an equation 


■i.rd a ?xn--$-iL 

Sc L 



“ ^m m m ~ 


(4.22) 


for the z-component of the vorticity in the mushy zone, a> m . After twice taking the curl of 
(4.21b), the vertical component of the resulting equation is 
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li]- ^ 

+ ^ Cm + p,(?>vi " w "> + 4; [ Ph(i) * v " - 


(4.23) 


The disturbance energy equation in the solid is 

LefZL-ri 




!Is__£ls.V 

dx dz J 


=k*v 2 t; ■ 


(4.24) 


The far-field disturbance boundary conditions are 


W L 


dw L 

dz 


9 9 


= © L - c L « t l -> 0 , 


(4.25a-e) 


as z -» bo, and 


dT, 


dz 


(4.26) 


as z -> The disturbance boundary conditions linearized about the nominal liquid/mushy- 
zone interface position (z m - 1) are 


4>'(x m .ym.i.*m) + ^^<(x m .y m .T m ) = 0, 


(4.27a) 


WL(x,y^,x)= w^(x m ,y m ,1,x m ) , 


(4.27b) 


,awL (x,y.4.T) _ dw m (x m ,y m ,l,T m ) 
s dz 


dz, 


m 


(4.27c) 


► Z a^X.y.^.T) d^X^.I.Q 


(4.27d) 


£ 3jJllx J ^ = a_w m (x m d&jQ aw m (x m ,y m .1,x m ) 

to 3 dz* & dz, dz n ■ (4 -^ 7e) 


L m 


HU 
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£ cojx.y^.x) ■ cD m (x mf y mi 1 ,x m ) , 

*2 9cQL( x »y tX) _ ^m(^m»ym^ » T ni? ^ 
^ dz az m 1 

T^X.y^x) * £ 9 

Ci(x t y f \ 9 x) * C m (x m ,y m ,1 ,x m ) , 

ai L (x y ^T) ^djM) ^ , tm) . 

9z dz 


9 z 


m 




dz &r 


(4.27f) 

(4.27fl) 

(4.27h) 
( 4 . 2 7 i ) 


, 2 — 

aT m (x m .y m J.x m ) + O m pi^ (XmtymtTm)t (4.27j) 


ac m< x nyym> 1 >?.flil + ^niil) 11 '(x m ,y m ,x m ), (4.27k) 

8z m dzi 


and those at the mushy-zone/solid interface are 

w m (x„,,y n ,.0,x m ) - 8w m< . Mm (x m . yn ,.0.i m ) - 0. <4.28a-c) 


<JT m ( 0 ) . 


Ts(x.y,0 1 x)+ 4 ^-^ns( x m*ym* x m) “ S T m( x m>ym«°- x m) + S — ^ x m«ym> x m)'( 4 - 28d ) 

CJx.y.O.x) + = °» (4.28e) 

♦'(x m .y m .o.t m )-f ♦(O) a V«m*"* s> = C(») 3Tm(X m^’°^ l 

_ k * 3T 9(x-y.P ^ll + k^*) ^^T,^(x m ,y m ,x m ) - k* d -^^Tis(x mi y m ,x m ) 

3z dz^ az 

+ (1-k*) ♦ ( x m'ynv°' x m) + ( 1_k ) ~ ^^^ x m»Ym» x m)- ( 4 - 28f ) 


ini i: 
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4.4.3. Fourier Decomposition 

We write the horizontal and temporal dependence of the dimensionless disturbance 
quantities in the liquid and solid as exp(cn + ia x x + ia y y), where a x and a y are the x- and 
y-components of the horizontal wavevector, and a Is the temporal eigenvalue. The 
dimensionless disturbance quantities in the mushy zone are similarly decomposed as z- 
dependent amplitudes multiplied by exp(a m t + ia^x + ia m y y). Substitution into (4.18c-d), 
(4.19), (4.20), (4.21 c-e), and (4.22)-(4.24) yields 


i fi i d i 0 n 0 do, rFy Fp 2 

H * L Xl ' 01 (4.29a) 


r, . 1 d 1 V 0 flo dW|_ „ 

^ Scdz Sc J L Sc dz 


(4.29b) 


0 


L 2 + Le ^ - ct Le'j © L - Le 5* W L = 0 , 




cb 


(4.29c) 


in the liquid, 


^*i-c)x L -^W L .0. 


{ L 2m| L 2m*|( ^ -o m )L 2m l-P v «) a 2 + 




dz, 


m 


[ p -<]} 


w, 


m 


Sc dz,,, y 4 Sc 


y’Sc 


(4.29d) 


_ii2i 3Sm. §!lfl &L a 2 x -o u 

Or, . 4 P . a m w m- -rr^ a m X m = 0 , (4.29e) 


[ l - + £♦ -k°-h + « • 


( 4 . 2 9 f ) 
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k m( < f>) L 2m + 


dk m (<t>) d _ + 1 iS-+ Le$c p * m ( 4 >) — Le5c pm (<{>) <r m 0 m 

dz m dz m dz,,, dz m ^ J 




^ ^ + ^ + 6 ( 1 ' k 1 ' 5 H I,r ' ^ Wm 

+ [^ +5(1 . k) C m i^-« , -K , C m o m ]^0. 


In the mushy zone, and 


ft — —X 

w m “ p A m » 


^k* L 2 + Le Le oj0s = 0 


(4.29h) 


(4.291) 


(4.29j) 


in the solid, where W L , Q L , X L , and 0 L are, respectively, amplitudes of disturbances to the 
vertical velocity, vorticity, solute distribution, and temperature in the liquid, W m , Q m , X m , 
0 m , and are, respectively, amplitudes of disturbances to the vertical velocity, vorticity, 
solute distribution, temperature, and volume fraction in the mushy zone, 0 S Is the amplitude 
of the disturbance temperature in the solid, and L 2 = d 2 /dz z - a 2 and L 2m «d 2 /dz 2 -a^ are 
linear operators. Here a and a m are the dimensionless horizontal wavenumbers in the liquid 
and solid, and mushy zone, respectively. Matching the solutions at the two interfaces requires 
that the dimensional horizontal wavenumbers be the same. This gives the ratio of the 
dimensionless wavenumbers as a m /a - % « H m /H L . A similar matching procedure for the 
temporal dependence gives the relationship o/a m - S - H m /H L . 

The dimensionless disturbance quantities satisfy the boundary conditions 

* ( 1 ) + &m p =o, (4.30a) 

dZfn 

W l (£) = W m (1) . 


(4.30b) 
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dW L (^) dW m (1) 
s dz dz m ' 

(4.30c) 

2 d 2 W L (^) d 2 W m (1) 

(4.30d) 

3 d 3 w L (5) d 3 w m (l) fsiiidw.ll] 

dz 3 dz 3 Sc dz,,, dz,,,, 

(4.30e) 

§«l(« -n m (l) . 

(4 .3 Of ) 

2 dQ L (S)_ dQ m (1) 

Ct dZm ’ 

(4.30g) 

®l(« - 4 e m (i) , 

(4.30h) 

* L ($) - ZmO) . 

( 4 . 3 0 i ) 

de L«>, 6 ^)a - de m (1), «*m<1>. 
* ; dz 3 dZm dz 3 Pl ' 

( 4 . 3 0 j ) 

* ^ dz 2 Pl dz, ' Cfe2 p «-' 

(4.30k) 


at the liquid/mushy-zone interface, and 


W m (0) = ^^21.Q m (0).0, 


®s(°)+ e m(0) + 


(4.31 a-c) 
(4.31 d) 


x m (0) + %f^p s = o. 


(4.31 e) 
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[ k " " dz 2 ' *m *m r P 

♦■*« M. k -^2) + [,i-O^x]*(0, =o, 


(4.31 f) 


where the constants p s and p L are the amplitudes of the disturbed mushy-zone/solid and 
liquid/mushy-zone interface positions, respectively. For the boundary conditions far from 
the interface, we set all disturbances to zero 
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dz 


asz-*«> , 


as z — » - . 


(4.32a-f) 

(4.32g) 


4.5. Numerical Solution 

Our objective is to find conditions for which infinitesimally small disturbances decay 
(Re(a) < 0) for all but a finite number of critical wavenumbers (typically one), and are 
neutral (Re(a) - 0) for the critical wavenumber(s). These conditions separate basic states 

that are linearly stable from those that are not. 

In the previous section, we formulated an eigenvalue problem for systems of ordinary 
differential equations on two semi-infinite intervals (in the liquid and solid) and a finite 
interval (in the mushy zone), coupled by boundary conditions at the deformable 
liquid/mushy-zone and mushy-zone/solid interfaces. For convenience, we truncate the 
domain and solve the problem on a finite computational domain. With regard to the far-field 
boundary conditions at z = ± °°, we set all disturbances to zero at z = % + h 


w= dw = x L = e L = n L = o 


in the liquid and at z = -h 


e s = o 


at z = \ + h , 


at z = -h . 


(4.33a-e) 


(4.33f) 
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in the solid. In this work, we have taken h to be 10. (We have checked a number of our 
results using larger values of h, and have found that eigenvalues do not change sensibly for h 
larger than 10.) Since we use Chebyshev polynomials in our numerical solution, we scale the 
liquid, mushy-zone, and solid coordinates by z 1 « ( 2z- 2£- h)/h, z 2 «2z m -h, and 
z 3 = (2z-h)/h respectively, so that each region lies between -1 and +1. The resulting 
system is solved using a Galerkin technique developed by Zebib (1987). 

We approximate the highest derivatives of the amplitudes of the disturbance velocity, 
vorticity, solute distribution, and temperature in the liquid, the velocity, vorticity, solute 
distribution, temperature, and porosity in the mushy zone, and the temperature in the solid, 
by truncated sums of Chebyshev polynomials of the form 


W[ 4 W = X A i T )< z i)- 

j=0 

(4.34a) 

N 

«P(2l) = X B j T j< Z l)- 
j=o 

(4.34b) 

N 

*P( z l) = X C J T j< 2 1>’ 

j=0 

(4.34c) 

N 

et 2 >( 2l ) = XDjTjtz,), 

j=0 

(4.34d) 

W m 4, (2 2 ) - XEjTj(z 2 ). 
j=o 

(4.34e) 

N 

fl ri 2) (2 2 ) = X F j T j( z 2)* 
j=o 

(4 . 3 4 f ) 

N 

Xji 2) (z 2 ) = X GjTj<2 2 ) , 

(4.34g) 


j=0 
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6m 2) ( z 2) = X HjTj(z 2 ) , 


j=0 


(4.34h) 


<I> (1) (z 2 ) = X K j T j (z 2 ) , (4 *34 i) 

j=0 

eg'W = I MjTjt*.) . (4 ’ 34j) 

j=0 

where Tj is the j-th Chebyshev polynomial and the coefficients Aj, Bj, Cj, Dj, Ej, Fj, Gj, Hj, K jt 
and Mj are to be found. Representations of lower order derivatives can be found by integrating 
(4.34a-j) and using standard properties of Chebyshev polynomials. The procedure described 
by Zebib (1987) reduces the problem to a matrix eigenvalue problem 

A$+oB$ = 0, (4.35) 

where o is the temporal eigenvalue, and the elements of the square matrices A and B depend on 
the square of the horizontal wavenumber a 2 = a£ + a* the bulk mole fraction C„, and the other 

dimensionless parameters. 

In §4.6 we characterize the stability of the nominal one-dimensional basic state solution 
in terms of the bulk concentration C„ and dimensionless solidification rate y, with all other 
parameters taken as fixed. For each value of C_, we seek one or more critical values of y 
(denoted by y) such that for y> y, disturbances decay for all wavenumbers, while for y<y‘, 
disturbances grow for all wavenumbers in some range. In order to determine the critical 
value(s) of y, we first determine the neutral curve (y versus the wavenumber a) separating 
those combinations of y and a for which all temporal eigenvalues o lie in the left half-plane (a 
stable basic state) from those for which at least one eigenvalue lies in the right half-plane 
(an unstable basic state). To determine the bulk concentration on the neutral curve for an 
arbitrarily chosen wavenumber a, we first choose a value of y and compute all eigenvalues c 
using (4.35). If all eigenvalues have negative real parts, y is decreased by half the previous 
value; otherwise the new value of y is chosen as twice the previous value. This process is 
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continued until we determine two values of y between which at least one sign change is 
obtained in the real part of the least stable temporal eigenvalue. The dimensionless 
solidification rate on the neutral curve is then determined using a bisection method. 

To compute the critical dimensionless solidification rate y* at which instability first 
occurs (i.e., the maximum on the neutral curve), we arbitrarily choose a wavenumber and 
compute the corresponding y on the neutral curve using the procedure described above. We 
then fix y at the value computed at the previous step and compute the eigenvalues a for a 
discrete set of wavenumbers in a chosen range. From this set we select the wavenumber 
corresponding to the a with largest real part. If this wavenumber is at an endpoint of the 
chosen range, we extend the range to include the wavenumber corresponding to the largest 
Re(o). We then select this wavenumber and determine the corresponding y on the neutral 
curve. We continue this process until the relative change in y is less than 10" 6 . Spectral and 
domain convergence results are shown in Table 3. 

4.6. Results 

The solutions of (4.29)-(4.33) depend on seventeen dimensionless parameters, as 
defined in §§4.2 and 4.3. (The approximate solutions on the finite interval also depend on h.) 
It is therefore not feasible to numerically explore the effects of more than a few combinations 
of these parameters on the onset of convection. In this work, we have thus restricted 
ourselves to the Pb-Sn system which, due to the low melting points of both components, has 
been the subject of several experimental studies. The thermophysical properties of the liquid 
and solid phases are as used by Coriell et al. (1980) at reference conditions corresponding to 
pure lead at its melting point. (The diffush/ity of Sn in Pb used corresponds to an infinitely 
dilute solution.) 

Figure 4.2(a,b) shows the dimensionless basic state solute mass fraction profiles in the 
liquid and mushy zone for C„ - 0.2, y « 2.59 x 10“ 2 , and r* 7.12 x 10~ 2 . (For Pb-Sn 
alloys, these values of the dimensionless solidification rate and liquid-side temperature 
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gradient correspond to V, - 80 psecH and G L = 50 K cnH.) Note that the solute profile 
varies almost linearly with the vertical coordinate z m in the mushy zone, within which 
almost all of the solute variation is confined. In the liquid region adjacent to the 
liquid/mushy-zone interface there is a compositional boundary layer, the thickness of which 
is very small compared to the mushy-zone thickness. (Note that we have used different 
length scales in the liquid and mushy zone. For the mushy zone the length scale is the 
thickness H m , the dimensional value of which is 1.685 cm for this combination of G L , Vj , and 

—3 

C„, and for the liquid the length scale is the diffusion length, which in this case is 3.75 x 1 0 
cm.) For the same values of y, r, and C_, figure 4.2(c) shows the basic state porosity 
profile. Figure 4.2(d-f) shows <t>(z m ) for three other values of C„. The volume fraction is 
largest (<|> - 1) at the liquid/mushy-zone interface, and decreases monotonically to a nonzero 
value at the mushy-zone/solid interface. These results show that as C„ increases, the 
porosity distribution becomes more linear and its value at the mushy-zone/solid interface 
Increases monotonically. Figure 4.2(f) shows that the mushy-zone thickness (H m - 0.527 
cm) for C« = 0.5 is much smaller (so that the liquid in the porous medium is more mobile) 

than for C„ * 0.1, 0.2 or 0.35. 

The dependence of the dimensional mushy-zone thickness, H m , on the bulk mass fraction 
and temperature gradient is shown in figure 4.3(a-c) for three solidification rates. The 
solid, dashed, and dotted curves in each figure show H m for y- 3.24x10 
(V,«= lOOpsec' 1 ). 1.62 x 10 -3 (V, « 50 p sec" 1 ), 9.72 x 10" 3 (Vj -SOpsec 1 ), 
respectively. We see that H m is a strong function of r (note that three different logarithmic 
scales have been used for H m ) and C„. depends relatively weakly on y, and assumes a 
maximum near C„ - 0.2 for each combination of solidification rate and temperature gradient. 
This maximum occurs because the composition difference AC m « C m (1)-C m (0) between the 
top and bottom of the mushy zone reaches a maximum near C„ - 0.2, and the mushy-zone 
composition varies almost linearly with z, as shown in figure 4.2. As AC m decreases to zero 
as C„ approaches C E , H m asymptotically decreases to zero. The decrease of H m as C« decreases 
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to the left of the maximum H m is again due to the fact that AC m decreases with decreasing C„. 
For r= 7.12 x 10“ 3 (G L « 5 K cnr 1 ), figure 4.3(a) shows that H m is a decreasing function of 
y. For r» 7.12 x 10 2 (Gl« 50 K cm -1 ), the variation of H m with y is much weaker (figure 
4.3b), and for r» 0.285 (G L « 200 K cm -1 ) is nearly independent of y (figure 4.3c). We 
also note that the mushy-zone thickness decreases strongly with increasing temperature 
gradient. 

For three values of C„, figures 4.4(a-c) show neutral curves (y versus a) for 
r= 7.1 2 x 10“ 2 (G l - 50 K cm- 1 ) and rotation rates ft 0 - 0 rpm (solid curve), 300 rpm 
(dashed-dot curve), and 500 rpm (dashed curve). These neutral curves correspond to steady 
onset of buoyancy-driven convection; we have found no evidence for the oscillatory onset of 
instability in the dendritic solidification of Pb-Sn alloys. For given values of C„ and r, the 
critical solidification rate (denoted by y*) is the maximum point on each neutral curve and 
determines one point on the stability boundary. For y > y* p the one-dimensional dendritic 
basic state solution is stable for all wavenumbers a, whereas for y < y* it is unstable for some 
range of a. 

Figure 4.4(a) shows neutral curves (y versus a) for C. - 0.1. At the extremum, the 
critical solidification rate is y « 2.20 x 10 -2 (Vj * 68 psec -1 ) and the critical wavenumber 
is a - 0.32 for the nonrotating case (solid curve). For Q 0 * 300 rpm the neutral curve 
(dashed-dot curve) is displaced downward for all wavenumbers, indicating the stabilization 
associated with rotation. Although the reduction In y* (toy» 1.82 x 10~ 2 , corresponding to 
V, - 56 p sec- 1 ) is significant (note the logarithmic y-scale), the effect is particularly 
strong at smaller wavenumbers (higher wavelengths). We also note that the critical 
wavenumber has been shifted to a higher value. These effects will be discussed in §4.7. For 
Q o = 500rpm, the critical solidification rate is reduced to y* ■ 1 .571 x 10 -2 
(Vj « 48.5 p sec -1 ) and the critical wavenumber has nearly doubled to a - 0.65. For 

0.2, figure 4.4(b) shows neutral curves (y versus a) analogous to those for - 0.1. 
The critical solidification rate assumes a maximum for this value of C m (nearly independent 
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0 f q q ) [y* . 2.59 x 10 2 (Vj - 80.0 p sec -1 ), 2.30 x 10 -2 (Vj « 71 p sec -1 ), and 
1.86x 10" 2 (Vj - 57.5 n sec -1 ) for Cl 0 - 0, 300, and 500 rpm, respectively) and the 
critical wavenumber increases (a - 0.27, 0.41, and 0.55 for Q 0 “ 300, anc * 5^0 r P m> 

respectively]. The Coriolis acceleration again suppresses the onset of buoyancy-driven 
convection with large wavelength as illustrated in figure 4.4(b). From the neutral curves 
for - 0.55 In figure 4.4(c) we see that aside from a slight shift downward to smaller 
solidification rates [y - 1.98 x 10" 2 <V t - 61.3 p sec -1 ), 1.75 x 10~ 2 (V, - 54 p sec -1 ), and 
1.42 x 10' 2 (V! = 43.9 psec -1 ) forQ 0 - 0, 300, and 500 rpm, respectively], the results 
are qualitatively similar to those for C„ - 0.1 and 0.2. We note that the smallest stable range 
of solidification rates is obtained for C„ * 0.2, which will be discussed below in connection 
with the stability boundaries. 

The vertical structure of the disturbance amplitudes of vertical velocity, temperature, 
and vorticity in the liquid and mushy zone, solute mass fraction in the liquid, and volume 
fraction in the mushy zone are shown in figures 4.5 and 4.6 for fl 0 - 0 and 500 rpm, 
respectively. The disturbances shown correspond to the least stable mode (i.e., that with 
least positive or most negative Fte(a)). For graphical clarity, the same length scale (based on 
mushy-zone thickness) is used in both layers. In these two figures, the mushy zone extends 
from z - -1 to 0, and the liquid lies above. Since temperature and composition in the mushy 
zone are linearly related, only the structure of the temperature disturbance is shown. 

For C„ = 0.55 and r = 7.12 x 10 -2 , figure 4.5 shows the disturbance amplitudes for the 
nonrotating case at the critical conditions y - 1.976 x 10 -2 (Vj - 61 p sec -1 ) and a - 0.32 
(as determined from the neutral curve in figure 4.4c). Figure 4.5(a) shows that the vertical 
velocity disturbance assumes a maximum just above the liquid/mushy-zone interface and 
decays strongly in both layers. We also note that the depth of penetration of the convective 
disturbance In each layer is of the same order. Figure 4.5(b) shows that the structure of the 
temperature disturbance is similar, except near the mushy-zone/solid interface where the 
temperature disturbance does not approach zero. This Is due to the fact that the boundary 


condition (4.31 e) (derived from the dimensional condition (4.7d)) requires that the solute 
mass fraction in the interdendritic fluid be equal to the eutectic mass fraction on the deformed 
mushy-zone/solid interface. Thus, the relationship (4.29j) between the temperature and 
solute mass fraction disturbances in the mushy zone requires that the disturbance 
temperature be nonzero whenever that interface is deformed. Figure 4.5(c,d) shows that 
disturbances to the liquid composition and mushy-zone volume fraction are largely confined 
to very thin layers adjacent to the liquid/mushy-zone interface 

For the same bulk mass fraction and temperature gradient as in figure 4.5, the vertical 
structure of the same disturbance amplitudes, along with the vertical vorticity in the liquid 
and mushy zone, are shown in figure 4.6 for the fl 0 - 500 rpm for the corresponding critical 
conditions [y « 1.458 x 10 2 (V, * 45 psec -1 ) and a « 0.66]. Aside from the considerable 
increase in the critical wavenumber, only quantitative differences are apparent in the 
velocity, temperature, composition, and volume fraction disturbances compared to the 
nonrotating case. The vertical component of the vorticity in the liquid and mushy zone is 
shown in figure 4.6(e) to be highly localized near the liquid/mushy-zone interface and to 
decay strongly within each layer. 

For ■ 0 and 500 rpm, we present our principal results in terms of a division of the 
T*C» plane in figures 4.7-4.10 for four values of r. These figures show the region (bounded 
by the dashed curve) in which the one-dimensional dendritic solution exists along with 
information regarding its stability and the stability of the one-dimensional plane-front 
solution (which exists for all combinations of C., y, r, and fl 0 ). Each region is characterized 
according to whether the one-dimensional dendritic solution exists (and if so, whether it is 
stable), and whether the plane-front solution is stable with respect to morphological and 
convective disturbances (Chapter 2; cf. Coriell et al. 1980). 

Figure 4.7(a) shows the division of the y-C„ plane for r - 0.285 (G L - 200 K cm- 1 ) and 
*= 0 r P m - The one-dimensional dendritic solution exists in regions 1 and 2, and is stable 
at sufficiently high solidification rates (region 1) and unstable for lower values of y (region 
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2). (The turning point on the dendritic stability boundary (curve D) near C„ - 0.2 is 
discussed in §4.7.) In regions 3-5, no one-dimensional steady solution of (4.8)-(4.13) 
exists. The plane-front solution exists in regions 1-5, but is unstable with respect to the 
morphological instability (curve M) in regions 1-3. The plane-front convective stability 
boundary (curve C) separates region 4, in which the plane-front solution is morphologically 
stable but convectively unstable, from region 5, in which the plane-front solution is stable 
with respect to both morphological and convective disturbances. In summary, the one- 
dimensional dendritic solution is stable in region 1 and the one-dimensional plane-front 
solution is stable in region 5; elsewhere the plane-front solution is unstable and the dendritic 
solution either does not exist or is unstable. 

For the same value of r (0.285), figure 4.7(b) shows a similar division of the y-C. 
plane for Q 0 - 500 rpm. The dendritic stability boundary is shifted to slightly smaller 
values of y relative to the nonrotating case, constituting a small enlargement of region 1 , in 
which the dendritic solution is stable. As for disturbances to the plane-front solution, the 
morphological branch is not sensibly changed, but the plane-front convective stability 
boundary has shifted upward (Chapter 2), and so the range of bulk mass fraction and 
solidification rate for which plane-front solidification is stable is significantly enlarged 
relative to the nonrotating case. The greater degree of suppression of buoyancy-driven 
convection that is achieved for the plane-front solution relative to the dendritic case is 

discussed in §4.7. 

Figures 4.8-4.10 show that the topology of the division of the y-C„ plane becomes more 
complex as the temperature gradient decreases. For r - 7.12 x 10" 2 (G L - 50 K cm- 1 ), 
figure 4.8(a-b) shows that for - 0 and 500 rpm. a new region (6) exists in which both 
the plane-front and dendritic solutions are stable. For Q 0 ■ 0 H 3 ™* fi 9 ure 4 -®( a ) shows that 
for the plane-front solution, the morphological and convective stability boundaries have 
shifted downward, and stable plane-front solidification is confined to a smaller part of the 
y-C„ plane than for r - 0.285. However, for this smaller temperature gradient, dendritic 


solidification is stable for a larger range of y and C_ relative lo the r - 0.285 case. For 

7 *'t 2 x 10 , the dendritic stability boundary intersects the dendritic existence curve 
just above the latter's intersection with the plane-front morphological stability boundary. 
Note also that the turning point on the dendritic stability boundary is more pronounced. For 
fto * 500 rpm, figure 4.8(b) shows that the regions in which only one one-dimensional 
solution is stable (regions 1 and 5 for dendritic and plane-front solidification, respectively) 
are considerably enlarged relative to the nonrotating case, whereas the size of the region in 
which both one-dimensional solutions are stable (region 6) is relatively unchanged, since its 
boundaries (the existence curve for dendritic solidification and the plane-front 
morphological stability boundary) are unaltered by rotation. The effect of rotation on the 
dendritic stability boundary is considerably stronger than for r - 0.285. 

For r m 3.56 x 10 2 (G L - 25 K cm -1 ) and « 0 rpm, figure 4.9(a) shows that the 
dendritic stability boundary intersects the dendritic existence curve well below the plane- 
front morphological stability boundary, so that the part of the y-C.. plane in which the plane- 
front solution is stable and the dendritic solution exists (region 6 in figure 4.8) has been 
divided according to the stability of the dendritic solution. In region 6, the dendritic and 
plane-front solutions are both stable, as for r - 7.12 x 10‘ 2 , while in region 7, the dendritic 
solution is unstable. We note that the mushy-zone thickness predicted in the lower part of 
region 6 is exceedingly small; one can conceive of more sophisticated solidification models 
(e.g., specifically predicting cellular solidification and imposing a threshold criterion on 
mushy-zone thickness for dendritic solidification) that would not predict the existence of 
dendritic solution in the lower part of region 6. For = 500 rpm, figure 4.9(b) shows 
that region 7 disappears at the expense of region 6, consistent with the stabilizing influence 
of rotation on dendritic solution manifested in the significant displacement of the dendritic 
stability boundary. 

For r ■ 7.12 x 10 3 (G L « 5 K cm -1 ) figure 4.10(a) shows for the nonrotating case that 
the dendritic solution exists in a part of the y-C m plane (region 8) in which the plane-front 


84 


solution is morphologically stable and convectively unstable. Here, as in region 2, both the 
dendritic and plane-front solutions exist, but neither is stable. The characteristics of the 
other regions are unchanged. (We have not computed the boundary separating regions 6 and 7 
for C„ < 1 x 1 0 -5 because, as discussed above, the basic state mushy-zone thickness predicted 
by our model becomes exceedingly small as the bulk mass fraction decreases further.) Note 
that the stable region of dendritic solidification becomes larger with decreasing temperature 
gradient. For - 500 rpm, figure 4.10(b) shows that region 8 has disappeared at the 
expense of regions 6 and 7, owing to the intersection of the dendritic existence curve and the 
plane-front morphological stability boundary at a larger value of y than the latter s 
intersection with the plane-front convective stability boundary. We also note that the 
dendritic stability boundary intersects the dendritic existence curve at a value of C.. in excess 
of 1 x 10 - ®. We also note that the degree of rotational stabilization for both the plane-front 
and dendritic solutions is considerably greater than shown in figures 4.7-4.9 for larger 
dimensionless temperature gradients. 

4.7. Discussion 

Our analysis of the conditions under which the governing equations (4.8-13) admit a 
steady one-dimensional dendritic solution for Pb-Sn alloys, and the conditions under which 
that solution and the one-dimensional plane-front solution are stable, allows for a more 
complete discussion of the solidification regimes for binary alloys than has heretofore been 
possible. 

That a one-dimensional steady dendritic solution exists for only some combinations of the 
temperature gradient, solidification rate, and bulk mass fraction is not surprising since this 
basic state is determined as the solution of a nonlinear ordinary differential equation system, 
for which a solution need not exist. The boundary in the parameter space within which the 
dendritic solution exists is determined by two constraints. One is that the liquid volume 
fraction at the liquid/mushy-zone interface is unity and cannot exceed that value within the 
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mushy zone. Hence, the vertical derivative of the volume fraction at the liquid/mushy-zone 
interface should be nonnegative. This condition is satisfied in the region above the lower 
branch of the dendritic existence curve. The other constraint is that the solute mass fraction 
in the interdendritic fluid cannot exceed the eutectic composition, C E . The interdendritic 
solute mass fraction in the basic state depends monotonically on the vertical coordinate z. 
Therefore, if the solute mass fraction at the liquid/mushy-zone interface is less than C E , this 
constraint can be satisfied everywhere In the mushy zone. Below the upper branch of the 
dendritic existence curve, the latter constraint is satisfied, and the mushy-zone thickness 
asymptotically approaches zero as this branch of the curve is approached from below. Hence, 
a one-dimensional dendritic solution exists only in the region between the two branches of the 
dendritic existence curve. 

When the plane-front morphological stability boundary is crossed, a cellular interface 
develops. As the solidification rate increases past the critical value, the cells becomes 
deformed and deepen. The formation and evolution of the cells have been studied theoretically 
by Ungar & Brown (1984) and Ungar et al. (1984). For still larger solidification rates, 
dendritic solidification occurs, with dendrites typically consisting of long vertical stems with 
secondary and tertiary branches. The development of the cellular interface and transition to 
dendritic solidification were observed experimentally by Trivedi (1984) and Venugopalan & 
Kirkaldy (1984). We thus conclude that the part of region 3 (cf. figures 4.7-4.10) below 
the lower branch of the dendritic existence curve is a region of cellular solidification, in 
which the plane-front solution is morphologically unstable and the dendritic solution does not 
exist. The experimental evidence showing that dendritic solidification occurs when either the 
bulk mass fraction or solidification rate is increased for fixed temperature gradient is 
consistent with the location and character of the dendritic regions 1 and 2 in our results. 

For a given value of the bulk solute mass fraction, there exists a critical solidification 
rate above which the one-dimensional dendritic solution is stable with respect to 
infinitesimal disturbances. For all values of r, and Q 0 * dendritic solidification is least stable 
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with respect to buoyancy -driven convection for solute mess fractions near 0.2 (the value at 
which the dendritic stability boundary has a turning point). This occurs because for this 
value of C„ the mushy-zone thickness H m and the composition difference AC m attain their 
maximum values independent of temperature gradient and solidification rate (cf. figure 4.3 
and accompanying description). If the Rayleigh number R^, we define for the mushy zone is 
proportional to AC m H m (Worster 1992), maximum destabilization for the Pb-Sn system 
occurs for C„ near 0.2. Fowler (1985) defined a Rayleigh number R- (a c -a T ) m L p 0 gK/Vj, 
where K is a monotonically Increasing function of the permeability, and showed that there is a 
critical value R c «i/(C E -CJ such that if R < R c . then dendritic solidification is stable with 
respect to small disturbances. Note that for Pb-Sn, R c assumes a minimum near C„ * 0.2. 
Thus, there is a critical solidification rate above which the condition R < R c is satisfied, and 
dendritic solidification is convectively stable. Our numerical results are thus consistent with 
the approximate analytical results of Fowler (1985) and Worster (1992). 

We also note that the stable range of y and C. grows as r decreases. This is due to the fact 
that liquid volume fraction (and hence permeability) decreases with decreasing temperature 
gradient, so that the mushy zone becomes more immobile and the onset of convection is 
suppressed. That is also consistent with Fowler's result, since the Rayleigh number 
decreases as permeability decreases. 

The amplitudes of the disturbance to the vertical velocity assume a maximum in the 
liquid ahead of the liquid/mushy-zone interface, and decay strongly in each region. 
Buoyancy-driven convection is largely confined to the liquid ahead of the interface and to the 
part of the mushy zone where the fluid is most mobile (near the liquid/mushy-zone 
interface). The interdendritic fluid farther from the liquid/mushy-zone interface is 
essentially stagnant. For r ■ 7.12 x 10 ^(G^* K cm -1 ), the critical wavenumber range 
is approximately 0.25 < a < 0.35 (depending on C„) for the nonrotating case (figures 4.4a- 
c). The corresponding convective wavelength of 0.06 cm is small compared to the mushy- 
zone thickness (H m - 0.311 cm for C„ * 0.55), but large compared to the diffusion length for 
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a typical solidification rate (4.92 x 10" 3 cm for Vj - 61 psec- 1 at the critical value of y). 
We believe that this convection corresponds to the mushy layer mode described by Worster 
(1992). Rotating the system at Q 0 - 500 rpm suppresses the convective Instability at 
relatively small wavenumbers and nearly doubles the critical wavenumber, so that the 
wavelength of the buoyancy-driven convection is reduced by almost half relative to the 
nonrotating case. 

That a one-dimensional dendritic solution exists under conditions for which plane-front 
solidification is linearly stable with respect to morphological disturbances might at first 
seem inconsistent. However, our result implies no such inconsistency, because the plane- 
front stability analysis (Coriell et al. 1980) is restricted to infinitesimal disturbances. In 
fact, in part of the region where the plane-front solution is linearly stable, it is known 
(Wollkind & Segel 1970; Alexander et al. 1986; Ungar & Brown 1984) to be unstable with 
respect to a subcritical ("finite amplitude") instability. Thus, in the region in which the 
plane-front solution is linearly stable, it is not surprising that we also find a one- 
dimensional dendritic solution. 

The remarkable suppression at modest rotation rates of plane-front convective 
instability in Pb-Sn binary alloys was discussed in Chapter 2 in terms of the well known 
Taylor-Proudman mechanism, described by Chandrasekhar (1961). The Coriolis 
acceleration also inhibits the onset of buoyancy-driven convection in the dendritic case, but 
to a lesser extent than for plane-front solidification. This is because the effective Taylor 
number (a dimensionless measure of the Coriolis acceleration) is smaller in the dendritic 
case, due to the larger "effective viscosity" of the interdendritic fluid in the porous medium. 
The far greater suppression of dendritic solidification at smaller wavenumbers than at larger 
wavenumbers, and the shift of the critical wavenumber to a higher value (both shown in 
figure 4.4) are expected on the basis of the character of the Taylor-Proudman mechanism. 


CHAPTER 5 
Conclusion 

For liquid Pb-Sn undergoing plane-front solidification, a remarkable suppression of the 
onset of buoyancy-driven convection is obtained at low solidification rates (more than a 
hundredfold increase in the critical value of C_ at V, = 1 nsetH by rotating the layers at 500 
rpm for the two values of G L considered). The mechanism of this stabilization is discussed in 
terms of the Taylor-Proudman theorem, in Chapter 2. It is not surprising that the 
morphological instability is unaffected by rotation. The morphological instability occurs at 
very short wavelengths, so the motion Is almost perpendicular to the solidification front (i.e., 
aligned with the axis of rotation). Hence, the Coriolis acceleration does not sensibly affect the 
morphological instability. The predicted inhibitory effects of rotation on convection are 
discussed in terms of previous experimental and theoretical studies of the effect of rotation on 
the onset of buoyancy-driven convection in single-component fluids heated from below and in 
binary fluids subject to thermal and solutal stratification. 

For mercury cadmium telluride H gi . x Cd x Te, the liquid density does not depend 
monotonically on temperature for some range of the bulk solute composition. For certain 
combinations of the operating parameters (solidification rate, nominal liquid-side vertical 
temperature gradient, and bulk solute concentration) a density maximum occurs at a 
temperature greater than the liquidus temperature and there exists a sublayer adjacent to the 
solid-liquid interface in which temperature gradient Is destabilizing. This has profound 
effects on the nature of neutral curves, eigenfunctions, and the stability boundaries. 

The morphological instability is completely insensitive to the net density stratification 
in the liquid, and is relatively insensitive to the variations of temperature and composition 
outside a relatively thin layer adjacent to the interface. This is evident from a comparison of 
the neutral curves and stability boundaries for H 9l . x Cd x Te to those for Pb-Sn. On the other 
hand, the onset of buoyancy-driven convection depends strongly on the nature of the 


stabilizing or destabilizing temperature and solute gradients in the liquid, as well as on the 
diffusivities of the stratifying agencies. For Pb-Sn, the onset of buoyancy-driven convection 
is monotonic in time (corresponding to the fingering regime of doubly-diffusive convection), 
while for Hg,_ x Cd x Te the motion is overstable or oscillatory in time (the diffusive regime). 
This is discussed in terms of the standard "parcel" argument for a linearly stratified doubly- 
diffusive fluid (Turner 1973, page 251) in Chapter 3. 

In contrast to the case where the density depends monotonically on temperature and 
composition (e.g., the lead-tin alloys considered by Coriell et al. (1980)), for Hg 1 . x Cd x Te 
there exists a critical value of the bulk mole fraction below which plane-front solidification 
is unstable at all solidification rates. This differs from the Pb-Sn case in that for Pb-Sn, 
plane-front solidification at any C„ is linearly stable for all sufficiently small solidification 
rates. Moreover, when the density varies nonmonotonically, there can exist a critical value 
of the dimensionless solidification rate y (denoted by y*) such that for y> y\ plane-front 
solidification is unstable for all values of the bulk mole fraction C„. In this case, for y < y*. 
there is a finite range of for which the plane-front solution is stable. This latter result 
differs from the Pb-Sn case, for which at all dimensionless solidification rates, plane-front 
solidification is stable for all values of C.. lying below some critical value. 

For Hg, _ x Cd x Te, the disturbances associated with the oscillatory onset of buoyancy-driven 
convection are largely confined to that thermally unstably stratified sublayer. The results 
should be expected to apply to the solidification of other binary or multicomponent liquids in 
which (3p L /3T L )(9T L /3z) changes sign within the layer. 

The effect of rotation on the stability of plane-front solidification of Hg^CdjJe 
pseudobinary liquid has also been investigated. Similar to the case for Pb-Sn alloys, the 
onset of buoyancy-driven convection is inhibited significantly by modest rotation rates, but 
the morphological instability is not sensibly altered. 

In Chapter 4, we show that a one-dimensional steady mushy-zone solution exists for only 
some combinations of the temperature gradient, solidification rate, and bulk mass fraction. 
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This Is not surprising since this basic state is computed from a nonlinear ordinary 
differential equation system for which a solution need to exist. The one-dimensional mushy- 
zone solution is stable for sufficiently large solidification rates for all wavenumbers. At 
lower solidification rates, it is unstable with respect to disturbances which grow 
monotonically in time. We also found that the dendritic basic state solution exists in part of 
the region where plane-front solidification is linearly stable with respect to morphological 
instability. This does not imply any inconsistency between our models of plane-front and 
dendritic solidification, since the plane-front stability analysis is restricted to Infinitesimal 
disturbances. In part of the region where the plane-front solution is linearly stable, 
previous studies have shown that it is unstable with respect to subcritical instabilities. 
Hence, it is not surprising that we find a one-dimensional dendritic basic state solution in the 
region where the plane-front solidification is linearly stable. 

We have also found that the Coriolis acceleration suppresses the onset of buoyancy- 
driven convection during dendritic solidification of Pb-Sn binary liquids. The relative degree 
of stabilization is less than for plane-front solidification. This is due to the fact that the 
effective Taylor number (a dimensionless measure of the Coriolis acceleration) is smaller in 
the dendritic case, due to the larger effective viscosity of the fluid in the porous medium. 

Even though the analysis in this dissertation is restricted to horizontally unbounded fluid 
layers, the predictions will be qualitatively correct for finite aspect ratios (ratio of mold 
radius to height) if the parameter Ro/g (a Froude number, where R 0 is the mold radius) 
is sufficiently small, as discussed in detail in Chapter 2. For the onset of thermal convection 
in rotating water or mercury layers heated from below, the excellent quantitative agreement 
between classical linear stability analysis for a horizontally unbounded layer and 
experimental work for finite aspect ratios provides a clear demonstration of the potential of 
the Coriolis acceleration to suppress buoyancy-driven convection in a rotating fluid. 

The relatively modest rotation rates required to significantly inhibit the onset of 
convection in the Pb-Sn and Hg^CdJe systems make the proposed method an interesting 
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candidate for a program of laboratory experiments. Experiments might be conducted using a 
completely filled cylindrical mold, thermally insulated on the vertical surface, and mounted 
axisymmetrically on a rotating horizontal turntable. If the liquid at the top of the rotating 
mold is in contact with a gas or vacuum, it will have a nearly paraboloidal free surface, on 
which the elevation above the point on the axis of rotation is Az - fl^r 2 /( 2g). Although this 
configuration is consistent with rigid-body rotation, it also leads to a nearly paraboloidal 
solid-melt interface. For high rotation rates or large mold radius, this will in turn lead to 
significant radial variations in the solidified alloy. Also, cooling at the radial boundary leads 
to a radial temperature gradient, which in turn leads to significant centrifugal effects. 
Therefore, the present results can not compared to the experimental works in which strong 
radial variations in macrosegregation are observed in the solid. 

In our models for both plane-front and dendritic solidification, the domain considered is 
horizontally unbounded and, the liquid and solid extend vertically to z - «. and z - -«», 
respectively. The interfaces in the basic state are assumed to be nominally planar. We also 
assume that solidification occurs sufficiently fast at the growing interface or in the mushy- 
zone so that the solid-liquid interface for plane-front solidification, or the entire mushy- 
zone for dendritic solidification can be considered to be in thermodynamic equilibrium. Thus 
nucleation and other kinetic effects are neglected in our models. 

We have also assumed that the solid is thermally isotropic, and that for plane-front 
solidification, the liquid-solid interface is isotropic. The effects of boundary thermal 
anisotropy have been dealt with in Rayleigh-B6nard problem by Pearlstein & Oztekin 
(1989), while anisotropic surface tension has been considered for plane-front solidification 
by Voorhees at at. (1 984). Both of these effects can be incorporated easily into the present 


work. 
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Appendix A 

The numerical solution of (2.17)-(2.19) can be obtained using a spectral Galerkln 
technique developed by Zebib (1987). We approximate the highest derivatives of the 
amplitudes of the disturbance velocity, vorticity, solute distribution, and temperature in the 
liquid, and temperature in the solid, by truncated sums of Chebyshev polynomials of the form 


J 


W (4) (z 2 ) = X K i T J< z z>’ 
i= o 

(Ala) 

j 

© (2) (Z 2 ) = S M j T J (2 2> * 
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(Alb) 
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II 

CM 
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SL. 

(Ale) 

e[ 2) (z 2 ) = 

j=0 

(Aid) 

j 

©i 2) (Zi) = X Z j T j (z l). 

(Ale) 


J=o 


where the coefficients Kj, M j( Pj. Q jp and Zj are to be found. Representations of lower order 
derivatives can be found by integrating (Ala-e) and using the following standard properties 
of Chebyshev polynomials 

? T M - ttn dT n+i( z > htZ dT n-1< z ^ (A2a) 

2Tn(z) n + 1 (fen-1 dz ’ 


2z T n (z) ■ a n T n+ i(z) - X, n .iT n _i(z) , 


(A2b) 
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where a n ■ X„ ■ 0 if n < 0, a 0 - 2, X 0 ■ 1, and a n -Xn-1 if n > 0. The Integration constants 
in the lower-order derivatives can be determined as functions of Kj, M jt Pj, Qj, and Zj using the 
boundary conditions (2.19), leading to 


J J+4-m 

W (m, (2 2 ) = X X [ w ji (m) Kj + y ! 2 ji (m) Mj + y 1 3 ji (m) P, + y u ,| (m) Qj + y 1 5ji (m) Zj] Tj( z 2 ) + 
ho j=o 


J+3-m 


X 9l| (m> P T l< z 2)- 


(A3a) 
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<o (n) (z 2 ) = I 5> j + y 2 1 jj (n) Kj + y 23 j| (n) Pj + Y24 ji <n) Q| + Y25 ji (n>z j] T j( z 2 > + 

W) i=0 


J+1-n 


X 9 2 j <n) P Tj(z 2 ), 
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X94) (n) PT j( z 2 ), 
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(A3d) 


(A3e) 


for m = 0,1 ,2,3 and n * 0,1. Substituting (A3a-e) into the governing equations and taking 
inner products with 


J+4 J+2 J+2 J+2 

Sw/O'r,. 5 > L|I <»'t i , 2 X 1 ,'% 

H> j-o j=o H> 


J+2 

Xe Sji (0) Tj, (A4) 

j-0 


and 
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respectively, we reduce the solution of (2.17)-(2.19) to a matrix eigenvalue problem 

As + aBs = 0, (A5) 

where c is the temporal eigenvalue, A and B are square matrices, and s= [Kj, Mj, Pj, Qj.Zj, P] 
is the eigenvector. The choice <A4) eliminates the spurious eigenvalues frequently 
encountered in the solution of eigenvalue problems for ordinary differential equations when 
using basis functions not satisfying the boundary conditions. 


95 


Appendix B 

To the best of our knowledge, the only data available for the density of liquid Hg l0( Cd x Te 
are those of Mokrovskii & Regel (1952) for pure HgTe, and Chandra & Holland (1983) for 
CdTe mole fractions in the range 0 £ C £ 0.2. On the basis of the latter, three equations of 
state (Chandra 1985; Antar 1987, 1988, 1991; Apanovich & Ljumkis 1991) have been 
proposed for liquid Hg^CdxTe. Unfortunately, for C - 0, Chandra & Holland's data differ 
systematically from those in Mokrovskii & Regel's well-known paper from Ioffe's Physico- 
Technical Institute in St. Petersburg. Chandra & Holland refer to two other papers by Regel, 
but make no reference to his work with Mokrovskii on the density of HgTe. 

Although the relative difference between the C - 0 data of Chandra & Holland and 
Mokrovskii & Regel is only about 0.65%, the average difference amounts to about 60% of the 
density variation reported by Mokrovskii & Regel over 963 KsTs 1173 K and to about 
75% of the variation reported by Chandra & Holland over 953 KsTs 1073 K. Moreover, 
the C - 0.1 data of Chandra & Holland show density extrema at three temperatures, whereas 
all other data vary either unimodally (for C £ 0.05) or monotonically (for C > 0.1) with 
temperature. Finally, Chandra & Holland reported data only for the Hg^CdxTe system, 
whereas Mokrovskii & Regel's (1952) density measurements for other liquids are in good 
agreement with those of Lucas & Urbain (1962) and Lucas (1964) for germanium, and of 
Glazov, Chizhevskaya & Evgen'ev (1969) for germanium, gallium antimonide, and indium 
antimonide, thus allowing an independent assessment of accuracy. 

We have thus chosen to regard the C - 0 data of Mokrovskii & Regel as accurate, and have 
undertaken to recalculate the data of Chandra & Holland on that basis. Chandra & Holland 
measured the liquid height A in a capillary tube, from which their reported liquid volumes 
V C h were inferred using a linear equation of the form V CH - 8jA + 8 2 . Based on discussion 
with Chandra (private communication), we hypothesize that the primary source of error Is 
in either the determination of A or the determination of V from a. Thus, we have recalculated 
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Chandra & Holland's data (excluding their measurements for C - 0.1 at the two artifactual 
extremal temperatures of 1053 K and 1083 K) by using 

V rc - s^ch + s 2 , (B1) 


and regarding s t and s 2 as unknowns, which we determine by minimizing the sum of the 
squared differences of the specific volumes 


E P-X [pMR(Ti) Prc(Ti)] . < B2 > 

at C - 0, where Pmr refers to the experimental HgTe density data of Mokrovskii & Regel, 
P rc - m HgTe /(s 1 V CH + s 2 ) is the recalculated density of Chandra & Holland for C - 0, the mass 
of pure HgTe used by Chandra & Holland for C - 0 is m HgTe - 62.6457 g, and V CH is the HgTe 
volume reported by Chandra & Holland at each temperature T, at which they measured the 
density for C - 0. The sum E p is minimized by s t - 0.9935 and s 2 - 4.553 x 10 10 cm 3 . 
Noting that the volumes reported by Chandra & Holland are on the order of 10 cm3. we see 
that (B1 ) is nearly equivalent to multiplying the data of Chandra & Holland by a constant. 

This allows recalculation of Chandra & Holland’s reported volumes (for all but the two 
data points cited above) according to (B1). Recalculated values of Chandra & Holland's data at 
all mole fractions are used to fit by least-squares a bivariate polynomial of the form (3.7a) 


in the temperature range of their density data. 
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5.885 x ^0~ 4 gcm^K' 1 

b 30 

- 
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The coefficients bjj are 

b 01 «= -1.578 x 10 -4 gcm^K * 1 
b 03 - 4.268 x 10" 8 gcm^K ' 3 
b 10 - -2.419 gem ' 3 

b 12 - 1.145 x 10 -5 gcnfV 2 
b 20 - -2.181 gem ' 3 

b 22 - -2.049 x 10" 4 gcm^K * 2 
b 31 - 4.436 x 10~ 2 gcm'V 1 . 
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Far from the liquid-solid interface, however, the liquid temperature exceeds the 
maximum temperature (T max (C)) at which p L data are available. For T > (C), the liquid 

density is approximated by a linear function of temperature (3.7b). The coefficients d (J in 
(3.7b) are computed by matching the density and its temperature derivative at T » T (C) 
and are found to be 


doo 

- 8.119 


g cm -3 

d-jQ ■ -2.435 

gem -3 

^20 

- -1.985 


g cm" 3 

dso ■ 10.48 

g cm' 3 

d 0 i 

- -7.703 x 

10“ 4 

g cm^K' 1 

d t1 - 2.753 x 10~ 3 

g cm'V 1 

d 2 i 

= -1.358 x 

10~ 2 

g cm'V 

d 31 - 1.302 x 10 -2 

g cm'V 1 


The coefficients r ( in (3.8) are computed by a least-squares fit of the maximum 
temperatures at which Chandra & Holland reported the density for different CdTe mole 
fractions. The values of r ( are 

r 0 = 1073 K r, = 291.7 K r 2 - -2250 K r 3 - 8333 K. 

The results are shown in figure B1. Solid curves represent liquid density determined 
from (3.7a, b) as functions of temperature for C - 0, 0.05, 0.1, and 0.2. Filled circles 
denote the C . 0 density data reported by Mokrovskii & Regel (1952). The other symbols 
denote density data reported by Chandra & Holland (1983). For 0 s C £ 0.16, we obtain a 
unimodal density variation with temperature. (At the value of C at which p becomes a 
monotonic function of. T, the maximum occurs at the liquidus temperature.) Although the 
C= 0 density data of Chandra & Holland are consistently lower than those of Mokrovskii & 
Regel, the density maxima occur at nearly the same temperature. 
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Appendix C 


The thermophysical properties of the liquid and solid (other than the density of the 
liquid) are taken as constants evaluated using the functional forms shown below at the bulk 
mole fraction and corresponding liquidus or solidus temperatures (Tu q and T SO |). 

We have taken the density of solid Hg^CdxTe as 

1_» _L_ + JLJL, (Cl) 

PS PCcfTe PHgTe 

where ; - CM CdTe M is the CdTe mass fraction, the molecular weights of CdTe and HgTe are 
M CdT G “ 240 0 /mol and M HgTe = 328.2 g/mol, M - CM cdTe + (1-C)M HgTe , and p CdTe and 
PHgTe are ^e s °l'd densities of pure CdTe and HgTe, respectively, represented by 

PCdTe = *0 + X 1 T ’ (C2a * 


and 


PHgTe = ♦o + 


(C2b) 


where the coefficients X, and 4 , are determined by least-squares fits to the experimental data 
of Glazov, Chizhevskaya & Evgen’ev (1969) and Mokrovskil & Regel (1952), respectively. 
The values of Xj and <t>j are 


X 0 - 5.820 gem 
<j> 0 - 8.201 gem' 


X! - -8.095 x 10' 5 g cm -3 K -1 
4 > 1 - -1.230 x 10 -4 g cm -3 K" 1 . 


The viscosity is represented over the entire range by 

v = 1 x 10' 2 exp(«i»o + ®i/t) cm 2 /sec. 


(C3) 


where the coefficients «6 0 - -3.401 and ft, - 3445 K are computed by least-squares fit to the 
experimental data of Glazov, Chizhevskaya & Glagoleva (1969) for pure CdTe. Due to the lack 
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of experimental viscosity data for HgTe, we use (C3) for the viscosity of Hg^Cd^Je, 
independent of composition. 

The thermal diffusivity of liquid and solid Hg^CdxTe are approximated by functions of the 
form 



k l - (co 0 + <o 1 C 1/2 + ( 0 2 C) ln(T/K) - 

(q 0 + n 1 c 1/2 + ft 2 C) , 

(C4a) 
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are determined by least-squares fit to the experimental data of Holland & Taylor (1983). 

The specific heat of liquid Hg^CdxTe is represented by 

2 2-1 

Cp.L( T * C ) =Z X q ij C l (T-943 K)i . (C5) 

l«0 j«0 


The coefficients qy are found by least-squares fit to the calculated specific heat of Su (1986). 
The values are 
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where we have used the relation 1 mol - 2 g-atom for the pseudobinary system Hg^CdJe to 
convert Su's data (in cal K -1 g-atom -1 ) to J K -1 mol -1 . 

The specific heat of solid Hg^CdxTe is taken as 

Cp,s * Cp.S.HgTe 0 "°) + Cp.S.CdT# C • * C 6 * 

where the temperature-dependence of the specific heats of pure HgTe and CdTe Is found from 
Mills (1974) to be 

Cp.S.HgTe = I 52 09 + 9 ’ 08 x10 " 3 < T/K >] J (C7a) 

and 

Cp.s.CtfTe = [40.0 + 3.3 X10 -2 (T/K)] J 1C 1 mol -1 . (C7b) 

The latent heat of fusion is approximated by 

L(C) - (1-C)L H g T e + CLcdTe . (C8) 

where data for the pure components are given by Mills (1974) as L HgTe « 3.6x10 J mol - 
and Lcjre - 4.48 x 1 0 4 J mol -1 . The interfacial energy is estimated by comparing Hg^CdJe 
and InSb (Seidensticker & Hamilton (1963) compared Ge and InSb to estimate the interfacial 
energy of InSb) 

°h»mW ‘ ° lnSb (^Isb ’ ( ° 9) 

where the latent heat of fusion per unit volume for Hg^CdxTe is taken to vary linearly with 
solute mole fraction (see (C8)), ao - 6.4797 A for InSb, and a 0 - 6.465 A is used for 
Hgi. x Cd x Te independent of composition (Glazov, Chizhevskaya & Glagoleva 1969). Here the 
surface tension is taken as a |nSb - 8.5 x 10 -6 J cm -2 , and the latent heat of fusion per unit 
volume of InSb is L « 1.20 x 10 3 J cm -3 . The capillary coefficient is then computed from 
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'P = a, 




H 9l-x Cd x T ® L 


(CIO) 


where T Hp is the liquidus temperature for a given bulk mole fraction of CdTe in the liquid. 
The liquidus and solidus temperatures 


T |jq - (943 + 681 C - 372 C 2 ) K (Cl 1 a) 

and 

T«ol - (9^3 + 202 C - 150 C 2 + 324 C 3 ) K , (C 1 1 b) 

liquidus slope 

m L * (681 - 745 C) K , (Cl 2) 

and segregation coefficient 

k - C (0.30 - 2.24 x 10“ 3 K" 1 T sol + 2.67 x 10" 6 K" 2 T s 2 , ) _1 for C > 0.1 , (Cl 3a) 

k ■ 374 for C < 0.1, (Cl3b) 

as functions of CdTe mole fraction were determined by least-squares fits to the experimental 
data of Szofran & Lehoczky (1981). 

Finally, the solute diffusion coefficient is taken as D L - 5.5 x 10' 5 cm 2 /sec, the value 
most commonly used in the literature (see e.g., Kim & Brown 1989). 
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Appendix D 

Our model of flow in the mushy zone is based on the observation that the primary 
dendrites formed during dendritic solidification of binary alloys are columnar and oriented 
parallel to the nominal solidification direction. We treat the dendritic region, saturated by 
interdendritic fluid to be a porous medium the permeability of which is taken to be 
orthotropic with horizontal isotropy. The functional dependence of the components and K v 
of the permeability is determined from the calculated drag coefficients of Sangani & Acrlvos 
(1982a) and Drummond & Tahir (1984) for flows normal and parallel to the axes of 
cylinders in hexagonal and triangular arrays, respectively. The asymptotic solutions for 
concentrated and dilute arrays (solid curves in figure D1) 

h = 1^2 (|n<p- 1 / 2 - 0.745+ q>-0.25<p 2 ) for <p < 0.3 , (D1) 

X 2 4k ' ' 

Kh _ (1 -q>)4V2T 1 _ (_SL_V /2 1 5/Z for <p > 0.7 . (D2 ) 

F ~ 27k [ W J 

are derived by Sangani & Acrivos (1982a), where <p = 1 - * is the solid volume fraction, <p max 
is the maximum value of the solid volume fraction when the cylinders touch, and X is the 
primary dendrite arm spacing. The open circles are permeability values calculated by 
Sangani & Acrivos (1982a) for flow through a hexagonal array of cylinders, and match the 
asymptotes (D1) and (D2) for sufficiently large and small 4>. Hence, for 0.3 < 9 < 0.7, we 
use a functional form 

_ ij=J£(| n <p- 1/2 -b + (p—0.25 <p 2 ) for <p < 0.3 , (D3) 

similar to (D1), and determine b - 0.7508 by least-squares fit to the calculated drag 

coefficients of Sangani & Acrivos (1982a). Our fit is shown for 0.25 < + < 0.725 by the 
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dashed curve in figure D1. The dotted curve of figure D1 shows the permeability calculated 
by Sangani & Acrivos (1982b) for simple cubic arrays of spheres. 

For the normal component of permeability, we use the results of Drummond & Tahir 
(1984) for flows parallel to the axes of cylinders in triangular arrays 


X? * ^ 1-4975+ 2 »-0.5<p 2 -0.002514 9 6 ) 




{■•“fear-’H 


for <p < 0.7 , ( D 4 ) 


for 9 > 0.7 . ( D 5 ) 


The primary dendrite arm spacing X for Pb-Sn systems was determined by Mason et al. 
(1982) using experimental data over a wide range of temperature gradient G L and 
solidification rate V t by the relationship 

X-AGl -0 , (D 6) 

where the coefficients A and B given by Mason et al. (in their table 2) are functions of 
solidification rate. In our calculations, we have determined A and B by interpolation. 
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Figure 2.1. Schematic depiction of plane-front solidification 
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Figure 2.4b. C„-a neutral curves for plane-front solidification of Pb-Sn with G L - 200 K cm* 1 and V, - 30 p sec* 1 for £2 0 
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Figure 35b. Dimensionless liquid-phase basic state density profiles for Hg^Cd-Te for several values of C„ with r - 8.2 x 10 
(G l - 25 K/cm) and y - 2.64 x 10" 4 (V, = 1 p/sec). 





Figure 3.3. 3 (p l /Plo)^ c - versus z for several values of C„ with y - 2.64 x 10 -4 and r - 8.2 x 10 (V, - 1 u/sec and 

G,- 25 K/cm). 
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Figure 3.5d. For Hg^CdxTe with C„ = 0.07, y ■ 2.64 x 10 , and r - 1.64 x 10 -3 (V, - 1 p/sec and G L « 50 K/cm), dimensionless 
basic state and disturbance profiles, i) Basic state density profile; ii) basic state density gradient (dp L /dz)/p L0 (solid curve) and the 
additive contributions from thermal (dotted) and solutal (dashed) stratification, (3p L /dT)(dT L /dz)/p L0 and (dp L /dC)(9C L /3z)/p LO , 
respectively; iii) amplitude of disturbance temperature; iv) amplitude of disturbance solute distribution; v) amplitude of disturbance 
vertical velocity. The disturbances shown are those which are the most unstable for a - 2.85. 
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Figure 3.8. Stable region in y-C„ plane for Hg^CdJe for Oq - 100 rpm and r - 8.2 x 10 -4 (G L - 25 K cm 1 ). The solid curve Is the 
morphological (Cl) stability boundary, and the dashed curves are the buoyancy-driven convective (Cl*) stability boundaries. 


























Figure 4.4a. y-a neutral curve for dendritic solidification of Pb-Sn with C„ - 0.1, r - 7.12 x 10~ 2 (G,- 50 K cm* 1 ) for 
several values of £2 0 . The solid, dashed-dot, and dashed curves in each figure are for £ 20 - 0, 300, 500 rpm, respectively. 
The stable combinations of y and a lie above the respective neutral curves. 
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Figure 4.4c. y-a neutral curve for dendritic solidification of Pb-Sn with C„ - 0.55, r - 7.12 x 10 -2 (G L - 50 K cm* 1 ) for 
several values of ty,. The solid, dashed-dot, and dashed curves in each figure are for fy, - 0, 300, 500 rpm, respectively. 
The stable combinations of y and a lie above the respective neutral curves. 
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Figure 4.6. For dendritic solidification of Pb-Sn with C„ - 0.55, r - 7.12 x 10" 2 (G L - 50 K cm- 1 ), and ft 0 - 500 rpm, 
amplitude of cBmensionless disturbances at critical conditions [y - 1.976 x 10~ 2 (V, - 61 p sec* 1 ) and a - 0.32] to 
c) composition in the liquid; d) volume fraction. 
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Figure B1. Liquid density versus T for C - 0 and 0.1. •: data of Mokrovskii & Regel (1952) for pure HgTe. o: data of Chandra & 
Holland (1983) for pure HgTe. 0 : data of Chandra & Holland for Hg 0 9 Cd 0 .iTe. — : least squares fit (3.7a,b) computed as per 
Appendix A. The low-temperature endpoint of each curve corresponds to the melting point determined from the phase diagram of 
Szofran & Lehoczky (1981). 
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